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This  dissertation  represents  an  effort  to  develop  a 
method  to  determine  the  effects  of  high  target  atom  directed 
velocity  and  temperature  on  interaction  rates  for  a flowing 
fissioning  gas.  Major  emphasis  is  placed  on  the  determination 
of  the  thermal  neutron  spectrum,  microscopic  neutron  cross 
section  behavior,  and  averaged  macroscopic  neutron  cross 
section  behavior. 

The  system  used  as  a model  in  this  study  is  a gas  cavity 
region  where  a low  density  mixture  of  fissioning  fuel  and 
working  fluid  is  in  thermodynamic  equilibrium  at  a high 
temperature  and  flowing  at  constant  velocity.  The  gas  cavity 
region  is  surrounded  by  a moderator  which  is  the  thermal 
neutron  source. 

A brief  review  of  gas  core  reactors  concepts  and  neutron 
thermal ization  is  presented.  The  development  of  these  topics 
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is  traced  historically,  and  the  scope  of  the  present  work  is 
established. 

The  exact  analytical  development  for  the  effective 
microscopic  cross  section  and  for  the  scattering  kernel  based 
on  the  works  of  Perkins  and  Oblow  is  presented. 

An  approximate  method  to  solve  the  problem  for  the 
conditions  of  completely  entrained  neutrons  and  uncollided 
neutrons  is  developed  and  is  tested  against  the  exact 
analytical  solution  for  extreme  conditions  of  temperature  and 
directed  velocity.  The  results  from  using  the  method  are  in 
agreement  with  the  analytical  solution.  The  method  uses  two 
different  velocity  frames  and  is  suitable  for  implementation 
in  a computer  code.  Results  for  the  effective  microscopic 
cross  section  calculated  using  this  approach  are  tested 
against  the  only  known  and  very  limited  results  which  are 
relevant  to  this  problem  and  show  good  agreement.  The 
developed  approach  is  implemented  in  the  computer  code  BRT-I, 
which  solves  the  integral  form  of  the  neutron  transport 
equation  in  the  thermal  energy  region,  following  its 
adaptation  to  run  in  a personal  computer.  The  modifications 
and  adaptations  included  in  BRT-I  led  to  the  creation  of  a new 
version  of  the  code  which  is  called  BRT-PCM. 

The  high  directed  velocity  and  temperature  effects  on  the 
thermal  neutron  spectrum,  on  the  averaged  macroscopic  cross 
sections,  and  on  the  average  inverse  neutron  velocity  are 
established  and  explained. 
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CHAPTER  1 
INTRODUCTION 

-It  1 Brief  Review  of  Gas  Core  Reactor  Studies 

The  first  idea  that  motivated  this  research,  was  the  need 
to  determine  the  reaction  rates  that  occur  throughout  a 
supersonic  magnetohydrodynamic  duct  in  a gas  core  reactor 
power  system  that  employs  a closed-cycle  disk 
magnetohydrodynamic  (MHD)  generator  for  energy  conversion.  In 
this  particular  system  [1],  the  fissionable  gas  flows  out  the 
core,  passes  through  a nozzle,  and  enters  a disk  MHD  generator 
that  is  configured  to  be  an  integral  part  of  the  reactor. 
Thus,  significant  fissioning  occurs  throughout  the  MHD  duct 
and  fission  fragment-induced  ionization  of  the  uranium-bearing 
fuel  gas/working  fluid  is  anticipated  to  yield  the  required 
nonequilibrium  electrical  conductivity  despite  the  relatively 
low  (for  an  MHD  device)  gas  temperature  [2].  Since  this  and 
other  previewed  applications  involve  gas  core  reactor 
concepts,  a brief  historical  review  of  gas  core  reactor 
development  is  presented. 

Gas  core  reactors  use  a gaseous  fissile  material  as  fuel 
for  power  generation.  The  fissioning  gas  is  inside  the  reactor 
core.  A wide  variety  of  physical  arrangements,  employing  ex- 
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core  or  in-core  moderation  of  neutrons  and  methods  of 
containing  the  gas  and  extracting  the  energy,  has  been 
suggested. 

One  of  the  first  reported  works  on  the  use  of  fissionable 
gases  in  nuclear  reactors  was  made  by  George  Bell  [3]  at  Los 
Alamos  in  1955.  His  reports  were  favorable  and  this  motivated 
more  work  on  the  same  concept.  These  works  intended  to 
determine,  fundamentally,  the  feasibility  of  using  gas 
as  a reactor  fuel,  because  this  gas  in  enriched  form,  was 
readily  available  from  gaseous  diffusion  plants. 

The  climax  of  this  early  work  happened  when  a Russian 
experimental  group,  headed  by  Kikoin  [4],  reported  that  it  had 
assembled  a critical  configuration  using  gas  as  a reactor 
fuel  in  1959. 

The  gas  core  reactor  concept  was  more  seriously 
considered  in  the  late  1950s,  with  the  newborn  space  age.  The 
gas  core  which  could  reach,  at  least  theoretically,  ultra-high 
temperatures  with  an  enormous  amount  of  heat  transfer 
capability  suggested  a highly  viable  alternative  for  space 
rocket  propulsion.  With  this  intent  in  mind,  many  studies  were 
performed  on  gaseous  core  reactor  concepts.  Gas  core 
propulsion  reactor  concepts  studied  in  the  '60s,  which 
promised  superior  heat  transfer  and  fluid  flow 
*^haracter ist ics , were  the  nuclear  light  bulb  [5]  , the  vortex 
[6],  and  the  coaxial  [5,7]  flow  cavity  reactors. 
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Of  the  many  reactor  concepts  proposed,  the  majority  was 
cavity  core  reactors,  characterized,  basically,  by  a large 
spherical  cavity  fueled  with  uranium  hexafluoride  gas  or 
uranium  metal  vapor.  In  the  propulsion  concepts,  the  vaporized 
uranium  fuel  stream  is  surrounded  by  a lighter,  faster  moving 
hydrogen  propellant  stream.  The  core  is  then  contained  inside 
a beryllium,  beryllium  oxide,  heavy  water  or  graphite 
moderating-reflector  shell.  Neutrons  born  in  the  cavity  as  a 
result  of  fission  go  into  the  moderator,  where  they  are 
thermalized,  and  then  diffuse  back  into  the  cavity  for  further 
fissions  characterizing  a thermal,  heterogeneous  reactor 
[5,8]  . 

Another  conceptual  application  of  gas  core  reactors  that 
appeared  in  the  middle  60s  was  their  use  in  direct 
magnetohydrodynamic  energy  conversion  systems,  not  only  for 
space  electric  propulsion,  but  also  for  commercial  power 
production  [9].  The  space  application  was  also  emphasized 
later  as  one  of  the  most  promising  ways  to  use  closed-cycle 
MHD. 

Some  other  new  concepts  [9]  were  also  developed  in  the 
70s,  like  the  Mixed-Flow  Gas  Core  Reactor  (MFGCR)  [10].  The 
main  feature  of  this  system  is  the  use  of  clusters  of  seven 
externally  Be  moderated  cavities,  surrounded  by  a molten  salt 
breeding  blanket  in  which  the  fuel  gas  mixture  passes  through 
the  cores  and  transports  the  heat  to  the  intermediate  heat 
exchangers.  The  use  of  multiple  cavities  leads  to  an  enhanced 
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moderation  as  compared  to  the  single  cavity  concept.  The  word 
"mixed"  stands  for  the  mixing,  at  the  outlet  of  each  cavity, 
of  the  hot  portion  of  the  gas  from  the  central  part  of  the 
cavity  with  the  cooler  flow  of  gas  that  cools  the  cavity  wall. 

Different  and  interesting  concepts,  also  from  the  70s, 
are  related  to  pulsed  systems,  which  have  the  advantage  of 
attaining  high  peak  values  of  variables  like  temperature  and 
pressure,  while  maintaining  relatively  low  cycle-averaged 
values.  Examples  of  these  systems  are  the  Pulsed  Nuclear 
Piston  (PNP)  [11,12,13],  known  as  the  nuclear  piston  engine, 
and  the  Pulsed  Gas  Generator  (PGG)  [14,15].  The  PNP  consists 
of  a small  pulsed  gaseous  core  reactor  with  variable  volume 
enclosed  by  a moderating-reflector  cylinder  and  piston 
assembly  which  operates  on  a thermodynamic  cycle  similar  to 
the  internal  combustion  engine.  The  PGG  is  like  the  PNP  but 
uses  a core  of  fixed  dimensions;  thus,  it  has  the  advantage  of 
mechanical  simplicity,  since  it  does  not  require  a moving 
piston.  Yet  it  also  has  lower  efficiency  due  to  the  lack  of 
mechanical  power.  The  PNP  and  PGG  typically  achieve  a few 
megawatts  of  power  generation  per  chamber. 

A more  recent  concept  [1,2]  developed  in  the  late  1980s 
that  combines  the  use  of  gas  core  reactors  and  closed-cycle 
MHD  generator  conversion  to  supply  energy  in  space  is  the 
innovative  Ultra  High  Temperature  Vapor  Core  Reactor  (UTVR) / 
Disk  MHD  generator  system.  The  UTVR/MHD  generator  does  not  use 
or  uranium  metal  vapor,  but  UF^,  which  has  extremely 
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desirable  saturation  vapor  pressure-temperature  behavior, 
insofar  as  direct  Rankine  cycle  operation  for  space  power  is 
concerned. 

The  gaseous  fuels  considered  earlier,  UF^  and  uranium 
metal  vapor,  present  problems  if  used  in  a direct  Rankine 
cycle  for  space  power.  The  UF^  is  limited  to  a Brayton  type 
cycle  when  employing  gas  turbines  or  MHD  generators,  because 
the  low  condensing  temperature  of  UF^  required  in  a direct 
Rankine  cycle  (approximately  400°K)  is  not  acceptable  for  heat 
rejection  in  space.  Uranium  metal  vapor  fueled  systems  also 
present  problems;  although  easily  liquified  at  the  heat 
rejection  or  condensing  end  of  the  Rankine  cycle  at  the 
pressures  required  for  criticality,  very  high  gas  core 
temperatures  (everywhere  above  6000°K)  are  needed  to  prevent 
condensation. 


1.2  Objectives 

The  objectives  of  this  research  are  to  determine  the 
effects  of  the  high  target  atom  directed  velocity  and 
temperature  on  the  thermal  neutron  spectrum,  effective 
microscopic  neutron  cross  sections,  and  averaged  macroscopic 
neutron  cross  sections  and  also  to  determine  how  these  effects 
change  when  selected  parameters  like  temperature  and  directed 
velocity  are  changed.  Another  objective  is  to  develop 
simplified  approximate  analytical  expressions  which  can  be 
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used  to  determine  the  combined  effects  of  high  target  atom 
directed  velocities  and  temperatures  on  interaction  rates  for 
a high  velocity,  flowing,  fissioning  gas.  The  completely 
entrained  and  uncollided  neutron  conditions  are  studied. 

By  high  directed  velocities  and  temperatures,  we  mean 
Mach  numbers  greater  than  one,  and  at  least  a few  thousands  of 
degrees  Kelvin.  The  way  in  which  very  high  temperatures  affect 
cross  sections  was  studied  before  [5],  but  the  combined 
effects  of  very  high  temperature  and  directed  velocity  of 
target  atoms  on  neutron  interaction  rates  are  a new  subject. 

The  usual  approach  of  trying  to  find  effective  cross 
sections  will  be  developed  in  a way  to  permit,  wherever 
possible,  the  use  of  these  effective  cross  sections  (which 
will  include  the  effects  of  both  directed  velocity  and 
temperature)  with  conventional  treatments,  to  obtain  correct 
reaction  rates. 

The  effects  of  very  high  fuel  gas  temperatures  have 
already  been  studied  by  Oblow  [5],  but  his  analysis  did  not 
include  the  effects  of  high  directed  velocity  target  atoms 
because  he  modeled  a cavity  reactor  where  only  thermal  motion 
was  considered.  Oblow 's  analysis  established,  for  selected 
fuel/moderator  systems,  the  regimes  where  the  effects  of  high 
fuel  gas  temperatures  and  the  resultant  target  atom 
distributed  velocity  become  significant  and  cannot  be 
adeguately  treated  by  conventional  algorithms. 
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This  research  is  intended  to  be  an  extension  of  Oblow's 
work,  but,  in  this  case,  the  new  feature  of  very  high  directed 
velocity  is  added  to  very  high  temperature. 

The  first  part  of  this  work  is  an  effort  to  determine 
expressions  for  the  effective  cross  sections  and  scattering 
kernel,  using  a rigorous  treatment.  Most  of  the  theoretical 
expressions  that  will  result  from  this  work,  when  the  rigorous 
treatment  is  applied,  will  not  be  analytically  solvable. 
Hence,  the  next  step  will  be  to  develop  approximations  to  the 
problem  in  order  to  make  possible  the  use  of  a pre-existing 
computer  code  to  solve  the  problem  when  both  the  directed 
velocity  of  the  target  atoms  and  high  temperatures  are 
present.  The  computer  code  will  have  to  be  modified  to  include 
the  developed  approximations  based  on  the  rigorous  analytical 
expressions  for  the  effective  cross  sections. 

Using  the  modified  computer  code,  it  is  proposed  to 
establish  regimes  where  high  target  atom  directed  flow 
velocity  effects  become  significant;  to  determine  the 
boundaries  where  high  directed  flow  velocity  effects  begin  to 
dominate  over  temperature  effects,  as  a function  of 
temperature;  and  to  determine  where  temperature  effects  begin 
to  dominate  over  directed  flow  velocity  effects,  as  a function 
of  flow  velocity. 

The  effective  cross  sections  expressions  to  be  developed 
will  allow  the  target  atom  total  velocity  distribution  to  be 
considered,  including  both  the  Maxwell-Boltzmann  thermal 
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velocity  component  and  the  directed  flow  velocity  component, 
and  used  in  the  relations  for  the  interaction  rates. 

The  effects  of  high  target  atom  directed  velocity  and 
temperature  will  be  studied  as  a function  of  gas  composition 
(the  only  fissionable  species  considered  are  U235  and  U23g)  , gas 
velocity,  and  temperature  by  using  the  developed  computer 
code. 

It  is  assumed,  in  general,  that  we  are  working  with  an 
homogeneous  medium  composed  of  an  ideal  gas  in  thermodynamic 
equilibrium  at  a specific  temperature,  with  a Maxwellian 
distribution.  In  addition  it  is  assumed  that  all  neutrons 
entering  the  system  are  thermalized  elsewhere  and  no  sources 
of  thermal  neutrons  are  present  within  the  system. 

The  system  modeled  is  a gas  cavity  region  where  a low 
density  mixture  of  fissioning  fuel  and  working  fuel  is  in 
equilibrium  at  a high  temperature  and  flowing  at  constant 
velocity.  The  gas  cavity  region  is  surrounded  by  a moderator 
which  is  the  thermal  neutron  source.  The  choice  of  this  system 
to  serve  as  a model  is  based  on  the  magnetohydrodynamic 
section  of  the  UTVR/Disk  MHD  generator  system  which  was  the 
original  motivation  for  this  work,  but  this  choice  does  not 
imply  that  the  present  effort  is  in  any  way  limited  to  the 
materials  and  densities  in  the  magnetohydrodynamic  section  of 
the  UTVR/Disk  MHD  generator  system. 

Considering  that  the  inclusion  of  the  directed  velocity 
in  the  effective  cross  section  can  make  evaluation  of  even  a 
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simplified  or  approximate  effective  cross  section  expression 
rather  cumbersome,  an  effort  will  be  made  to  develop  a simple 
method  that  can  be  used  to  quickly  determine  whether  or  not  it 
is  necessary  to  include  the  directed  velocity  in  the 
calculations. 


1.3  Overview 


Almost  all  the  work  done  in  the  past,  although 
recognizing  the  dependence  of  the  cross  sections  on  the 
relative  neutron-target  velocity,  usually  considered  the 
velocity  of  the  target  to  be  much  smaller  than  the  velocity  of 
the  neutron,  and  simply  ignored  the  former  in  the  calculation 
of  cross  sections  in  the  slowing  down  region. 

The  velocity  of  the  target  due  to  thermal  motion  is 
relevant  in  the  thermal  region,  whose  upper  limit  is  usually 
considered  to  be  around  1 eV  [16],  and  is  included  in  the 
calculations  by  means  of  relative  velocities.  Also,  in  the  low 
energy  resonance  region  the  target  atom  velocity  can  have 
substantial  importance  because  it  gives  rise  to  a Doppler 
shift  effect.  The  treatment  normally  employed  in  this  case  is 
to  use  the  relative  velocity  in  the  Breit-Wigner  [17]  single 
level  resonance  formula  (which  uses  resonance  parameters  to 
determine  resonance  cross  sections  as  a function  of  energy) , 
and  the  Bethe-Placzek  [16]  approximation  for  evaluation  of 
averaged  resonance  cross  sections  when  temperature  effects  are 


10 


taken  into  account.  This  can  be  considered  a reasonable 
approximation  for  almost  all  cases  in  practice,  when 
considering  current  commercial  power  reactors,  but  doesn't 
apply  for  the  particular  case  of  a high  directed  velocity, 
high  temperature,  flowing  fissioning  gas  because  the  Bethe- 
Placzek  approximation  does  not  give  reasonable  results  for 
high  target  atom  velocities  and  low  energy  resonances.  For  the 
case  of  a high  directed  velocity,  high  temperature,  flowing, 
fissioning  gas,  even  the  treatment  for  the  thermal  energy 
range,  where  the  target  atoms  are  assumed  to  have  a Maxwell- 
Boltzmann  type  distribution,  will  require  significant 
modification. 

In  most  applications  of  nuclear  energy,  parts  of  a system 
are  in  translational  motion  with  respect  to  each  other.  This 
motion  can  introduce  variations  in  the  behavior  of  the  system, 
affecting  its  neutronics  parameters.  For  example,  material 
speeds  as  low  as  3 m/s  in  DjO  and  25  m/s  in  H2O  increase  the 
thermal  diffusion  length  by  10%  over  that  of  the  stationary 
case  [18].  Although  this  effect  is  usually  neglected  in 
neutronics  analysis,  under  certain  conditions  it  is  not 
negligible  [18].  This  is  certainly  the  case  for  high  target 
atom  directed  velocity. 

The  high  target  atom  directed  velocity  and  very  high 
temperature  of  the  flowing  fissionable  gas  can  cause  the 
thermal  neutron  spectrum  to  be  shifted  to  higher  energies, 
when  compared  with  usual  energy  spectra  found  in  conventional 
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reactor  systems  where  temperatures  and  velocities  are  usually 
much  smaller. 

To  study  the  effects  of  high  target  atom  directed 
velocity  and  temperature  this  research  is  strongly  based  on 
the  classical  theory  of  neutron  thermal ization,  related  to  an 
ideal  gas  model. 


1.4  Historical  Review  of  Neutron 
Slowing  Down  and  Thermalization 

As  the  present  work  is  focused  primarily  on  the  neutron 
thermalization  problem  under  special  circumstances,  a brief 
historical  review  of  neutron  thermalization  is  presented  here. 

In  the  year  1930,  Bothe  and  Becker  discovered  artificial 
nuclear  gamma-radiation  by  bombarding  Be  with  alpha-particles. 
The  reaction  Be+alpha  was  subsequently  investigated  by  Joliot 
and  Curie.  From  the  results  of  those  experiments,  Chadwick 
drew  the  conclusion  that  besides  the  gamma-radiation  already 
found,  there  was  another  neutral  particle  produced,  with  mass 
approximately  equal  to  that  of  a proton,  which  he  called  a 
neutron.  According  to  this,  the  history  of  the  neutron  can  be 
considered  as  having  as  a starting  point  its  discovery  by 
James  Chadwick  in  1932  and  was  radically  altered  after  the 
discovery  of  fission  by  Hahn  and  Straussmann  in  1939. 

The  subject  of  neutron  thermalization  and  slowing  down 
can  be  considered  as  new,  when  compared  with  modern  physics  as 
a whole,  but,  the  mathematics  involved  is  very  old  since  the 
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treatment  of  particle  diffusion  and  dynamical  behavior  under 
various  scattering  laws  has  remained  unchanged,  since  its 
development  by  Maxwell  and  Boltzmann  (1859,1895).  Actually, 
Boltzmann's  original  equation  can  describe  very  accurately  the 
transport  of  neutrons,  and  due  to  the  very  low  density  of  the 
"neutron  gas"  in  matter,  the  nonlinear  equation  can  be 
simplified,  allowing,  in  a few  special  cases,  analytical 
solutions  [19], 

Fermi  et  al.  [19],  in  1934,  experimentally,  studying 
neutron  diffusion  in  paraffin  wax,  proved  that  neutrons  can 
lose  a large  quantity  of  energy  when  colliding  with  hydrogen 
nuclei.  Fermi  concluded  that  provided  the  neutron-proton 
elastic  scattering  cross  section  is  much  greater  than  the 
corresponding  absorption  cross  section,  collisions  against 
protons  would  reduce  the  neutron  energy  to  a value  much 
smaller  than  the  initial  one;  he  repeated  his  experiments 
using  other  materials  without  hydrogen  and  the  results 
indicated  that  neutron  slowing  down  was  occurring,  but  slower 
than  in  the  hydrogenous  media.  Despite  his  early  work  in 
slowing  down,  Fermi  could  not  show,  experimentally,  the 
effects  of  temperature  on  the  neutron  slowing  down  process. 
The  experimental  proof  of  the  so  called  "thermal  effect"  was 
presented  by  Moon  and  Tillman  in  1935  [19].  They  showed  that 
the  ratio  of  the  activity  of  a silver  cylinder  in  paraffin  wax 
at  90°K  to  the  activity  when  the  paraffin  was  at  room 
temperature  was  approximately  1.26. 
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The  first  calculations  involving  neutron  slowing  down 
were  made  by  Fermi  [19]  who  presented  the  results  in  the 
"Ricerca  Scientifica"  in  1936.  In  this  paper,  that  even  today 
cannot  be  considered  outdated,  Fermi  deduced  the  probability 
of  energy  loss  during  neutron  slowing  down,  and  formulated  an 
infinite  medium  balance  equation  for  the  neutron  energy 
spectrum.  He  established  the  basis  for  solving  many  related 
problems,  like  the  calculation  of  slowing  down  length  by 
probability  methods,  the  diffusion  and  leakage  of  neutrons 
into  a vacuum,  and  the  meaning  of  albedo.  He  also  proved  that 
for  hydrogen,  the  thermal  absorption  varies  with  the  inverse 
of  the  neutron's  velocity  and  showed  that  the  spectrum  of 
epithermal  neutrons  varies  inversely  with  energy. 

Ornstein  and  Uhlenbeck  appear  to  have  been  the  first  to 
make  use  of  the  linearized  Boltzmann  transport  equation  to 
solve  slowing  down  problems,  like  the  slowing  down  length  of 
neutrons  from  a point  source  in  an  infinite  medium,  in  1937 
[19].  Their  results  agreed  very  well  with  Fermi's  results 
which  had  been  derived  from  quite  different  considerations. 

Although  the  thermal ization  problem  was  first  considered 
by  Fermi  [19]  with  his  scattering  laws  in  193  6,  the  first 
significant  step  in  this  development  is  historically 
considered  to  be  the  paper  of  Wigner  and  Wilkins  [20]  in  1944. 
Before  this,  Arley  [19]  included  the  effect  of  thermal  motion 
on  the  scattering  by  hydrogenous  molecules  and  neutron 
interactions  with  an  isotropic  oscillator,  Sachs  and  Teller 
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[19]  treated  the  scattering  of  slow  neutrons  by  molecular 
gases,  treating  the  gases  as  classical  rotators,  and  many 
other  authors,  like  Pomeranchuk,  and  Weinstock  [19]  developed 
and/or  improved  methods  to  treat  the  problems  involved  in 
neutron  thermalization;  however,  despite  the  importance  of 
their  investigations,  they  were  all  qualitative. 

Actually,  the  neutron  thermalization  problem  is 
essentially  the  solution  of  Boltzmann  equation,  subject  to 
adequate  boundary  conditions,  and  coupled  with  the  scattering 
kernel.  Even  today,  the  solution  of  this  problem  is  a 
formidable  one. 

Wigner  and  Wilkins  [20],  considering  an  infinite,  perfect 
monoatomic  gaseous  medium  containing  an  inverse  velocity 
absorber  and  a constant  high  energy  spatially  distributed 
reduced  the  Boltzmann  equation  from  seven  independent 
variables  to  only  one,  namely,  energy.  By  making  some  other 
assumptions,  they  were  able  to  solve  the  equation  and,  most 
important,  from  the  solution  they  could  state  the  physical 
fact  that  the  rate  of  thermalization  depends  only  on  the  ratio 
of  the  2200  m/s  absorption  cross  section  to  the  free  atom 
scattering  cross  section  at  ambient  temperature.  The 
preferential  absorption  at  low  energies  explained  the 
distortion  of  the  neutron  spectrum,  from  the  expected 
Maxwellian,  in  this  region. 

The  extension  of  this  model  to  gases  of  higher  mass, 
that  would  allow  the  analysis  of  solids  and  liquids  by  means 
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of  an  "effective  mass",  was  also  obtained  by  Wilkins  [21]  in 
1944.  The  so  called  "heavy  gas  model"  permitted  the  Boltzmann 
equation  to  be  solved  exactly  for  other  variables,  like  time 
and  space  [19,22] . 

The  first  practical  solution  of  a neutron  thermal ization 
problem  was  obtained  by  Brown  and  St.  John  [23],  in  1954. 
Using  an  effective  mass  given  by  Sachs  and  Teller,  they  solved 
the  exact  gas  integral,  which  is  the  time  and  space 
independent  Boltzmann  equation,  using  a monoatomic  gas  model. 
The  equation  was  solved  numerically,  to  be  applied  in  the 
calculation  of  the  neutron  spectrum  for  the  Savannah  River, 
DjO-moderated,  reactor.  They  also  tabulated  the  effective 
fission  and  absorption  cross  section  of  plutonium  that 
presents,  as  an  interesting  aspect,  a strong  resonance  at  0.3 
eV.  The  principal  motivation,  at  that  time,  for  the  research 
on  neutron  slowing  down  and  thermalization  was  the  accurate 
calculation  of  cross  sections  that  would  determine  the 
behavior  of  neutrons  in  the  system. 

Cohen  [24],  in  1955,  presented  a method  to  determine  the 
upper  limit  of  the  thermal  energy  group  and  made  calculations 
on  the  modification  of  neutron  age  by  thermalization  effects. 

Von  Dardel  [19]  discovered  the  effect  of  "diffusion 
cooling,"  that  brought  the  related  problem  of  "diffusion 
heating",  and  developed  a theory  about  it  based  on  an  energy 
balance  equation.  "Diffusion  cooling"  is  the  softening  of  the 
neutron  spectrum  caused  by  the  preferential  leakage  of  high 
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energy  neutrons  from  regions  where  the  spatial  variation  of 
the  flux  has  a positive  curvature.  "Diffusion  heating"  is 
another  effect  that  causes  the  neutron  flux  to  deviate  from 
Maxwellian  behavior  because  neutrons  of  higher  energies 
diffuse  into  a given  volume  element  more  readily  than  do  those 
of  lower  energies.  The  first  effect  was  also  studied  by 
Beckurts  [25],  Hurwitz  and  Nelkin  [22],  Nelkin  [19], 
Kazarnovsky  [19],  and  others,  until  1960,  when  Nelkin  [26] 
presented  an  exact  transport  theory  expression  for  the 
diffusion  cooling  coefficient  which  showed  how  the  diffusion 
cooling  coefficient  is  connected  with  that  part  of  the  neutron 
spectrum  which  represents  the  deviation  from  equilibrium 
resulting  from  leakage  effects. 

The  first  remarkable  try  to  solve  the  problem  of  the 
spatial  dependence  of  the  neutron  flux  near  boundaries  was 
made  by  Conkie  [27],  in  1960.  His  method  was  based  on  a 
approximation  to  the  Boltzmann  equation,  employing  the 
hydrogen  gas  scattering  kernel,  and  showed  the  type  of 
convergence  that  could  be  expected  from  the  polynomial 
representation.  In  the  same  year,  Honeck  [19]  and  Honeck  and 
Kaplan  [19]  used  the  integral  equation  formulation  of  the 
Boltzmann  equation  to  calculate,  numerically,  the  fine 
structure  flux  in  a cylindrical  water  cell  lattice,  using  a 
gas  model  with  an  effective  mass,  and  their  results  can  be 
considered  reasonable  when  compared  to  experiments.  These 
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calculations  formed  the  basis  for  the  THERMOS  code  developed 
by  Honeck  [ 28 ] . 

Corngold  [29]  solved  the  problem  of  slowing  down  in  a 
space  independent  infinite  medium,  and  obtained  an  exact,  but 
asymptotic,  solution  for  the  distribution  of  neutrons  slowing 
down  from  a high  energy  source,  which  is  valid  down  to  about 
0.2  eV. 

Corngold,  Michael  and  Wollman  [30]  showed  that,  in  the 
case  of  hydrogen  gas,  the  time  eigenvalues  have  a limit  point 
plus  a continuum  region. 

In  1966,  Blackshaw  and  Murray  [31]  investigated,  within 
a classical  physics  framework,  the  elastic  scattering  of  low- 
energy  neutrons  by  the  nuclei  of  a monoatomic  gas,  which  have 
an  isotropic  Maxwellian  velocity  distribution.  Using  a unified 
nis^themat  ical  treatment,  and  also  unigue  single— integral 
expressions  valid  for  any  analytical  or  numerical 
representation  of  the  microscopic  scattering  and  absorption 
cross  sections  as  functions  of  the  relative  neutron— nuclear 
speed,  they  developed  a new  form  of  the  scattering  probability 
function  in  velocity  space;  the  velocity  scattering  kernel  and 
its  moments;  and  total  scattering  and  absorption 
probabilities.  This  permitted  the  exploration  of  a broader 
class  of  low-energy  scattering  problems. 

Many  other  authors  have  presented  relevant  works  in  this 
area.  Unfortunately,  it  is  not  possible  to  name  all  of  them  in 
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this  brief  review;  this,  however,  does  not  mean  that  their 
contribution  is  less  than  the  ones  presented. 


CHAPTER  2 

THEORETICAL  DEVELOPMENT 


2 . 1 Introduction 


The  particular  system  from  which  the  idea  for  this  work 
arose  is  presented  in  Figure  2.1.  The  figure  represents  a 
section  of  the  UTVR/MHD  Generator  Rankine  Cycle  Power  System 
where  the  gaseous  core,  the  moderator,  the  nozzle  and  the  MHD 
section  can  be  identified. 

In  this  system,  the  velocity  of  the  target  atoms  in  the 
gaseous  core  is  relatively  small  (approximately  50  m/s)  and 
can  be  neglected,  as  is  usually  done  in  commercial  power 
reactors  in  the  neutronics  calculations.  However,  after 
passing  through  the  nozzle,  the  velocity  of  the  target  atoms 
is  substantially  increased,  reaching  Mach  numbers  of  around 
three,  and  the  mean  gas  density  is  significantly  decreased  in 
the  MHD  duct  section  relative  to  the  density  in  the  gas  core 
(a  factor  of  approximately  20) . 

The  specific  subject  that  motivated  this  research  is  the 
desire  to  determine  as  accurately  as  possible  the  interaction 
rates  in  the  region  inside  the  MHD  duct,  which  is  expected  to 
have  a combination  of  high  target  atom  velocity  and 
temperature.  The  reason  for  focusing  on  this  particular  region 
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DUCT 


Figure  2.1  - Representation  of  UTVR/DISK  MHD. 
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is  that  the  concept  depends  upon  being  able  to  obtain  high 
electrical  conductivity  through  fission— induced  noneguilibrium 
ionization;  knowledge  of  the  interaction  rates  in  the  MHD  duct 
is  therefore  critical.  It  should  be  noticed  that  the  region 
inside  the  nozzle  is  also  expected  to  have  a combination  of 
high  target  atom  directed  velocity  and  temperature  and  a 
higher  mean  gas  density  than  the  MHD  duct,  but  since 
fissioning  in  this  region  does  not  contribute  to  enhanced 
electrical  conductivity  in  the  MHD  duct,  accurate  knowledge  of 
the  fission  (interaction)  rate  in  this  region  is  not  crucial 
to  concept  viability. 

Although  it  is  anticipated  that  the  target  atom  will 
present  a velocity  profile  inside  the  MHD  duct,  which  is 
certainly  not  flat,  a constant  velocity  profile  for  the 
medium,  representing  an  average  velocity,  is  assumed  in  this 
work. 

In  the  UTVR/Disk  MHD  generator  concept,  the  radial 
dimension  of  the  duct  is  comparable  to  the  core  dimensions 
(approximately  1 m)  while  the  duct  height  (approximately  0.2 
m)  is  reasonably  smaller  than  the  core  dimensions  and  the  mean 
gas  density  in  the  MHD  generator  is  typically  about  5%  of  the 
density  inside  the  core  [1].  Conseguently , essentially  all  the 
neutrons  generated  by  fission  inside  the  duct  escape  to  the 
moderator  without  undergoing  any  interaction  in  the  duct  while 
they  are  fast;  the  neutrons  are  then  thermalized  in  the 
moderator.  Thus,  it  is  assumed  that  all  the  neutrons  entering 
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the  duct  are  thermalized  and  that  there  are  no  sources  of 
thermal  neutrons  in  the  gas  medium.  In  addition,  it  is  assumed 
that  the  gas  medium  is  homogeneous,  in  thermal  equilibrium  at 
a specific  temperature  and  at  steady  state. 

The  purpose  of  this  analytical  work  is  to  determine  an 
sppi^opriate  formulation  for  the  collision  frequency  when 
studying  neutron  interactions  with  target  atoms  subjected  to 
directed  velocity  in  a medium  at  very  high  temperature, 
as  in  the  MHD  duct  described  above  and  then,  using  the 
obtained  formulation,  to  determine  the  effective  neutron  cross 
sections  and  the  scattering  kernel  in  the  thermal  energy 
region. 


2.2  Collision  Frecfuency 

A good  starting  point  is  the  derivation  of  the 
expressions  that  represent  the  collision  frequency  of  neutrons 
with  atoms  having  a high  directed  velocity  in  a high 
temperature  medium.  An  homogeneous  medium  composed  of  a ideal 
gas  in  thermodynamic  equilibrium  at  a specific  temperature 
will  be  assumed. 

In  a stationary  medium,  where  the  target  atom  velocity, 
V,  is  equal  to  zero,  the  collision  frequency  F(v) , that  is, 
the  frequency  with  which  reactions  occur,  for  a single  neutron 
moving  with  velocity  v is  defined  by  the  expression: 

F(v)=vE*(v)  . 


(2-1) 


Considering  that 
S‘(v)=Nct^(v) 
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(2-2) 

the  collision  frequency  can  be  expressed  as 

F(v)=vNct‘(v)  . (2-3) 

In  the  case  of  a target  moving  with  a velocity 
distribution,  D(V) , the  cross  sections  are  a function  of  the 
relative  velocity,  v^,  between  the  neutron  and  the  target, 
such  that, 

Vp=v-V.  (2-4) 

Then,  in  this  case,  the  collision  frequency  needs  to  take 
into  account  the  target  atom  velocity  distribution. 

Considering  only  the  thermal  motion,  and  assuming  that 
the  target  atom  velocity  distribution  is  described  by  the 
Maxwell-Boltzmann  distribution  characterizing  an  ideal  gas  in 
thermal  equilibrium  at  a temperature  T,  the  target  atom 
velocity,  V,  is  now  equal  to  V^,  the  velocity  due  to  the 
thermal  motion,  and  the  velocity  distribution  is  as  follows: 


D (V)  =D^  (V^)  = (m/27rkT)  ^''^exp  (-mV^V2kT)  or, 
D(V)=D,(V,)  = (BV7r2^^)exp(-B2v,2)  , where 
B=(m/2kT) 

The  collision  frequency  is  then  formulated  as 
F(v,V^)=J  |v-vj ‘a^dv-vj  ) •N»D(V^)dV^,  and 


(2-5) 

(2-6) 

(2-7) 


Vp=v-V,. 


(2-8a) 
(2-8b) 

If  the  target  atom  is  now  considered  as  moving  at  a specific 
directed  velocity,  without  thermal  motion,  the  velocity  of  the 
target  atom  will  be  V^,  such  that: 
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V=V„,  and  (2-9) 

v=v-V„.  (2-10) 

In  this  case,  the  collision  frequency  can  be  written  as 

F(v,V„)  = |vjNaMvJ  ) , or  (2-11) 

F(v,V„)=  v^Nct^(v^)  . (2-12) 

The  above  two  collision  frequency  formulations  take  into 
account,  separately,  the  thermal  motion  and  the  directed 
velocity  of  the  target  atom.  Now  consider  a system  where  the 
target  atoms  are  subjected  simultaneously  to  the  combined 

effects  of  thermal  motion  and  directed  velocity.  The  geometry 
that  represents  this  system  is  shown  in  Figure  2.2,  where  v 
and  V are  the  neutron  and  target  velocities  in  the  laboratory 
system,  related  to  a static  frame  and  is  the  relative 
neutron-nuclear  velocity.  Also  in  this  figure,  V„  is  the 

directed  velocity  of  the  medium,  is  the  thermal  velocity 
and  is  the  relative  velocity  between  the  neutron  and  the 
directed  velocity  of  the  medium. 

If,  as  in  Figure  2.2,  the  medium,  as  a whole,  moves  with 
directed  velocity  V„,  a fixed  external  observer  will  see  that 
the  target  velocity,  V,  is 


v=v,+v„ 


(2-13a) 


and  the  neutron-nucleus  relative  velocity  in  this  case  is 


Vr=v-(V^+V„)  . (2-13b) 

The  collision  frequency  for  the  case  where  the  medium  is 
moving  with  directed  velocity  and  the  target  atoms  are  also 
submitted  to  thermal  motion  can  be  described  by  an  equation 
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LEGEND: 

v^:  Neutron-target  relative  velocity  (cm/s) 

Vjyj:  Target  atom  directed  velocity  (cm/s) 

Vjyj:  Relative  velocity  between  neutron  and 

target  atom  directed  velocity  component  (cm/s) 
Vrp:  Target  atom  thermal  velocity  (cm/s) 

V : Neutron  velocity  (cm/s) 

V : Target  atom  total  velocity  (cm/s) 

0 : Polar  angle  of  Vj^  (rad) 

: Azimuthal  angle  of  Vj^  (rad) 

0q:  Polar  angle  of  Vjyj  related  to  (rad) 


Figure  2.2  - Geometry  of  Neutron-nucleus  Interaction. 
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similar  to  Equation  (2-8)  , by  replacing  with  V,  or,  in  terms 
of  relative  velocity, 

F(V,V)=  / v^»a^(v^)«N»D(V)dV.  (2-14) 

In  the  above  expression,  D(V)dV  is  the  probability  that 
the  target  resultant  velocity  lies  within  dV,  at  V. 

Equation  (2-14)  can  be  changed  to  the  domain  of  the 
relative  velocity  to  obtain: 

F(V,V)=  J D(V(v^)  ) v^CT^(v^)Ndv^.  (2-15) 

S ince  dv^=v^2^j^^jj  ^ ^ 
the  above  can  be  written  as 


F(v,  V)  =//  D(V (v^)  ) Vj.»  (Vj.)  ‘N*  Vj,^dVj,d(n) 
where. 


(2-16) 


D(V(vJ)=v^2j  D(V(vJ)d(n)^^.  (2-17) 

Then,  substituting, 

F(v,V)=J  D(V(v^)  )v^a^(v^)  »Ndv^.  (2-18) 

Thus,  it  now  remains  to  determine  the  distribution  function 
D(V(v^)).  A similar  problem  was  studied  by  Perkins  [18]  and 
his  solution  was  modified  for  the  present  case.  The  obtained 
results  are  presented  below: 

The  velocity  distribution,  D(V(v^)),  can  be  written  as 
D(V(vJ  )dv^=/D„(v„'  )dv„'D^(v„'  ,v^)dv^  (2-19) 

where  the  distributions  are  defined  as: 

DH(v„')dv„':  Probability  that  the  directed  velocity 

component  V„ ' yields  v„  ' within  dv„  ' . 
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Probability,  for  a fixed  v„',  that  the 
temperature  velocity  component,  V^,  yields  within  dv^  which 
in  the  case  of  a Maxwellian  distribution  can  be  expressed  as 
/Vp)  = [BV7r^^^]exp[-B2(v^2+^„'2-2v^v„'/X(j,)  . (2-20) 

case  of  constant  speed,  monodirectional  motion 
for  the  directed  velocity,  has  a speed  in  the  direction 
0,  0 (see  Fig  2.2).  This  corresponds  to  a relative  speed  v„ 
between  the  neutron  and  the  medium  in  the  direction  0^  with  an 
azimuthal  angle  a (both  relative  to  v ) . Then,  D (v  ' ) dv  ' is 
(v„ ' ) dv„ ' = [ 5 (v„ ' -vj  <S  (Mo ' -Mo)  <5  (a  ' -a) /v„ ' 2] . 

v„'2dv„'dMo'da«  (2-21) 

where  Mo  is  the  cosine  of  0p. 

Substituting  Equations  (2—20)  and  (2—21)  into  Equation 
(2-19)  we  have 

D(V(vJ  ) = [BV7r^^^]exp[-B2(v^2+v„2-2v^v„Mo)  • (2-22) 

Then,  substituting  Equation  (2—22)  into  Equation  (2—17) 
and  integrating  yields 

D(V(v^)  )=[Bv^/(7T^/2^„)  ] {exp[-B2(v^-v„)2]- 

exp[-B2(v^+v„)2]}  (2-23) 

where 

Vh'=v/+v2-2V„vm  (2_24) 

and  M is  the  cosine  of  6. 

A careful  inspection  of  Equation  (2-23)  and  Equation  (2- 

13b)  shows  that  in  the  case  of  V„  equal  to  zero,  the 

distribution  is  the  one  expected  for  thermal  motion  only  and 
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in  the  limit  of  zero  temperature,  the  relative  velocity  is 
equal  to  v^. 

Substituting  Equation  (2-23)  into  Equation  (2-18),  yields 
a final  formulation  for  the  collision  frequency  that  includes 
the  effects  of  both  temperature  and  directed  velocity  of  the 
medium: 

F(v,V)=J  v^2a^(v^)N»B/(7r^/2v„)» 

(exp[-B2(v^-v„)2]-exp[-B2(v^+v„)2]  }dv^.  (2-25) 

As  can  be  observed.  Equation  (2-25)  is  a function  of  the 
neutron-nucleus  relative  velocity  (which  is  a function  of  both 
the  directed  velocity  and  temperature  of  the  medium) . This 
equation  can  now  be  used  in  the  derivation  of  effective  cross 
sections . 


2.3  Effective  Cross  Section 

The  formulation  of  the  effective  cross  section  is  derived 
from  the  concept  of  the  effective  collision  frequency. 

For  a stationary  target,  V=0,  and  v^=v.  Then,  the 
collision  frequency  is  as  given  in  Equation  (2-3) , 

F(v,0)=vNct^(v)  . (2-26) 

To  define  an  effective  collision  frequency  for  a medium 
with  moving  targets  we  can  equate  Equation  (2-14)  to  an 
equation  where  an  effective  collision  density  is  used  and  the 
target  is  assumed  to  be  at  rest. 

F(v,V)=  J v^»a*(v^)  •N*D(V)dV=Fg^^(v,0) 


(2-27) 
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Then,  this  concept  and  Equation  (2—26)  can  be  used  to  define 
Fgff  (v,0)=vNa’^g^^(v)  , and  (2-28) 

c^'eff(v)  = [F(v,V)  ]/vN.  (2-29) 

Thus,  the  effective  cross  section,  is  defined  to  yield 

the  correct  ensemble  collision  frequency  and  can  be  expressed 
as : 

(v)  = (l/v)  J v^«a^(v^)*D(V)dV.  (2-30) 

Using  the  formulation  for  F(v,V)  presented  in  Equation 
(2-25)  and  substituting  into  Equation  (2-29)  yields: 
t^'eff(v)=  (1/v)/  V^2^'(v^).B/(7t1/2v„)* 

(exp[-B2(v^-v„)2]-exp[-B2(v^+v„)2]  }dv^.  (2-31) 

This  Equation,  (2-31) , represents  the  general  definition 

of  the  total  effective  cross  section  for  the  case  studied, 
i.e.,  when  there  is  constant  speed,  monodirectional  motion  for 
the  directed  velocity,  takes  into  account  the  effects  of 

temperature  and  directed  velocity  of  the  target  atoms  and  all 
that  is  needed  for  its  evaluation  is  the  behavior  of  a^(v^)  and 
the  value  of  v^^.  This  equation  will  be  used  throughout  this 
work  in  the  derivation  of  other  expressions.  The  superscript 
t is  dropped  to  permit  the  generalization  of  the  development 
to  any  kind  of  cross  section. 

For  the  simple  case  where  aS(v^)=a®,  a constant.  Equation 
(2-31)  can  be  evaluated  to  get 

(v)  = (v„/v)  (aVBv„)  {exp[-(Bv„)2]/7ri/2 
+ [Bv„+l/  (2Bv„)  ]erf  (BVj,)  ) . 


(2-32) 
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If  the  target  atoms  are  not  submitted  to  a directed 
velocity  (limit  of  zero  directed  velocity)  and  only  thermal 
effects  are  present,  and  the  resultant  equation  is 

^^%ff  (V)  = (aVBv)  { exp  [-  (Bv) 

+ [Bv+1/ (2Bv)  ] erf  (Bv)  ) . (2-33) 

This  equation  is  equivalent  to  the  one  obtained  by  Oblow  [5] 
for  temperature  motion  only  and  has  also  been  used  in  many 
classical  references  to  describe  the  effective  scattering 
cross  section  in  the  thermal  region  for  a monoatomic  gas  with 
an  isotropic  Maxwellian  distribution  of  nuclear  velocities 
[32]  when  an  energy-independent  microscopic  scattering  cross 
section  is  used.  This  demonstrates  that  Equation  (2-31) 
predicts  the  correct  behavior  for  this  particular  case. 

A few  interesting  results  can  be  obtained  by  deriving  an 
equation  for  the  effective  cross  section  by  making  the  neutron 
velocity  equal  to  v,^,  and  considering  the  target  atoms  to  be 
subjected  only  to  thermal  motion.  For  this  case,  the  effective 
cross  section  is 

‘^eff(V„)=J  V^2^(v^).B/(7rl/\2j. 

{exp[-B^(v^-v„)2]-exp[-B2(v^+v„)2]  }dv^.  (2-34) 

Comparing  Equations  (2-31)  and  (2-34)  it  can  be  observed 
that  these  equations  are  related  as  follows: 

VCTgf^(v)=v„a^^f  (v„)  , and  consequently,  (2-35) 

(^eff  ('^)  = (v^/v)  (v„)  . (2-36) 

Equation  (2-36)  shows  that  the  actual  effective  cross 
section  for  a neutron  with  velocity  v,  interacting  with  target 
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atoms  subjected  to  both  directed  velocity  and  thermal  effects, 
is  equal  to  the  effective  cross  section  for  target  atoms 
subjected  only  to  thermal  motion,  scaled  by  the  ratio  v^/v, 
for  a neutron  with  velocity  v„.  This  shows  that  if  v„  is  known, 
the  problem  of  finding  the  effective  cross  section  for 
neutrons  interacting  with  target  atoms  with  both  thermal  and 
directed  velocity  is  equivalent  to  the  problem  of  finding  a 
scaled  effective  cross  section  for  neutrons  of  different  (but 
related)  energies  interacting  with  target  atoms  subjected  to 
only  thermal  motion. 

Equation  (2-34)  can  be  rewritten  by  using  the  definition 
of  an  effective  cross  section  for  neutrons  interacting  with  a 
Maxwellian  distribution  of  moving  targets  as 

^'eff(v„)  = (l/v„)J  M(V,)v^a‘^(v^)dV,,  (2-37) 

where  the  superscript  T stands  for  temperature  effects  only. 

Considering  that  the  integral  of  M(V^)  is  normalized  to 

f it  can  be  verified  that  if  <7^  (v^.)  =c/v^,  the  effective 
cross  section  is 

«7'eff(v„)  = (c/v„)  , (2-38) 

and  substituting  Equation  (2-38)  into  Equation  (2-36)  we  get 

«^eff(v)  = (c/v)  . (2-39) 

Equation  (2-39)  shows  that  1/v  cross  section  behavior  is 
not  affected  by  the  combination  of  directed  velocity  and 
thermal  velocity  for  the  target  atoms. 

Another  interesting  situation  to  be  examined  is  the  limit 
where  v„  tends  to  zero  (v  approximately  equal  to  V„)  and. 
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consequently,  from  Equation  (2-13b) , tends  to  V^.  In  this 
limit,  Equation  (2-37) , substitutinq  by  V^,  is  expressed  as 

^'eff  (Vm)  = ( l/v„)  J M(V,)  (V,)  dV, , (2-40) 

and  the  integral  that  appears  in  the  equation  can  be  evaluated 
as  a function  of  only.  The  resulting  effective  cross 
section  is 

^'eff(v„)  = (C/v„)  , (2-41) 

where  c'  is  the  result  of  the  integration.  This  expression  is 
similar  to  Equation  (2-38)  and,  by  substituting  this  result 
into  Equation  (2-36)  , it  is  obtained,  like  in  the  previous 
situation, 

f^eff(v)  = (c'/v)  . (2-42) 

It  can  be  observed  that  the  actual  effective  cross  section 
although  presenting  a discontinuity  for  v,,  equal  to 
zero  (from  Equation  (2-41)),  will  tend  to  be  proportional  to 
1/v  when  v„  is  very  small  (approaching  to  zero)  . 

Physically,  the  neutron  will  "see"  only  the  thermal 
motion  of  the  target  atom  and  the  actual  effective  cross 
section  will  depend  only  on  the  thermal  motion.  This  result  is 
equivalent  to  the  result  that  would  be  obtained  for  the  case 
of  temperature  motion  only  when  the  velocity  of  the  neutron  is 
much  smaller  than  the  velocity  of  the  target. 

To  demonstrate  that  the  derived  equation  (2-31)  for  the 
effective  cross  section  is  consistent  and  general,  it  must  be 
examined  in  the  limits  of  zero  directed  velocity  and 
temperature. 


zero 
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The  test  in  the  limit  of  zero  directed  velocity,  can  be 
done  by  making  V^=0,  and  in  consequence  v„=v,  in  Equation  (2- 
31)  to  get 

aeff(v)=  (1/V)J  •B/(7T^/2^)  • 

{exp[-B2(v^-v)2]-exp[-B2(v^+v)2]  }dv^.  (2-43) 

The  above  expression  for  is  equal  to  the  exact  equation 

derived  by  Oblow  considering  temperature  only. 

In  the  limits  of  zero  temperature,  Vj  tends  to  zero, 
B tends  to  infinity.  Applying  limits,  and  the 
following  properties  of  the  delta  function 

(1)  lini(B^)(B/y7r)  [exp(-BV)  ]=<S(x) 

(2)  5 (a-b) -5  (a+b)=2b<S  (a^-b^) 
to  Equation  (2-31)  yields 

"^eff(v)(T.o)=-^ [v,V(w„)]a(v^)  [«S(v^-v„)-5(v^+v„)  ]dv^  (2-44) 

=/ [v,V(wJ  ]a(v^)2v„6(v^2_^„2j^^^^ 

Letting  y=v^^  in  the  above  equation,  changing  variables 
and  solving  gives 

^eff  (T-0)^  ^ • (2-45) 

It  can  be  verified  that  in  this  case,  when  (T-^O)  , since  v =v 

r M ^ 

va^,,(v)=v^a(v^)  . (2-46) 

It  can  be  observed  that  if  Equation  (2-46)  is  multiplied 
by  N,  it  is  equivalent  to  Equation  (2-11)  , which  is  the 
expression  for  the  collision  frequency  with  a target  atom 
subjected  only  to  directed  velocity. 
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2.4  Scattering  Kernel 

In  order  to  complete  the  development  of  a generalized 
thermal ization  model  that  accounts  not  only  for  the  high 
temperature  but  also  for  the  directed  velocity  of  the  medium, 
it  is  necessary  to  derive  a new  scattering  kernel  that 
includes  these  two  features. 

The  first  thought  would  be  to  follow  the  path  of  Murray 
[31]  and  derive  an  exact  expression  for  the  new  scattering 
kernel  starting  with  the  derivation  of  the  velocity  transfer 
probability,  P(v'-^v)  . This  conventional  approach  results  in  an 
equation  for  the  scattering  kernel  which  has  two  additional 
degrees  of  freedom,  when  compared  to  the  equation  derived  for 
thermal  motion  only,  because  two  new  variables  appear  when  the 
target  atom  directed  velocity  is  included  in  the  derivation. 
One  variable,  the  cosine  of  the  angle  between  the  directed 
velocity  of  the  target  and  the  velocity  of  the  neutron,  is 
related  to  the  anisotropy  that  results  from  the  directed 
velocity;  the  other  variable  is  the  directed  speed  of  the 
target.  These  two  new  degrees  of  freedom  make  solution  of  the 
equation  rather  difficult. 

A solution  to  the  problem  of  finding  the  new  scattering 
kernel  is  suggested  by  a careful  observation  of  Equation  (2- 
36)  . From  Equation  (2-36)  it  was  deduced  that  the  actual 
effective  cross  section  when  both  thermal  and  directed 
velocities  effects  are  considered  is  equal  to  the  effective 
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cross  section  for  a neutron  with  velocity  when  only 
thermal  effects  are  included,  scaled  by  the  ratio  Vj^v.  This 
suggests  the  use  of  a new  reference  frame  that  moves  with 
velocity  instead  of  a static  frame.  If  it  is  recognized 
that  the  expressions  on  the  left  and  right  hand  sides  of 
Equation  (2-35) , when  multiplied  by  N represent  the  definition 
of  collision  frequency.  Equation  (2-35)  can  also  be  seen  as  a 
statement  that  the  collision  frequency  is  an  invariant  in  all 
frames . 

In  this  new  reference  frame  the  neutron  velocity  is 


(v-V^)-v^,  and  the  target  velocity  is  (2-47) 

(V-VJ=V,.  (2-48) 

The  effective  scattering  kernel  in  the  static  frame  can 

be  written,  by  definition,  as 

<?eff  (v'-v)=a^^^(v')P(v'->v)  . (2-49) 

Substituting  Equation  (2-36)  into  Equation  (2-49)  we  get 
^eff  ( V '-v)  = ( v„  '/V ' ) (V„  • ) P (V  • -V)  . (2-50) 

To  have  the  effective  scattering  kernel  (or  the  effective 


velocity  transfer  cross  section)  written  as  a function  of  v 

M ' 

it  now  remains  to  change  the  coordinates  of  the  velocity 
transfer  probability  from  the  static  frame  to  the  moving 
frame.  Considering  that  dv'=dV|^*,  in  spherical  velocity  space 
the  transformation  of  the  static  frame  velocity  transfer 
probability  to  the  moving  frame  yields 
P(v'-+v)  = (v'/v„'  )2p(v„'^v„)  , 

where  (v'/v^^')^  is  the  Jacobian  of  the  transformation. 


(2-51) 
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Substitution  of  Equation  (2-51)  into  Equation  (2-50) 
yields 

t^eff  (V ' ^v)  = ( V /v„  ' ) (v„  ' ) P (v„  ' -v„)  . (2-52) 

Equation  (2-52)  can  also  be  written  in  terms  of  the 
scatterinq  kernel  in  the  movinq  frame  as  follows 

^^eff  (V  ’-V)  = (V  '/v„  ' ) (v„  . (2-53) 

The  scattering  kernel  in  the  moving  frame,  a ,, (v„ '-^v„)  , 
depends  only  on  thermal  motion  and  an  exact  expression  for 
this  case  was  derived  by  Oblow  [5]  for  a static  frame.  The 
equation  derived  by  Oblow  can  be  modified  to  be  applied  to  the 
present  case  and  the  resultant  equation  is 

^eff(VM'^VM)  = (vVv„')2(B//7T)3((A+l)/2A)  (V„Vv„  ' ) 

(l/p)exp(-(BVA)  [ (A+1)v„v„'Ml"V„'2]  } 

^ (A.Dp/ZA^r^^'  (Vp)  exp  (-b2v^2)  i^{  2bV„v„  • 

[ (1-Ml^)  ( (v^p)2-(  (A+1)/2A)2)  (2-54) 

Substituting  Equation  (2-54)  into  Equation  (2-53)  we  get 
«^eff  ( V ' -V)  = ( V ' /V„  ' ) 2 ( B/Ztt)  2 ( ( A+1 ) /2A)  ( v//v„  ' ) 
(l/p)exp{-(BVA)  [ (A+1)v„v„'/Xl-v„'2]  } 

•^(A.1)p/2AVr‘^'  exp  (-B^V^^j  I^(  2b2v„V„  ' 

[(1-Ml^)  ((v^/p)2-((A+1)/2A)2)  (2-55) 

where, 

P=[v„'2+v„2-2v„v//iL]'/^  , (2-56) 

P=Vm-v„'  is  the  change  in  neutron  velocity, 

imaginary  Bessel  function  of  the  second  kind  of 
order  zero,  and 

A=M/m,  for  m equal  to  the  neutron  mass. 
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Equation  (2-55)  is  complete  as  it  stands  and  represents 
the  scattering  kernel  in  the  static  frame.  It  should  be 
noticed  that  the  anisotropy  added  by  the  directed  velocity  of 
the  target  atoms  is  included  in  and,  again,  there  are  two 
more  degrees  of  freedom  in  the  equation,  the  anisotropy  and 
the  speed  related  to  the  directed  velocity. 

If  the  neutron  distribution  in  the  moving  frame  is 
assumed  to  be  a Maxwellian,  the  condition  of  detailed  balance 
applies  in  this  frame.  This  can  be  shown  by  reversing  the 
scattering  direction  in  Equation  (2—54) , noticing  that  the 
integral  in  this  equation  is  invariant  under  such 
transformation.  Multiplying  each  scattering  kernel  by  the 
respective  Maxwellian  distribution  yields 
M ( v„ , T)  ( v„-v„  ■ ) =M  ( v„  ' , T)  ( v„  ' -v„)  . 


(2-57) 


CHAPTER  3 

SIMPLIFIED  NUMERICAL  DEVELOPMENT 
3.1  Introduction 

The  exact  derivation  of  the  equations  that  represent  the 
effective  cross  sections  and  the  scattering  kernel  in  a system 
where  the  target  atoms  are  subjected  to  both  thermal  velocity 
and  directed  velocity  are  presented  in  Chapter  2 . Except  under 
special  conditions,  these  equations  are  impossible  to  solve 
analytically  and  even  their  solution  by  numerical  methods  is 
very  difficult,  suggesting  possibly  a solution  employing  Monte 
Carlo  methods. 

The  development  of  a new  computer  code,  based  on  Monte 
Carlo  methods,  can  prove  to  be  very  time  consuming  and  the 
application  of  this  code,  even  for  a simple  case,  to  obtain  an 
estimate  of  the  effects  on  thermal  neutron  interaction  rates 
of  high  directed  velocity  and  temperature  of  the  target  atoms 
can  require  a great  amount  of  computer  time. 

Without  deviating  from  the  main  objective  of  this  work, 
i.e.,  the  deteirmination  of  the  effects  of  high  directed 
velocity  and  temperature  on  the  reaction  rates  of  the  target 
atoms  interacting  with  thermal  neutrons,  it  is  possible  that 
a more  straightforward  and  simpler  method  can  be  followed. 
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Unfortunately,  these  simpler  methods  usually  depend  on 
assumptions  that  make  the  final  solution  less  accurate  than 
the  rigorous  solution;  nevertheless,  for  many  cases,  if  the 
approximate  solution  is  reasonably  acceptable,  the  loss  in 
accuracy  and  generality  is  more  than  outweighed  by  the 
increased  simplicity  in  the  calculations  and  the  resultant 
reduction  in  computer  time.  What  is  sought  here  is  the 
development  of  a simplified  method  to  obtain  a reasonable 
approximate  solution  to  the  problem  using  an  existing  computer 
code. 


2-^2 The  Choice  of  a Computer  Code 

According  to  Section  3.1,  what  is  sought  is  an  existing 
computer  code  that  could  be  modified  to  include  the  effects  of 
both  thermal  velocity  and  directed  velocity  of  the  target 
atoms  in  the  computation  of  thermal  neutron  interaction  rates. 

The  chosen  code  is  the  BATTELLE  REVISED  THERMOS  (BRT-I) 
[33]  (see  Appendix  A).  This  code  incorporates  revisions  and 
modifications  to  the  original  THERMOS  [28]  code;  significant 
enhancements  include  a new  and  improved  library  tape, 
transport  kernel  modifications  for  the  treatment  of  void 
regions,  transport  kernel-cosine  current  calculations,  an 
increase  in  the  number  of  space  points  in  the  cylindrical  and 
slab  geometry  calculations,  and  an  option  for  cell  smear. 
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The  original  THERMOS  computer  code  solves  the  integral 
neutron  transport  equation  for  the  scalar  thermal  neutron 
density  and,  to  eliminate  the  angular  dependence  in  the 
thermal  neutron  density,  the  code  makes  the  assumption  that 
the  birth  rate  density  is  isotropic.  BRT-I  was  written  with 
emphasis  on  the  accurate  prediction  of  thermal  neutron 
reaction  rates  (which  makes  its  choice  appropriate  to  the 
problem  being  studied) . 

Some  special  features  in  BRT-I  (see  App.A)  include:  1) 
flexibility  in  the  way  that  data  are  input,  particularly  the 
ability  to  use  the  Brown  and  St.  John  approximation  to 
internally  calculate  the  scattering  cross  section;  when  this 
option  is  activated,  the  scattering  kernel  is  also  internally 
calculated  using  the  built-in  free  gas  scattering  kernel,  that 
includes  a feature  that  allows  the  user  to  input  the  desired 
temperature  for  each  isotope;  and  2)  the  ability  to  change  the 
built-in  default  energy  range.  These  special  features  combined 
with  the  other  basic  features  described  above  led  to  the 
choice  of  BRT-I  as  the  code  to  be  modified  to  include  the 
formulations  developed  in  this  work  to  account  for  the  effects 
of  the  directed  velocity  of  the  target  atoms  on  thermal 
neutron  interaction  rates. 
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Development  of  Approximations  to  the  Problem 

In  an  effort  to  develop  a simplified  method  to  solve  the 
according  to  what  is  stated  above,  some  assumptions 
and  approximations  must  be  made  to  the  system  that  involves 
neutrons  interacting  with  target  atoms  submitted  to  high 
directed  velocity  and  temperature  in  order  to  make  possible 
the  solution  of  the  developed  analytic  expressions  within  the 
BRT-I  code.  The  approach  suggested  in  Chapter  2 of  using  two 
velocity  frames  is  pursued  here. 

The  geometry  related  to  these  frames  is  shown  in  Figure 
3.1.  The  static  frame  is  the  frame  that  shows  the  neutron  and 
target  velocities  in  the  laboratory  system  and  the  moving 
frame  is  the  one  that  moves  with  the  medium's  directed 
velocity.  The  neutron  and  target  atom  velocities  in  both  the 
static  and  the  moving  frame  are  defined  in  Chapter  2.  The 
development  of  the  approach  has  its  starting  point  with 
Equation  (2-36) , where  v is  the  neutron  velocity  in  the  static 
frame  and  Vj^  is  the  neutron  velocity  in  the  moving  frame. 

In  the  development  of  the  method  it  must  be  recognized 
that  in  analyzing  the  effects  of  target  atom  directed  velocity 
on  neutron— nucleus  interactions,  two  extreme  conditions  can  be 
easily  identified:  1)  the  case  when  the  neutrons  have  not 

suffered  any  collision  yet  (uncollided  neutrons) ; and  2)  the 
case  when  the  neutrons  are  completely  entrained.  Since  these 
conditions  are  very  different  in  nature,  they  require 
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MOVING  FRAME 


Figure  3.1  - Neutron-nucleus  Geometry. 


43 


different  treatments.  The  method  of  solving  the  problem  for 
each  condition  is  as  follows: 

3.3.1  Uncollided  Neutrons 

Although  the  uncollided  neutron  flux  presents  a 
Maxwellian  distribution  in  the  static  frame,  from  Figure  3.1 
it  can  be  seen  that  the  neutron  distribution  in  the  moving 
frame  is  anisotropic.  The  variable  related  to  the  directed 
velocity  that  is  responsible  for  the  anisotropy  of  v in  the 
moving  frame,  /x,  can  be  clearly  identified  in  the  defining 
equation  of  v„.  Equation  (2-24) . The  target  atoms  present  a 
Maxwellian  distribution  in  the  moving  frame.  However,  the 
directed  velocity  of  the  medium  (target  atoms)  makes  both  the 
target  atom  velocity  distribution  and  the  collided  neutron 
flux  distribution  anisotropic  in  the  static  frame. 

It  can  be  seen  that  different  and  antithetical  conditions 
occur  in  the  different  frames.  The  selected  method  to  deal 
with  the  problem,  using  the  computer  code  BRT-I,  can  be 
divided  in  two  steps.  The  first  step  is  to  solve  the  integral 
neutron  transport  equation  to  determine  the  neutron  density 
and  interaction  rates  in  the  moving  frame  (where  the  target 
atoms  have  only  thermal  movement  and  the  neutron  velocity  is 
represented  by  v„)  . The  second  step  is  to  obtain  the  data 
(averaged  macroscopic  cross  sections)  for  the  case  when  the 
target  atom  directed  velocity  is  included  by  suitably 
converting  the  data  obtained  for  the  case  of  zero  target  atom 
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directed  velocity  (first  step).  In  the  second  step,  all  the 
averaged  cross  sections  obtained  from  the  output  of  BRT-I  are 
converted  back  to  the  domain  of  the  static  frame  (v) . In  this 
second  step  suitable  approximations  for  the  neutron  velocity 
related  to  the  moving  frame  need  to  be  obtained  since  v„  is  a 

M 

function  of  which  for  specific  values  of  v and  V„  includes 
one  more  degree  of  freedom  than  can  be  handled  by  BRT-I. 

In  order  to  accomplish  the  second  step  some 
approximations  need  to  be  made  with  respect  to  the  equation 
which  relates  the  neutron  speed  in  the  moving  frame,  v^^,  to 
the  neutron  speed  in  the  static  frame,  v,  the  directed 
velocity  of  the  target  atom,  V„,  and  n. 

From  Equation  (2-24) , 

v„=(V„2+v2-2V„vm)1/^  (3_1) 

In  the  above  equation,  only  the  positive  value  of  the  square 
root  is  retained  since  it  refers  to  a speed;  the  dependence  of 
v^  on  /i  is  clearly  seen. 

To  use  Equation  (3-1)  in  the  above  form  to  convert  the 
data  obtained  from  BRT-I  presents  some  difficulty;  this  is 
because  the  variation  of  /i  leads  to  many  different  values  for 
(speed)  for  each  speed  v,  even  when  keeping  constant. 
What  is  needed  is  a one-to-one  relationship  between  v and  v 

M 

in  order  to  make  the  change  from  the  moving  frame  to  the 
static  frame  (or  vice-versa)  straightforward.  As  mentioned 
before,  fj.  can  introduce  an  anisotropy  into  the  system.  Thus, 
one  approach  might  be  to  expand  Vj^  in  Legendre  polynomials  and 
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retain  at  least  the  first  two  terms;  the  first  term  (zeroeth 
moment)  giving  an  average  (isotropic)  value  of  v„  and  the 
second  term  (first  moment)  giving  a measure  of  the  linear 
anisotropy  involved.  However,  this  expansion  cannot  be  used 
since,  again,  there  is  a ^-dependence  in  the  second  term. 

Recognizing  that  BRT-I  calculates  the  scalar  neutron 
density  and,  hence,  works  with  speeds  instead  of  velocities 
suggests  the  development  of  an  approach  replacing,  for  each 
value  of  v,  the  corresponding  velocity  distribution  of  v„  over 
/X  by  a single  velocity  vector  averaged  over  n and  then  working 
with  an  averaged  speed;  this  averaged  vector,  v^,  is  assumed 
to  be  representative  of  the  neutron  velocity  distribution  in 
the  moving  frame.  Through  this  process,  the  averaged  speed  v^ 
in  the  moving  frame  can  be  directly  related  to  the  neutron 
speed  V in  the  static  frame. 

The  averaged  speed  corresponding  to  the  averaged  velocity 
Vg  is  obtained  by  averaging  Equation  (3-1)  over  /x.  Considering 
that  there  is  azimuthal  symmetry  in  the  neutron  velocity 
distribution  in  the  static  frame  and  that  /x  can  vary  from  —1 
to  1,  the  averaged  speed  is  given  by 

V3=(1/2)  (v2+V„2-2W„/x)'/^dM  (3-2) 

which  yields 

^a=v+(V„V3v)  (3_2a) 

for  V greater  than  or  equal  to  V„,  and 
^a=(vV3V„)+V„ 

for  V smaller  than  or  equal  to  V^. 


(3-2b) 
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Figure  3.2  shows  the  modified  geometry  that  represents 
the  velocity  distribution  of  over  ji  replaced  by  a single 
velocity  vector  averaged  over  It  should  be  noticed  that, 
because  of  the  azimuthal  symmetry  of  v,  for  specific  values  of 
V,  Vg  and  V„,  many  velocity  vectors  exist  but  the  important 
point  here  is  that  the  speed  v^  is  always  the  same,  allowing 
a one-to-one  relationship  between  v^  and  v. 

It  is  interesting  to  notice  that,  according  to  Equation 
(3-2a),  for  V much  greater  than  V„,  v^  tends  to  v.  At  this 
point,  an  analysis  of  shows  that  if  is  comparable  to  v, 
the  problem  turns  out  to  be  the  one  expected  when  the  target 
atoms  are  submitted  only  to  thermal  movement  and  in  this  case 
the  problem  can  be  solved  directly  in  the  static  frame. 

From  Equation  (3-2b)  it  can  be  observed  that  if  is 
much  greater  than  v,  v,  tends  to  . 

Another  conclusion  can  be  drawn  from  Equations  (3-2a)  and 
(3-2b):  that  v^  is  always  greater  than  v and  V„.  This  shows 

n 

that  the  neutron  spectrum  in  the  moving  frame  is  always  harder 
than  in  the  static  frame. 

After  this  approximation  is  made  for  v„,  the  problem  in 
the  moving  frame  is  now  simplified  to  the  case  where  for  each 
neutron  speed  (velocity),  v,  in  the  static  frame,  neutrons 
with  a corresponding  speed,  v^,  are  interacting  with  a 
Maxwellian  distribution  of  target  atoms  characterized  by  a 
temperature  T.  The  replacement  of  which  is  /l4— dependent,  by 


\ 

\ 
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Figure  3.2  - Modified  Neutron-nucleus  Geometry. 
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Vg,  which  is  ii- independent , provides  the  desired  one-to-one 
relationship  between  v and  mentioned  above. 

The  second  step  in  the  above  method  can  be  accomplished 
if  it  is  observed  from  Figure  3.2  that  the  relative  velocity 
in  the  moving  frame  (between  the  average  neutron  velocity  in 
the  moving  frame,  v^,  and  the  thermal  velocity  component  of 
the  target  atom  velocity,  V^)  and  the  relative  velocity  in  the 
static  frame,  (between  the  neutron  velocity,  v,  and  the 
nucleus  velocity,  V)  is  represented  by  the  same  velocity 
vector,  v^g.  This  provides  the  required  relationship  between 
the  effective  cross  section  in  the  moving  frame  and  in  the 
static  frame,  as  shown  in  Chapter  2,  Section  2.3.  In  this  case 
Vj.g,  the  relative  speed,  is  given  by 

Vpa=(V32+V,2-2v^V,Mi)'/^  (3-3) 

where  is  the  cosine  of  the  angle  between  and  V^. 

Using  this  approximation  in  Equation  (2-34), v„  is 
replaced  by  v^  and  v_,  by  v^^  yielding 
‘^eff  (^a)  =/  (V^J  . B/  (7t1/2v^2)  . 

(exp[-B2(v^g-vj2]_exp[-B2(v^g+vj2]  jdv^^,  (3-4) 
and  Equation  (2-36)  can  be  rewritten  as 

^eff(v)  = (Vv)CJeff(V3)  . (3-5) 

Using  Equation  (3-5),  the  cross  sections  obtained  from 
the  solution  of  the  problem  in  the  moving  frame  (first  step) 
can  now  be  transferred  to  the  static  frame  (v  domain) . 

Until  this  point,  approximations  were  made  to  facilitate 
the  transfer  of  data  from  the  moving  frame  to  the  static  frame 
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and  vice-versa.  Now,  the  validity  of  using  BRT-I  to  solve  the 
problem  in  the  moving  frame  must  be  established.  One  of  the 
primary  conditions  upon  which  the  BRT-I  computer  code  solution 
to  the  integral  neutron  transport  equation  is  based,  is  that 
there  be  a nearly  isotropic  thermal  neutron  birth  rate 
density.  The  birth  rate  density  is  nearly  isotropic  if  either 
the  scattering  process  is  nearly  isotropic  or  if  the  neutron 
number  density  distribution  is  nearly  isotropic.  For  the  case 
of  uncollided  neutrons,  even  if  the  birth  rate  density  is 
nearly  isotropic  in  the  static  frame,  it  should  be  apparent 
that  the  condition  of  a nearly  isotropic  birth  rate  density  is 
not  fulfilled  in  the  moving  frame  for  the  case  when  a large 
directed  velocity  of  the  target  atoms  is  present  since  this 
directed  velocity  introduces  an  anisotropy  in  the  distribution 
of  thermal  neutrons  in  the  moving  system.  Nevertheless,  it 
should  also  be  recognized  that  what  is  needed,  according  to 
the  selected  approach  (see  Equation  (3-5)),  is  the  effective 
cross  section  in  the  moving  frame  for  neutrons  with  speed  v^, 
interacting  with  target  atoms  having  a Maxwellian  distribution 
characterized  by  the  temperature  T,  and  having  a relative 
speed  v^g. 

In  the  /i-averaging  process,  the  distribution  is 
replaced  by  the  monodirectional  averaged  neutron  velocity  v^. 
It  should  be  noticed  that  for  each  given  magnitude  of  v there 
^ (distribution  which  is  replaced  by  the  corresponding  v^, 
all  of  them  in  the  same  direction. 
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The  replacement  of  the  distribution  by  leads  to  a 
moving  frame  in  which  neutrons  with  monodirectional  velocities 
are  interacting  with  target  atoms  which  present  a Maxwellian 
velocity  distribution.  It  is  known  that  the  results  for  the 
neutron  reaction  rate  when  the  velocity  distribution  of  the 
neutrons  is  isotropic  and  there  is  monodirectional  target  atom 
boundary  motion  are  the  same  as  when  the  boundary  motion 
velocity  distribution  of  the  target  atoms  is  isotropic  and 
there  is  monodirectional  motion  for  the  neutrons  [18],  This 
leads  to  the  creation  of  a eguivalent  moving  frame  where  the 
have  isotropic  distribution  and  the  target  atoms  have 
monodirectional  boundary  motion.  Then,  the  reaction  rate 
obtained  for  the  above  condition  is  the  same  as  the  one 
*^^i-^iried  in  the  moving  frame  in  which  neutrons  present 
monodirectional  velocities.  It  should  be  recognized  that  for 
both  the  moving  frame  and  the  equivalent  moving  frame,  the 
boundary  velocity  is  zero  and,  the  equivalent  moving  frame 
being  analyzed  by  BRT-I  has  a Maxwellian  distribution  of 
target  atoms  and  an  isotropic  distribution  of  neutrons.  The 
condition  of  isotropic  neutron  distribution  automatically 
leads  to  a condition  in  which  the  BRT-I  isotropic  birth  rate 
requirement  is  satisfied  and,  thus,  the  use  of  BRT-I  for  the 
equivalent  moving  frame  is  satisfied. 

Considering  the  above  approach,  it  needs  to  be  recognized 
that  after  the  neutrons  have  collided  once,  they  will  have  an 
anisotropic  distribution  in  the  static  frame  and  in  the  moving 
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frame,  and  this  approach  cannot  be  used  because  the  averaged 
speed  Vg  is  calculated  under  the  assumption  of  an  isotropic 
distribution  of  neutron  velocities  in  the  static  frame. 


3.3.2  Completely  Entrained  Neutrons 

As  already  mentioned  (see  Section  3.3,  Item  3.3.1),  the 
uncollided  neutrons  have  a Maxwellian  distribution  in  the 
static  frame.  This  distribution  corresponds  to  a velocity 
in  the  moving  frame  which  resembles  a Maxwellian- 
shape  distribution  shifted  to  higher  velocities  (this  can  be 
verified  by  making  the  change  of  variables  from  the  static 
frame  to  the  moving  frame) . 

As  the  neutrons  suffer  scattering  collisions  they  become 
more  and  more  entrained  in  the  medium,  as  represented 
schematically  in  Figure  3.3.  In  this  figure  the  neutrons  are 
shown  to  start  with  an  isotropic  distribution  in  the  static 
frame,  to  assume  intermediate  distributions,  and  after  being 
completely  entrained,  to  have  an  isotropic  distribution  in  the 
moving  frame. 

The  isotropic  distribution  assumed  by  the  neutrons  in  the 
moving  frame,  after  being  completely  entrained,  is  similar  to 
the  one  that  they  would  have  in  the  static  frame  in  the 
absence  of  the  directed  velocity  of  the  medium.  The  knowledge 
that  the  neutrons  assume  an  isotropic  distribution  in 
themoving  frame  permits  the  analysis  of  the  case  when  the 
neutrons  are  completely  entrained. 


MOVING  FRAME 


Figure  3.3 


Schematic  Geometry  for  Neutrons 
Undergoing  Entrainment. 
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As  observed  before,  the  velocity  distributions  for  the 
intermediate  states  are  not  readily  obtained  because  the 
determination  of  an  averaged  speed  in  the  moving  frame,  v , is 
very  difficult.  However,  the  condition  where  the  neutrons  are 
completely  entrained  can  be  treated  using  an  approach 
different  from  the  one  used  in  Item  3.3.1  above  for  the  case 
of  uncollided  neutrons.  This  can  be  accomplished  if  it  is 
noted  that  the  relative  velocity  v„  in  the  laboratory  system 
(which  is  also  the  neutron  velocity  in  the  moving  frame) , 
between  the  neutron  velocity  v in  the  laboratory  system  and 
the  medium's  directed  velocity,  V„,  is  isotropic  in  the  moving 
frame. 

To  be  completely  entrained  the  neutrons  must  suffer  a 
number  of  scattering  collision  where  they  are  subjected  to 
both  the  effects  of  the  directed  velocity  of  the  target  atoms 
and  the  thermal  movement  of  the  target  atoms.  The  inclusion  of 
the  directed  velocity  of  the  target  atoms  makes  the  scattering 
process  linearly  anisotropic,  with  the  directed  velocity  of 
the  target  atoms  dictating  the  preferred  direction.  When  the 
neutrons  become  completely  entrained  they  are  in  dynamical 
equilibrium  with  the  target  atoms  in  the  moving  frame,  as  far 
as  the  directed  velocity  of  the  target  atoms  is  concerned. 
Then,  in  the  moving  frame,  as  already  mentioned,  they  behave 

they  were  subjected  only  to  the  thermal  movement  of  the 
target  atoms,  without  considering  the  effects  of  the  directed 
velocity  of  the  target  atoms.  In  this  situation  of  thermal 
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neutrons  subjected  only  to  thermal  motion,  the  scattering 
process  can  be  reasonably  assumed  to  be  nearly  isotropic. 

The  first  step  in  the  method  is  again  performed  using 
BRT-I  in  the  moving  frame.  The  requirement  of  isotropic  (or  at 
least  linearly  anisotropic)  birth  rate  density,  for  the 
situation  of  completely  entrained  neutrons  is  accomplished  if 
it  is  recognized  that  the  scattering  process  can  be  assumed  to 
be  no  worse  than  linearly  anisotropic,  as  explained  before. 
This  situation  should  present  an  exact  result  since  it  is  no 
longer  necessary  to  replace  the  now  isotropic  v„  with  an 
averaged  v^.  In  the  second  step,  when  the  results  from  the 
moving  frame  are  transported  back  to  the  static  frame,  it  is 
again  necessary  to  obtain  a one-to-one  relationship  between  v 
since  V is  not  isotropic.  Now,  however,  instead  of 
averaging  v^,  v is  averaged  over  the  angle  that  defines  its 
anisotropy,  In  so  doing,  it  should  be  noticed  that  for 

each  v„  there  are  many  corresponding  neutron  velocities,  v, 
each  one  characterized  by  a speed  v and  an  angle,  according  to 
the  values  assumed  by  leading  to  a neutron  velocity 

distribution  for  each  v„.  The  above  described  procedure 
replaces  the  speed  distribution  (range)  corresponding  to  the 
neutron  velocity  distribution  in  the  static  frame  by  a single 
averaged  speed,  v*.  Although  this  approximation  can  introduce 
some  errors,  it  is  common  practice  to  replace  speed  ranges  by 
averaged  speeds  with  reasonable  results. 
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The  averaging  of  v over  ii^  is  performed  using  a procedure 
similar  to  the  one  employed  in  the  averaging  of  in  Item 
3.3.1  above.  Figure  3.3  shows  that  for  each  set  of  velocity 
vectors  v^^  and  there  can  be  many  different  values  for  the 
neutron  speed  v which  is  defined  as 

''=(V+V-2v„V,Mj)’«  (3_6) 
where  is  the  cosine  of  the  angle  between  v„  and  V . 

The  averaged  value  of  v to  be  considered  in  the 
approximation  in  order  to  have  a one-to-one  relationship 
between  the  neutron  speed  in  the  moving  frame  and  the  neutron 
speed  in  the  static  frame,  v*,  is  obtained  by  averaging 

Equation  (3-6)  over  /Xj,  with  varying  from  -1  to  1.  This 
yields 


(3-6a) 


(3-6b) 


v„+(V„2/3v„) 

greater  than  or  equal  to  V„,  and 
v*=  Vm+(v„V3V„) 
for  V|^  smaller  than  or  equal  to  V„. 

The  above  equations  show  that  v*  is  always  greater  than 
both  v„  and  V„  (a  result  similar  to  the  one  obtained  in  Item 
3.3.1  above  for  v^,  v and  V„)  . 

this  case  the  relative  speed,  v^,,  is  given  by 

or  (3^, 

Vr.=  (v"+V^-2VVM,)’'^  ,3^, 

where  ^3  is  the  cosine  of  the  angle  between  v„  and  and  is 

the  cosine  of  the  angle  between  v*  and  V,  as  shown  in  Figure 
3.4. 
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Figure  3.4  - Modified  Geometry  for  Completely 
Entrained  Neutrons. 
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As  in  the  preceding  section,  the  relative  velocity  is  the 
same  for  both  frames  and  Equation  (3-5)  , modified  for  this 
case,  yields 

= • (3-7) 

Summarizing,  to  handle  the  problem  when  the  neutrons  are 
completely  entrained,  the  integral  neutron  transport  equation 
will  be  solved  in  the  moving  frame  using  v^  and  the  results 
will  then  be  transferred  back  to  the  static  frame  where  the 
neutron  speed  is  represented  by  the  averaged  speed  v*. 

3.4  Analysis  for  Selected  Limiting  Conditions 

An  analysis  for  selected  limiting  conditions  is  performed 
separately  for  the  two  extreme  cases,  uncollided  neutrons  and 
completely  entrained  neutrons,  as  follows: 

3.4.1  Uncollided  Neutrons 

Before  performing  any  analysis  for  selected  particular 
conditions  it  is  necessary  to  recall  the  equations  for  the 
neutron  speed  in  the  static  frame,  as  a function  of  v and  V„. 
This  is  because  the  problem  will  be  first  solved  in  the  moving 
frame  in  terms  of  v^  and  and  then  converted  back  to  the 
static  frame  in  terms  of  v. 

The  analysis  must  compare  the  results  in  the  static  frame 
obtained  with  the  developed  approach  with  the  results  obtained 
for  similar  cases  in  Chapter  2 where  the  exact  approach  was 
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formulated.  Two  cases  can  be  distinguished  from  Equations  (3- 
2a)  and  (3 -2b)  . One  when  v is  greater  than  V„  and  another  when 

V is  smaller  than  V„.  These  conditions  are  shown  in  Figures 
3.5  and  3.6,  respectively. 

The  equation  for  v equal  to  or  smaller  than  V,,  is  given 

by: 

v=(3V„(v^-V„))V2 

and,  in  this  case,  Vg>V„>v;  the  positive  sign  of  the  square 

root  is  retained  since  there  is  no  meaning  to  a negative 
speed. 

When  V is  equal  to  or  greater  than  V„,  the  equation  for 

V is  given  by: 

v=(l/2) 

and,  in  this  case,  v^>v>V„.  Since  Equation  (3-9)  must  hold  for 
V„=0,  in  which  case  v^=v,  the  negative  value  for  the  square 
root  is  discarded. 

It  is  now  possible  to  examine  the  behavior  predicted  by 

the  developed  approach  for  some  selected  particular 
conditions . 

3.4. 1.1  The  limiting  condition  in  which  V„  tends  to  7.e>rn 


From  Equation  (3-9) , it  can  be  seen  that  v^  tends  to  v 
and  the  problem  reduces  to  the  one  of  thermal  movement  only. 
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Figure  3.5  - Distribution  of  in  the  Moving  Frame  for 
V Greater  than  Vj^. 


MOVING  FRAME 


Figure  3.6  - Distribution  of  in  the  Moving  Frame 
for  V Smaller  than 
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3. 4. 1.2  For  much  greater  than  v and  V.. 

0 ' M 

In  this  case,  the  neutron-nucleus  relative  velocity  can 
be  approximated  by  in  both  frames.  The  effective  cross 
section  in  the  moving  frame  is  proportional  to  l/v^  and,  from 
Equation  (3-5) , the  effective  cross  section  in  the  static 
frame  is  proportional  to  1/v;  this  is  the  expected  behavior 
without  making  any  assumptions,  as  demonstrated  in  Chapter  2. 

3. 4. 1.3  The  limiting  condition  when  is  much  greater  than 
If  Vm  is  much  greater  than  V^,  from  either  Equation  (3 -2b) 

or  Equation  (3-8)  v^  is  always  greater  than  V„,  and  it  is 
apparent  that  v^  is  also  much  greater  than  V^.  Consequently, 
in  this  case,  the  neutron-nucleus  relative  speed  in  the  moving 
frame  is  expected  to  approach  being  temperature  independent. 
In  fact,  in  the  limit  of  very  large  V„,  the  neutron-nucleus 
relative  speed  approaches  V„  and  from  Equation  (3-2b) , v^  also 
approaches  V^.  The  effective  cross  section  in  the  moving 
frame,  thus  proportional  to  l/v^  (or  1/V„)  and, 

from  Equation  (3-5),  a^ff(v)  is  proportional  to  1/v.  Equation 
(3-5)  IS  thus  equivalent  to  Equation  (2-46)  and  the  behavior 
predicted  using  the  approach  developed  for  uncollided  neutrons 
is  the  same  behavior  predicted  in  Chapter  2,  without  making 
any  assumptions,  under  the  same  conditions. 


3.4.2  Completely  Entrained  Neutrons 
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3. 4. 2.1  The  limiting  condition  in  which  V,,  t,P,nds  to  tipt-o 
From  Equation  (3-6a)  , it  can  be  seen  that  v*  tends  to  v,,, 

which  by  definition  tends  to  v when  V„  tends  to  zero.  Thus, 
there  is  no  need  to  change  frames  and  the  problem  reduces  to 
the  one  of  thermal  movement  only. 

3. 4. 2. 2 For  much  greater  than  v* 

Equations  (3-6a)  and  (3-6b)  demonstrate  that  v*  is  always 
greater  than  both  v„  and  V„.  Hence,  in  this  case  is  also 
much  greater  than  both  v„  and  V„.  The  effective  cross  section 
in  the  moving  frame  is  proportional  to  l/v„  and,  from  Equation 
(3-7) , the  effective  cross  section  in  the  static  frame  is 
proportional  to  1/v*,  as  expected. 

3.,  4. 2. 3 The  limiting  condition  when  is  much  areatP.r  th^n 

In  this  case  the  relative  velocity  in  the  moving  frame  is 
expected  to  be  temperature-independent.  In  fact,  v^*  tends  to 
v„  which  in  turn  tends  to  V„.  The  effective  cross  section  in 
the  moving  frame,  cr^^^Cv^),  is  proportional  to  l/v„  (or  1/V„) 
and,  from  Equation  (3-7),  a^,,(v*)  is  proportional  to  v*  as 
expected.  Equation  (3-7)  is  thus  equivalent  to  Equation  (2-46) 
and  the  behavior  predicted  using  the  approach  developed  for 
completely  entrained  neutrons  is  the  same  behavior  predicted 
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in  Chapter  2,  without  making  any  assumptions,  under  the  same 
conditions. 

The  above  shows  that  for  these  selected  limiting 
conditions  the  behavior  predicted  under  the  assumptions  made 
for  the  cases  of  uncollided  neutrons  and  completely  entrained 
neutrons  is  consistent  with  the  behavior  predicted  by  the 
exact  approach. 

Although  Items  3.4.1  and  3.4.2  above  represent  different 
physical  conditions,  the  similarity  in  the  conclusions  for 
analogous  situations  should  be  noted. 

^ — Preliminary  Test  of  the  Developed  Appr-n;.r^h 

The  effects  of  both  high  directed  target  atom  velocity 
and  temperature  on  the  interaction  rates  with  neutrons  hasn't 
been  treated  in  the  past  and,  hence,  only  a few  related 
theoretical  works  were  found.  The  only  data  available  for 
comparison  are  from  the  work  of  Perkins  [18].  in  this  work  the 
computer  code  SOPHIST  III  was  utilized  to  calculate  the 
effective  neutron  scattering  cross  section  for  hydrogen  in  a 
medium  subjected  to  a speed  (VJ  of  lo^  m/s.  Perkins  first 
examined  hydrogen  at  0 °K,  assuming  a constant  relative  cross 
section  of  20.3  barns  and  then  he  repeated  the  calculations, 
for  the  same  V^,  for  hydrogen  at  2500  °K. 

Perkins  analyzed  two  conditions:  1)  the  case  in  which  the 
medium  (target  atoms)  is  moving  isotropically  with  constant 
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speed  (expanding  or  contracting)  about  the  point  in  space 
where  the  collision  occurs;  and  2)  the  case  in  which  the 
medium  is  moving  with  monodirectional  velocity  (a  single 
directed  velocity  for  the  target  atoms) . In  both  cases  thermal 
effects  are  considered  and  the  velocity  of  the  neutrons  is 
assumed  constant.  For  the  case  of  monodirectional  boundary 
motion,  he  demonstrated  that  the  results  for  the  neutron 
reaction  rate  when  the  velocity  distribution  of  the  neutrons 
is  isotropic  and  there  is  monodirectional  target  atom  motion 
are  the  same  as  when  the  velocity  distribution  of  the  target 
atoms  is  isotropic  and  there  is  monodirectional  motion  for  the 
neutrons. 

To  demonstrate  his  approach  Perkins  presented  in  two 
examples  the  results  from  calculations  for  the  effective 
scattering  cross  section  of  hydrogen,  at  0 °K  and  2500  °K.  He 
compared  the  resulting  cross  sections  for  the  cases  of 
temperature  motion  only,  isotropic  boundary  motion,  and 
monodirectional  boundary  motion;  in  this  last  case  he 
presented  the  results  for  the  conditions  where  the  cosine  of 
the  angle  between  the  neutron  velocity  and  the  target  atom 
velocity  is  1 and  -l. 

His  results  are  presented  in  Figures  3 and  4 of  reference 
18.  From  the  results  it  can  be  inferred  that  the  condition  of 
uncollided  neutrons  was  applied,  since  the  medium's  directed 
velocity  is  lO'^  m/s  and  Perkins  computed  neutron  energies  far 
below  the  energy  corresponding  to  this  velocity. 
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For  the  uncollided  neutron  condition,  the  approach 
developed  in  this  work  is  for  the  case  of  a single  directed 
velocity  for  the  target  atoms  (also  subjected  to  thermal 

effects)  and  an  isotropic  distribution  for  the  neutrons.  The 
results  are  to  be  compared  with  Perkins'  results  for  the  case 
of  isotropic  boundary  (target  atom)  motion  and  a single 

directed  velocity  of  the  neutrons  since  this  is  equivalent  to 
the  situation  where  there  is  isotropic  distribution  of 
neutrons  and  a single  directed  velocity  of  the  target  atoms  as 
mentioned  above. 

To  perform  the  comparison,  some  neutron  energies  were 
first  selected,  then  for  each  selected  energy  the  respective 
neutron  velocity  and  the  target  atom  (or  medium)  directed 
velocity  were  used  to  calculate  v^  (see  Equations  (3-2a)  and 
(3-2b) ) . Using  v^,  the  effective  cross  section  in  the  moving 
frame,  considering  only  the  thermal  motion,  was 

obtained.  Then,  using  Equation  (3-5)  the  effective  cross 

section  in  the  laboratory  system  for  a neutron  with  speed  v, 

^ eff  ' considering  both  the  target  atom  directed  velocity 
and  thermal  velocity  distribution  was  calculated.  The  results 
are  presented  in  Table  3.1  and  Table  3.2,  where 
E^  (eV) : neutron  energy  in  the  static  frame 
V (m/s) : neutron  speed  in  the  static  frame 

(^/®) • neutron  speed  (averaged)  in  the  moving  frame 
(barns):  effective  neutron  scattering  cross 

section  in  the  moving  frame  (thermal  movement) 
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Table  3.1  - Comparison  between  the  Effective  Scattering  Cross 
Sections  for  Hydrogen  at  2500  °K,  with  V =10^  m/s 
(uncollided  neutrons) . 


En 

(eV) 

V 

(m/s) 

/a 

(m/s) 

^'eff 

. (vj 

^%ff 

(V). 

(V)n 

V^V 

0.01 

1386 

10064 

26.0 

190.0 

188.8 

7.261 

0.03 

2400 

10192 

23.3 

100.0 

98.9 

4.247 

0.05 

3098 

10320 

23.2 

80.0 

77.3 

3.331 

0.10 

4382 

10640 

23.1 

58.0 

56.1 

2.428 

0.30 

7589 

11920 

22 . 0 

36.0 

34.6 

1.571 

0.50 

9798 

13200 

21.1 

29.0 

28.4 

1.347 

1.00 

13856 

16261 

21.0 

25.0 

24.7 

1.174 

2.00 

19596 

21297 

20.3 

21.5 

22.0 

1.087 

5.00 

30984 

32060 

20.3 

20.7 

21.0 

1.034 

10.00 

43818 

44578 

20.3 

20.4 

20.6 

1.017 

Note:  Table  3.1  compare  results  obtained  from  Perkins 
[18]  and  obtained  using  the  developed  approach. 
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Table  3.2  - Comparison  between  the  Effective  Scattering  Cross 
Sections  for  Hydrogen  at  0 °K,  with  V =10^^  m/s 
(uncollided  neutrons) . 


En 

(eV) 

V 

(m/s) 

(m/s) 

(vj 

(v)o 

(V), 

v^/v 

0.01 

1386 

10064 

20.3 

150.0 

147.4 

7.261 

0.03 

2400 

10192 

20.3 

86.0 

86.2 

4.247 

0.05 

3098 

10320 

20.3 

67.0 

67.6 

3.331 

0.10 

4382 

10640 

20.3 

49.0 

49.3 

2.428 

0.30 

7589 

11920 

20.3 

32.0 

31.9 

1.571 

0.50 

9798 

13194 

20.3 

28.0 

27.3 

1.347 

1.00 

13856 

16261 

20.3 

23.0 

23.8 

1.174 

2 . 00 

19596 

21297 

20.3 

21.0 

22.1 

1.087 

5.00 

30984 

32060 

20.3 

20.5 

21.0 

1.034 

10.00 

43818 

44578 

20.3 

20.4 

20.6 

1.017 

Note:  Table  3.2  compare  results  obtained  from  Perkins 
[18]  and  obtained  using  the  developed  approach. 
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(barns):  effective  neutron  scattering  cross 
section  in  the  static  frame,  obtained  from  Perkins 

(barns);  effective  neutron  scattering  cross 
section  in  the  static  frame,  obtained  using  the  approximation 
developed  in  this  work. 

Figures  3 . 7 and  3 . 8 are  graphical  representations  of  the 
results.  Since  the  results  from  this  work  and  Perkins'  work 
are  in  such  good  agreement,  they  appear  as  the  single  dashed 
line  in  Figures  3.7  and  3.8. 

It  should  be  mentioned  that  the  presented  values  in 
Tables  3.1  and  3.2  from  Perkins'  work  may  not  exactly 
represent  his  results  since  they  were  picked  off  or  extracted 
from  small  log-log  plots. 

From  the  results  presented  in  the  tables  it  can  be 
concluded  that  the  static  frame  effective  cross  sections 
obtained  by  using  the  effective  cross  section  for  uncollided 
neutrons  in  the  moving  frame  (thermal  movement  only)  and  then 
translating  these  values  back  to  the  static  frame  are  in  good 
agreement  with  the  results  presented  by  Perkins. 

Although  there  are  no  known  data  available  to  compare  for 
the  case  of  completely  entrained  neutrons,  in  order  to  observe 
the  predicted  behavior  using  the  developed  approximate  method, 
a few  calculations  were  performed  for  the  effective  neutron 
scattering  cross  section  for  hydrogen,  at  0 °K  and  2500  °K, 
respectively,  with  a target  atom  directed  velocity  of  10^  m/s. 
This  velocity  is  chosen  to  allow  a later 
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Figure  3.7  - Effective  Scattering  Cross  Section  for  Hydrogen  at 
(uncollided  neutrons) . 
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Figure  3.8  - Effective  Scattering  Cross  Section  for  Hydrogen  at  2500 
(uncollided  neutrons) . 
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the  results  obtained  for  the  condition  of  uncollided  neutrons 
and  completely  entrained  neutrons.  The  results  of  those 
calculations  are  presented  in  the  Table  3.3  and  Table  3.4,  and 
the  variables  that  appear  in  these  tables  are  defined  as 
follows : 

(^/®) • neutron  speed  in  the  moving  frame. 

V*  (m/s) : neutron  speed  (averaged)  in  the  static  frame 
(barns);  effective  neutron  scattering  cross 
section  in  the  moving  frame  (thermal  movement) . 

^®eff(v*)  (barns):  effective  neutron  scattering  cross 
section  in  the  static  frame  obtained  using  the  developed 
approach. 

In  this  case,  for  selected  values  of  v„  (i.e.,  for  some 
V and  V^) , the  data  for  the  effective  cross  section,  a®  Cv  ) 
in  the  moving  frame  (for  thermal  movement  only)  is  obtained 
from  the  log-log  plots  presented  by  Perkins  for  hydrogen  at 
2500  K and  0 °K  considering  only  temperature  motion.  These 
data  were  obtained  from  Perkins'  work  in  order  to  have 
consistency  in  the  comparison  between  the  condition  of 
uncollided  neutrons  and  completely  entrained  neutrons  that 
will  be  carried  out  in  the  development  of  this  work. 

According  to  the  established  approach,  the  neutron  speed 
in  the  static  frame,  v*,  is  calculated  from  Equations  (3-6a) 
and  (3-6b)  and  the  effective  cross  section  in  the  laboratory 
frame,  CTeff(v*)p,  is  calculated  using  Equation  (3-7). 
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Table  3.3  - Effective  Scattering  Cross  Section  for 
Hydrogen  at  0 °K,  with  V„=10^  m/s 
(completely  entrained  neutrons) . 


En 

(eV) 

(m/s) 

V* 

(m/s) 

^%ff 

(v„) 

ef  f 

(V*)n 

v^v* 

0.53 

1386 

10064 

20.3 

2.79 

0.1377 

0.54 

2400 

10192 

20.3 

4 . 78 

0.2355 

0.56 

3098 

10320 

20.3 

6.09 

0.3002 

0.59 

4382 

10640 

20.3 

8.36 

0.4119 

0.74 

7589 

11920 

20.3 

12.92 

0.6365 

0.91 

9798 

13194 

20.3 

15.07 

0.7424 

1.38 

13856 

16261 

20.3 

17.29 

0.8518 

2 . 37 

19596 

21297 

20.3 

18.68 

0.9200 

5.36 

30984 

32060 

20.3 

19.63 

0.9671 

10.37 

43818 

44578 

20.3 

19.96 

0.9833 

Note:  Table  3.3  is  calculated  using  the  developed 
approach,  from  thermal  data  obtained  from  Perkins  [18] 
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Table  3.4  - Effective  Scattering  Cross  Section  for 
Hydrogen  at  2 500  °K,  with  V„=10^  m/s 
(completely  entrained  neutrons) . 


En 

(eV) 

(m/s) 

V* 

(m/s) 

^%ff 

(Vm) 

^%ff 

(V*), 

V„/V* 

0.53 

1386 

10064 

110.0 

15.15 

0.1377 

0.54 

2400 

10192 

65.0 

15.31 

0.2355 

0.56 

3098 

10320 

52.0 

15.61 

0.3002 

0.59 

4382 

10640 

40.0 

16.48 

0.4119 

0.74 

7589 

11920 

28.0 

17.82 

0.6365 

0.91 

9798 

13194 

25.0 

18.56 

0.7424 

1.38 

13856 

16261 

22.0 

18.74 

0.8518 

2.37 

19596 

21297 

20.5 

18.86 

0.9200 

5.36 

30984 

32060 

20.4 

19.73 

0.9671 

10.37 

43818 

44578 

20.3 

19.96 

0.9833 

Note:  Table  3.4  is  calculated  using  the  developed 
approach,  from  thermal  data  obtained  from  Perkins  [18]. 
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From  Table  3.3,  it  can  be  observed  that  at  0 °K  the 
values  of  the  effective  cross  section  in  the  laboratory  frame 
are  initially  very  small  and  increase  as  the  neutron  speed 
increases,  tending  to  the  asymptotic  value  of  20.3  barns.  This 
can  be  explained  if  it  is  recognized  that  once  the  neutrons 
are  fully  entrained  they  acquire  dynamical  equilibrium  with 
the  system  and  then,  for  neutron  speeds  comparable  to  the 
speed  of  the  medium,  the  relative  velocity  and  interaction 
probability  are  very  small  (in  the  limit  of  zero  relative 
velocity,  no  reaction  occurs).  As  the  neutron  energy 
increases,  the  relative  velocity  also  increases  and  an 
increase  in  the  effective  cross  section  can  be  observed.  When 
the  neutron  velocity  is  much  higher  than  the  velocity  of  the 
medium,  v„  is  approximately  equal  to  v*  and  the  actual 
effective  cross  section  in  the  laboratory  system  is 
approximately  equal  to  the  effective  cross  section  as  a 
function  of  only  the  neutron  velocity  in  the  laboratory  system 
only. 

Table  3.4  presents  similar  results  but  for  hydrogen  at 
2500  °K  instead  of  at  0 °K.  Comparing  Table  3.3  and  Table  3.4, 
it  can  be  observed  that  the  addition  of  the  thermal  movement 
of  the  medium  has  the  effect  of  increasing  the  effective  cross 
sections  for  neutrons  with  lower  velocities.  As  the  velocity 
of  the  neutrons  increases,  this  effect  tends  to  be  reduced 
since  the  neutron-nucleus  relative  velocity  tends  towards  the 
value  of  the  neutron  velocity,  as  expected.  Physically,  the 
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neutron  is  still  entrained  with  the  directed  velocity  but  the 
temperature  component  of  the  target  velocity  leads  to  a higher 
probability  of  interaction. 

For  very  high  neutron  velocities.  Tables  3.3  and  3.4 
predict  the  same  values  for  the  effective  cross  section  in  the 
static  frame.  The  results  presented  in  Table  3.3  and  Table 

3.4  are  graphically  shown  in  Figures  3.9  and  3.10, 
respectively. 


3-^ — Adaptations  in  the  BRT-I  Computer  Code 

By  making  the  assumptions  mentioned  in  Section  3.5  it  is 
now  possible  to  use  the  computer  code  BRT-I  to  solve  the 
problem  for  selected  conditions  in  the  moving  frame  and  to 
then  convert  the  obtained  results  back  to  the  static  frame. 

To  use  BRT-I  for  this  purpose,  an  option  of  the  code  is 
invoked  which  permits  the  manual  input  of  the  cross  section 
and  the  energy  or  speed  mesh  data.  Since  the  code  will  be  used 
to  solve  the  integral  transport  equation  in  the  moving  frame, 
the  velocity  (speed)  input  must  be  the  velocity  (speed)  of  the 
neutrons  in  the  moving  frame,  v^  for  the  case  of  uncollided 
neutrons  or  for  the  case  where  the  neutrons  are  completely 
entrained.  Consequently,  the  obtained  results  are  for  the 
moving  frame  and  there  is  thus  a need  to  modify  the  code  in 
order  to  obtain  results  related  to  the  speed  of  the  neutron  in 
the  static  frame,  v,  for  the  completely  entrained  case. 


iff  (v  ) (barn) 


76 


Figure  3.9  Effective  Scattering  Cross  Section  for 
Hydrogen  at  0°K  (completely  entrained 
neutrons) . 


(V  ) (barn) 
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Figure  3.10  - Effective  Scattering  Cross  Section  for 

Hydrogen  at  2500°K  (completely  entrained 
neutrons) . 
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According  to  the  method  selected  to  deal  with  the 
problem,  the  BRT-I  computer  code  will  have  as  input  data 
related  to  the  moving  frame  (v^  or  v„)  . The  code  will  be 
modified  by  inserting  a new  subroutine  DIRVE,  for  the  case  of 
completely  entrained  neutrons.  This  routine  will  be  invoked 
after  the  solution  of  the  integral  equation  is  completed,  but 
before  subroutine  EDIT  is  called,  in  order  to  have  the  output 
related  not  only  to  the  moving  frame,  but  also  to  the  static 
frame.  The  new  structure  of  the  code,  with  the  inclusion  of 

DIRVE,  and  the  source  listing  of  the  new  subroutine  is  shown 
in  Appendix  D. 

Since  the  problem  now  includes  the  directed  velocity  of 
the  target  atoms,  the  code  also  has  to  be  modified  to  include 
the  speed  of  the  medium,  V„,  as  an  input.  Although  the 
adaptations  made  in  the  computer  code  structure  are  small,  the 
inputs  must  be  carefully  selected,  particularly  the  input  of 
the  speed  mesh. 


CHAPTER  4 

EFFECTS  OF  HIGH  TARGET  ATOM  DIRECTED 
VELOCITY  AND  TEMPERATURE 


4 . 1 Introduction 

The  purpose  of  this  Chapter  is  to  present  and  discuss  the 
effects  resulting  from  analysis  of  selected  combinations  of 
temperature  and  directed  velocity  of  target  atoms  using  the 
BRT-PCM  code  (see  Appendix  D)  . The  results  for  the  completely 
entrained  neutron  condition  in  the  static  frame  are  obtained 
from  BRT-PCM.  The  results  for  the  uncollided  neutron  condition 
in  the  moving  frame  are  obtained  from  BRT-PCM  and  then 
converted  to  the  static  frame  using  the  approximate  equations 
derived  in  Chapter  3 for  this  condition. 

The  cases  examined  include  combinations  of  temperature 
and  directed  velocity  of  the  target  atoms  which  are  presented 
in  Table  4.1. 


Table  4.1  - Selected  Target  Atoms  Temperatures 
and  Directed  Velocities. 


T (°K) 

2000 

3000 

4000 

5000 

7000 

0.0 

1.0 

2 . 0 

3.0 

4.0 

79 
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It  should  be  mentioned  that  in  Table  4.1,  and  for  all 
cases  in  this  Chapter,  the  target  atom  velocity  is  expressed 
in  units  of  2200  m/s. 

All  cases  studied  use  the  material  composition  shown  in 
Table  4.2  and  the  system  modeled  is  similar  in  geometry  to  the 
one  presented  in  Appendix  B,  Section  B.l. 


Table  4.2  - Composition  and  Number 
Density  of  Materials. 


NUCLIDE 

N 

(atoms/b.cm) 

Core 

nitrogen 

1.000  E-12 

hydrogen 

1.287  E-3 

uranium  235 

4.283  E-5 

uranium  238 

3.224  E-5 

helium 

6.764  E-4 

fluorine 

2.763  E-4 

Moderator 

beryllium 

1.200  E-1 

Although  the  choice  of  temperatures  and  directed 
velocities  of  the  target  atoms  is  influenced  by  the  problem  of 
the  UTVR/MHD  Disk  Generator  (see  Chapter  1),  the  conclusions 
related  to  the  effects  of  the  directed  velocity  of  the  target 
atoms  are  not  tied  to  any  particular  system  and  should  be 
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valid  for  different  configurations,  possessing  similar  ranges 
of  temperatures  and  directed  velocities  of  the  target  atoms. 

BRT-  PCM  was  used  with  20  space  points  and  speed  points 
for  all  cases.  The  lower  limit  of  the  speed  structure  was 
selected  according  to  each  particular  case  studied  and  the 
upper  limit  was  set  to  be  approximately  2.0  eV. 

4_,2  Results  Obtained  Using  BRT-PCM  for  Selected  Cases  of 

High  Target  Atom  Directed  Velocity  and  Temperature 

According  to  the  procedures  described  in  the  preceding 
Chapter,  the  analysis  is  performed  separately  for  the 
conditions  of  completely  entrained  neutrons  and  uncollided 
neutrons.  Results  examined  include  the  averaged  macroscopic 
neutron  scattering,  absorption,  and  fission  cross  sections  for 
the  uncollided  condition.  For  the  completely  entrained 
condition  the  results  examined  include  the  neutron  flux 
distribution,  the  averaged  macroscopic  neutron  scattering, 
absorption,  and  fission  cross  sections,  and  also  the  inverse 
of  the  average  neutron  velocity. 

A graphical  representation  of  the  neutron  flux 
distribution  as  a function  of  energy,  for  zero  target  atom 
directed  velocity  and  for  target  atom  and  neutron  source 
temperatures  of  2 000  °K,  4000  °K  and  7000  °K  is  shown  in  Figure 

4.1. 
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Figure  4.1  - Plot  of  Flux  (E)  Versus  E for  V = 0.0  (neutron  source 
temperature  equals  gas  temperature) . 
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4.2.1  BRT-PCM  Results  for  the  Combined  Effects  of  High 

Directed  Velocity  and  Temperature  of  the  Target  Atoms  for 
Case  of  Uncollided  Neutrons 

The  results  for  the  averaged  macroscopic  cross  sections 
are  obtained  using  the  data  from  BRT-PCM  for  zero  target  atom 
directed  velocity.  Since  the  case  examined  is  for  uncollided 
neutrons,  all  the  isotopes  are  at  the  same  temperature. 

The  approach  to  determine  the  averaged  macroscopic  cross 
sections  is  divided  in  steps  as  follows:  l)  the  average 
neutron  velocity  for  each  system  temperature,  v,  is  obtained 
from  BRT-PCM  (the  neutron  flux  distributions  for  the 
uncollided  cases  are  those  given  in  Figure  4.1);  2)  according 
to  the  approach  developed  in  Chapter  3 , the  speed 
corresponding  to  the  average  neutron  velocity,  v,  and  the 
target  atom  directed  velocity,  V„,  are  used  to  determine  v^ 
(see  Equations  3.2a  and  3.2b)  and  the  scaling  factor  v^/v;  3) 
a procedure  similar  to  the  one  presented  in  Chapter  3,  Section 
3.5,  to  determine  the  hydrogen  effective  microscopic 
scattering  cross  section,  including  the  target  atom  directed 
velocity,  is  used  to  determine,  for  each  isotope,  first  the 
effective  microscopic  cross  sections  at  v^  (thermal  motion 
only),  cjg^^(Vg),  using  the  raw  data  from  the  NITAWL  library;  4) 
then,  using  Equation  3.5,  the  effective  microscopic  cross 
sections  at  v,  including  the  directed  velocity  effects  is 
determined  (see,  for  example,  the  results  for  hydrogen  in 
Tables  3.1  and  3.2);  and  5)  the  effective  macroscopic  cross 
sections,  including  the  target  atom  directed  velocity. 


are 
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obtained  using  the  atom  number  densities  and  the  effective 
microscopic  cross  sections  determined  in  step  4. 

To  examine  the  target  atom  directed  velocity  effects,  all 
the  results  are  presented  for  the  average  neutron  velocity 
corresponding  to  the  system  temperature. 


Table  4.3  - Inverse  of  the  Average  Neutron  Velocity  as  a 
Function  of  Target  Atom  Temperature  (neutron  source 
temperature  equals  gas  temperature) . 


2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

o 

o 

0.3187 

0.2635 

0.2322 

0.2128 

0.1912 

'i'hs  results  for  the  "uncollided”  cases  are  also  presented 
for  the  condition  when  the  neutron  source  is  kept  fixed  at 
2000  °K  while  the  target  atom  or  gas  temperature  takes  on 
values  presented  in  Table  4.1.  The  reason  for  this  is  to 
provide  an  alternate  comparison  of  the  results  obtained  for 
the  uncollided  neutron  cases  with  the  results  obtained  for  the 
completely  entrained  neutron  cases,  where  the  source 
temperature  is  also  kept  fixed  at  2 000  °K  while  the  target 
atom  or  gas  temperature  takes  on  values  presented  in  Table 
4.1.  The  inverse  of  the  average  neutron  velocity  used  in  the 
computations  for  these  cases  is  shown  in  Table  4.4. 
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Table  4.4  - Inverse  of  the  Average  Neutron  Velocity  as  a 
Function  of  Target  Atom  Temperature  (neutron  source 
temperature  at  2000  °K)  . 


Vm 

2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

o 

o 

0.3187 

0.2924 

0.2742 

0.2613 

0.2450 

4. 2. 1.1  Averaged  macroscopic  neutron  scattering  cross  section 
The  results  for  the  averaged  macroscopic  scattering  cross 
section  are  presented  in  Table  4.5  and  Figure  4.2. 

The  results  in  Table  4.5  show  that  for  a fixed  directed 
velocity  of  the  medium,  as  the  temperature  increases,  the 
cross  section  decreases,  as  expected.  This  is  because 
hydrogen,  the  most  significant  isotope,  as  far  as  scattering 
is  concerned,  presents  a scattering  cross  section  which 
decreases  with  increasing  energy  in  the  low  energy  region  (see 
Chapter  3,  Section  3.5). 

Table  4.5  - Averaged  Macroscopic  Neutron  Scattering  Cross 

Section  (1/cm)  as  a Function  of  Target  Atom  Temperature 
and  Directed  Velocity  (uncollided  neutrons:  neutron 
source  temperature  equals  gas  temperature) . 


Vm 

2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

o 

• 

o 

3.278  E-2 

3.008  E-2 

2.800  E-2 

2.644  E-2 

2.429  E-2 

1.0 

3.371  E-2 

3.063  E-2 

2.836  E-2 

2.671  E-2 

2.447  E-2 

2.0 

3.644  E-2 

3.221  E-2 

2.942  E-2 

2.749  E-2 

2.498  E-2 

3.0 

4.067  E-2 

3.466  E-2 

3.105  E-2 

2.869  E-2 

2.576  E-2 

4.0 

4.572  E-2 

3.765  E-2 

3.304  E-2 

3.013  E-2 

2.670  E-2 
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For  all  the  cases  it  can  be  observed  that  the  averaged 
macroscopic  scattering  cross  section  increases  continuously 
with  increasing  V^. 

From  the  table  it  can  be  observed  that  the  ratio  between 
the  cross  sections  at  Vj^=4 . 0 and  Vj^=0.0  decreases  continuously 
with  increasing  temperature  (1.395  at  2000  °K  and  1.099  at 
7000  °K)  . This  demonstrates  that  as  the  temperature  increases 
the  directed  velocity  effects  are  less  pronounced.  The  ratio 
between  the  averaged  macroscopic  scattering  cross  section  at 
2 000  °K  and  7 000  °K  increases  continuously  with  increasing 
This  is  because  at  higher  temperatures  the  average  neutron 
velocities  are  also  higher  and  the  effect  of  increasing  target 
atom  directed  velocity  is  less  pronounced  than  at  lower 
temperatures  which  present  lower  average  neutron  velocities 
(see  Chapter  3,  Table  3.1).  These  effects  are  clearly  shown  in 
Figure  4.2. 

The  continuous  increase  in  the  averaged  macroscopic 
neutron  scattering  cross  section  with  increasing  directed 
velocity  can  be  explained  if  it  is  remembered  from  Chapter  3, 
Section  3.5,  that  the  directed  velocity  of  the  target  atoms 
has  the  effect  of  increasing  the  hydrogen  effective 
microscopic  scattering  cross  section,  the  most  significant 
isotope,  as  far  as  scattering  is  concerned,  for  the  case  of 
uncollided  neutrons.  The  effects  of  the  target  atom  directed 
velocity  on  the  other  isotopes  depend  upon  the  behavior  of  the 


( 1/cm) 
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Figure  4.2  - Plot  of  Averaged  Macroscopic 

Scattering  Cross  Section  Versus 
Temperature  (uncollided  neutrons) . 
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respective  effective  microscopic  cross  sections  with 
increasing  energy  (directed  velocity) . 

Table  4.6  shows  the  results  for  the  averaged  macroscopic 
scattering  cross  section  when  the  neutron  source  is  kept  fixed 
at  2000  °K. 


Table  4.6-  Averaged  Macroscopic  Neutron  Scattering  Cross 

Section  (1/cm)  as  a Function  of  Target  Atom  Temperature 
and  Directed  Velocity  (uncollided  neutrons:  neutron 
source  temperature  at  2 000  °K)  . 


2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

O 

• 

o 

3.278  E-2 

3.242  E-2 

3.186  E-2 

3.128  E-2 

3.023  E-2 

1.0 

3.371  E-2 

3.318  E-2 

3.250  E-2 

3.194  E-2 

3.069  E-2 

2 . 0 

3.644  E-2 

3.539  E-2 

3.436  E-2 

3.345  E-2 

3.201  E-2 

3.0 

4.067  E-2 

3.880  E-2 

3.722  E-2 

3.594  E-2 

3.404  E-2 

4.0 

4.572  E-2 

4.295  E-2 

4.072  E-2 

3.899  E-2 

3.653  E-2 

4. 2. 1.2  Averaged  macroscopic  neutron  absorption  cross  section 
The  results  obtained  for  this  variable  are  presented  in 
4 ^ for  the  various  combinations  of  directed  velocity 
and  temperature  of  the  target  atoms. 

Table  4 . 7 shows  that  the  averaged  macroscopic  neutron 
absorption  cross  section  tends  to  decrease  as  the  directed 
velocity  of  the  medium  increases,  an  effect  similar  to  the  one 
observed  for  a temperature  increase  only.  This  reflects  the 
energy  dependent  behavior  of  the  effective  microscopic 
absorption  cross  section  for  all  the  involved  isotopes. 
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Table  4.7  - Averaged  Macroscopic  Neutron  Absorption  Cross 

Section  (1/cm)  as  a Function  of  Target  Atom  Temperature 
and  Directed  Velocity  (uncollided  neutrons:  neutron 
source  temperature  equals  gas  temperature) . 


Vm 

2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

o 

• 

o 

8.066  E-3 

6.133  E-3 

4.988  E-3 

4.267  E-3 

3.458  E-3 

1.0 

8.029  E-3 

6.010  E-3 

4.864  E-3 

4.167  E-3 

3.403  E-3 

to 

• 

o 

7.761  E-3 

5.597  E-3 

4.499  E-3 

3.890  E-3 

3.260  E-3 

3.0 

6.885  E-3 

4.854  E-3 

3.963  E-3 

3.527  E-3 

3.098  E-3 

4.0 

5.434  E-3 

3.998  E-3 

3.463  E-3 

3.235  E-3 

2.967  E-3 

It  can  be  observed  that  as  the  temperature  increases,  the 
effect  of  the  directed  velocity  of  the  medium  decreases  (the 
ratio  between  the  cross  sections  at  V^=0.0  and  Vj^=4.0  is  1.484 
at  2000  °K,  and  1.165  at  7000  °K)  . 

The  presented  changes  in  cross  sections  with  temperature 
also  show  that,  as  the  velocity  of  the  medium  is  increased, 
the  temperature  effect  is  less  pronounced. 

The  effects  presented  for  the  averaged  macroscopic 
absorption  cross  section  show  that  if  the  directed  velocity  of 
the  medium  is  kept  constant  and  the  temperature  of  the  medium 
is  continuously  increased,  a point  will  be  reached  where  the 
effects  of  the  directed  velocity  of  the  medium  are  so  small, 
compared  to  the  effects  of  temperature,  that  they  can  be 
neglected.  The  same  reasoning  applies  for  the  case  when  the 
temperature  of  the  medium  is  kept  constant  and  the  directed 
velocity  of  the  target  atoms  is  increased. 
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Table  4.8  shows  the  results  for  the  averaged  macroscopic 
absorption  cross  section  when  the  neutron  source  is  kept  at 
2000  °K. 


Table  4.8  - Averaged  Macroscopic  Neutron  Absorption  Cross 

Section  (1/cm)  as  a Function  of  Target  Atom  Temperature 
and  Directed  Velocity  (uncollided  neutrons:  neutron 
source  temperature  at  2000  °K)  . 


V. 

2000 

°K 

3000 

°K 

4000 

°K 

5000 

°K 

7000 

°K 

o 

o 

8 . 066 

E-3 

7.172 

E-3 

6.528 

E-3 

6 

.060 

E-3 

5 

.464 

E-3 

1.0 

8.029 

E-3 

7.084 

E-3 

6.415 

E-3 

5 

.935 

E-3 

5 

.335 

E-3 

2 . 0 

7.761 

E-3 

6.715 

E-3 

6.008 

E-3 

5 

.521 

E-3 

4 

.937 

E-3 

3 . 0 

6.885 

E-3 

5.868 

E-3 

5.217 

E-3 

4 

.790 

E-3 

4 

.303 

E-3 

4 . 0 

5.434 

E-3 

4.689 

E-3 

4.239 

E-3 

3 

.957 

E-3 

3 

. 657 

E-3 

The  results  from  Table  4 . 7 are  represented  graphically  in 
Figure  4.3. 

4. 2. 1.3  Averaged  macroscopic  neutron  fission  cross  section 

The  results  for  the  averaged  macroscopic  fission  cross 
section  are  presented  in  Table  4.9. 

The  results  in  Table  4.9  present  a behavior  very  similar 
to  that  of  Table  4.7;  the  averaged  macroscopic  fission  cross 
section  decreases  continuously  with  increasing  target  atom 
directed  velocity. 


( 1/cm) 
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9.0  E-3 


Figure  4.3  - Plot  of  Averaged  Macroscopic 

Absorption  Cross  Section  Versus 
Temperature  (uncollided  neutrons). 
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Table  4.9  - Averaged  Macroscopic  Neutron  Fission  Cross 

Section  (1/cm)  as  a Function  of  Target  Atom  Temperature 
and  Directed  Velocity  (uncollided  neutrons:  neutron 
source  temperature  equals  gas  temperature) . 


.Vm...  .1 

2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

o 

o 

6.568  E-3 

4.997  E-3 

4.070  E-3 

3.493  E-3 

2.857  E-3 

1.0 

6.583  E-3 

4.914  E-3 

3.964  E-3 

3.392  E-3 

2.785  E-3 

2 . 0 

6.488  E-3 

4.604  E-3 

3.633  E-3 

3.099  E-3 

2.591  E-3 

3 . 0 

5.876  E-3 

3.959  E-3 

3.091  E-3 

2.676  E-3 

2.357  E-3 

o 

• 

4.524  E-3 

3.062  E-3 

2.497  E-3 

2.289  E-3 

2.209  E-3 

The  directed  velocity  effects  decrease  with  increasing 
temperature,  as  expected.  The  explanations  presented  for  Item 
4.2. 1.2  above  also  apply  for  this  case.  However,  the  effect  of 
the  directed  velocity  of  the  target  atoms  in  decreasing  the 
averaged  macroscopic  neutron  fission  cross  section  is  slightly 
larger  than  for  the  averaged  macroscopic  neutron  absorption 
cross  section.  This  is  because  the  microscopic  fission  cross 
section  for  falls  off  faster  with  energy  than  the 
microscopic  absorption  cross  section  for  the  same  isotope. 

Results  for  the  averaged  macroscopic  fission  cross 
section  when  the  neutron  source  is  kept  constant  at  2000  °K 
are  presented  in  Table  4.10. 
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Table  4.10  - Averaged  Macroscopic  Neutron  Fission  Cross 

Section  (1/cm)  as  a Function  of  Target  Atom  Temperature 
and  Directed  Velocity  (uncollided  neutrons:  neutron 
source  temperature  at  2000  °K)  . 


Vm 

2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

0.0 

6.568  E-3 

5.842  E-3 

5.318  E-3 

4.940  E-3 

4.461  E-3 

1.0 

6.583  E-3 

5.804  E-3 

5.249  E-3 

4.854  E-3 

4.360  E-3 

2 . 0 

6.488  E-3 

5.583  E-3 

4.965  E-3 

4.540  E-3 

4.026  E-3 

3 . 0 

5.876  E-3 

4.924  E-3 

4.306  E-3 

3.898  E-3 

3.428  E-3 

4.0 

4.524  E-3 

3.774  E-3 

3.311  E-3 

3.020  E-3 

2.704  E-3 

4.2.2  BRT-PCM  Results  for  the  Combined  Effects  of  High 
Directed  Velocity  and  Temperature  of  the  Target  Atoms  for 
the  Case  of  Completely  Entrained  Neutrons 

Since  the  main  focus  of  this  work  is  the  analysis  of  the 
gaseous  fuel  region  (see  Chapter  2,  Section  2.1),  the  neutron 
source  temperature  is  kept  constant  at  2 000  °K  in  the  analysis 
for  the  completely  entrained  neutron  condition.  This  is 
because  physical  limitations  restrict  the  external  moderator/ 
reflector  region  (which  is  likely  to  consist  of  BeO  and/or 
C^^)  , where  the  neutron  source  temperature  is  established,  to 
around  2000  °K.  Note  that  such  a limitation  should  also  be 
imposed  for  the  uncollided  condition.  However,  if  this  is 
done,  then  the  assumptions  upon  which  the  developed  solution 
method  is  based  would  be  violated.  Specifically,  in  using  BRT- 
PCM  to  determine  interaction  rates  it  is  assumed  that  the 
neutrons  and  target  atoms  are  in  thermal  equilibrium.  If  the 
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moderator/reflector  region  temperature  is  fixed  at  2000  °K, 
then,  for  higher  gas  temperatures,  the  uncollided  neutron  flux 
and  gas  target  atoms  would  no  longer  be  in  thermal 
equilibrium.  Thermal  equilibrium  in  the  uncollided  case  is 
only  possible  if  both  the  moderator/reflector  and  gas 
temperature  are  the  same. 

A graphical  representation  of  the  neutron  flux 
distribution  as  a function  of  energy,  for  zero  target  atom 
directed  velocity  and  for  target  atom  temperatures  of  2000  °K, 
4000  °K  and  7000  °K  is  shown  in  Figure  4.4;  the  neutron  source 
temperature  is  2000  °K. 

4 . 2 . 2 . 1 Neutron  flux 

The  results  obtained  from  BRT-PCM  for  the  neutron  flux 
distribution,  as  a function  of  energy,  for  target  atom  (gas) 
temperatures  of  2000  °K,  4000  °K,and  7000  °K,  respectively, 
with  the  directed  velocity  of  the  target  atoms  varied 
parametrically  are  graphically  represented  in  Figures  4.5,  4.6 
and  4.7.  The  neutron  source  temperature  in  all  cases  is  kept 
fixed  at  2000  °K. 

The  purpose  of  these  figures  is  to  show  the  general 
behavior  of  the  neutron  flux  for  the  cases  when  the  target 
atoms  are  subjected  to  the  effects  of  both  high  temperature 
and  directed  velocity  and  to  compare  the  results  for  different 
combinations  of  temperature  and  directed  velocity.  The 
selected  combinations  of  temperatures  and  directed  velocities 
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Figure  4.4  - Plot  of  Flux  (E)  Versus  E for 
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Figure  4.5  - Plot  of  Flux  (E)  Versus  E for  Tg^g  = 2000°K  (com 
pletely  entrained  neutrons) . 
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Figure  4.6  - Plot  of  Flux  (E)  Versus  E for  Tg^g  = 4000°K  (completely 
entrained  neutrons) . 
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Figure  4.7  - Plot  of  Flux  (E)  Versus  E for  Tg^g  = VOOQOk  (completely 
entrined  neutrons) . 
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presented  in  these  figures  is  made  in  order  to  clearly  show 
the  difference  in  the  effects. 

All  the  figures  show  the  same  general  behavior:  the 
addition  of  the  directed  velocity  of  the  target  atoms  shifts 
the  neutron  flux  to  higher  energies,  as  expected.  As  the 
directed  velocity  of  the  target  atoms  is  increased  this  effect 
becomes  more  pronounced  (see,  e.g..  Fig  4.5).  It  can  be 
observed  that  the  inclusion  of  the  directed  velocity  results 
in  a shifting  of  the  low  energy  end  of  the  curve  to  an  energy 
that  is  much  higher  than  the  low  energy  end  for  the 
corresponding  case  which  includes  only  temperature  effects.  It 
can  also  be  observed  that  the  maximum  of  the  neutron  flux 
distribution  curve  increases  with  increasing  target  atom 
directed  velocity  and  is  also  shifted  to  a higher  energy  (see 
Fig.  4.5) . 

The  behavior  of  the  neutron  flux,  in  the  presence  of  the 
directed  velocity  effect  of  the  target  atoms,  which  shows 
values  in  the  low  energy  end  of  the  distribution  that  are  much 
smaller,  at  the  same  energies,  than  the  corresponding  values 
when  only  temperature  effects  are  included,  helps  to  explain 
the  behavior  of  the  averaged  macroscopic  cross  sections. 

Figures  4.5,  4.6  and  4.7  show  that  for  sufficiently  high 
neutron  energies  the  effect  of  the  directed  velocity  of  the 
target  atoms  is  small  and  the  distribution  curve  exhibits 
essentially  the  shape  it  would  have  for  temperature  effects 
only,  with  only  a small  deviation  in  magnitude.  It  should  also 
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be  noticed  that  for  a specific  directed  velocity  of  the  target 
atoms,  the  energy  corresponding  to  the  low  energy  end  of  the 
curves  is  always  the  same,  independent  of  the  medium's 
temperature  (see  Figures  4.5,  4.6,  and  4.7). 

The  behavior  presented  for  the  neutron  flux  with  the 
inclusion  of  the  directed  velocity  of  the  target  atoms  is 
consistent  with  the  approach  developed  in  Chapter  3 and  also 
makes  sense  physically  if  it  is  realized  that  neutrons  in  the 
static  frame  initially  with  energies  (velocities)  smaller  than 
the  directed  velocity  of  the  target  atoms  will  have  to  suffer 
upscatterings  to  become  completely  entrained  and  these 
upscatterings  will  induce  a hardening  of  the  neutron  flux 
distribution.  As  the  energy  of  the  neutrons  increases,  the 
effect  of  the  directed  velocity  on  the  effective  microscopic 
cross  sections  and  on  the  neutron  number  density  distribution 
is  diminished;  the  upscattering  is  less  pronounced  than  for 
small  neutron  energies  and  this  explains  the  behavior  of  the 
neutron  flux  at  high  energies. 

4. 2. 2. 2 Averaged  macroscopic  neutron  scattering  cross  section 

The  results  for  the  averaged  macroscopic  scattering  cross 
section  obtained  from  BRT-PCM  are  presented  in  Table  4.11. 

Table  4.11  demonstrates  that  for  the  case  of  completely 
entrained  neutrons,  for  a fixed  temperature,  an  increase  in 
the  directed  velocity  of  the  target  atoms  yields  a decrease  in 
the  averaged  macroscopic  neutron  scattering  cross  sections. 
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The  results  presented  in  Table  4 . 11  are  shown  graphically 
in  Figure  4.8. 


Table  4.11  - Averaged  Macroscopic  Neutron  Scattering  Cross 
Section  (1/cm)  as  a Function  of  Target  Atom  Temperature 
and  Directed  Velocity  (completely  entrained  neutrons) . 


Vm 

2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

o 

• 

o 

3.277  E-2 

3.242  E-2 

3.186  E-2 

3.128  E-2 

3.023  E-2 

1.0 

3.101  E-2 

3.091  E-2 

3.053  E-2 

3.008  E-2 

2.922  E-2 

2 . 0 

2.749  E-2 

2.773  E-2 

2.760  E-2 

2.734  E-2 

2.673  E-2 

3.0 

2.392  E-2 

2.437  E-2 

2.443  E-2 

2.431  E-2 

2.391  E-2 

4 . 0 

2.076  E-2 

2.135  E-2 

2.153  E-2 

2.152  E-2 

2.128  E-2 

When  the  increase  in  temperature,  for  a fixed  directed 
velocity  of  the  target  atoms,  is  examined,  it  is  observed 
that:  1)  at  V„=0.0,  the  averaged  macroscopic  neutron 
scattering  cross  section  continuously  decreases,  as  should  be 
expected  for  temperature  effects  only  (see  the  discussion 
under  Item  4. 2. 1.1;  2)  at  V|^=1.0,  the  averaged  macroscopic 
neutron  scattering  cross  sections  are  smaller  than  for  Vj^=0.0, 
due  to  the  flux  hardening  effect  of  the  directed  velocity  of 
the  target  atoms  and,  again,  continuously  decrease  with  an 
increase  in  temperature,  but  at  a lower  rate  than  when  V„=0; 
3)  at  V|^=2.0,  the  averaged  macroscopic  cross  sections  have 
smaller  values  than  at  V^=1.0,  as  expected,  but  they  increase 
with  temperature  up  to  T=3000  °K  and  then  decrease 
continuously;  4)  at  V^=3 . 0 , the  averaged  macroscopic  cross 


( 1/cm) 
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- Plot  of  Averaged  Macroscopic 
Scattering  Cross  Section  Versus 
Temperature  (completely  entrained 
neutrons) . 


Figure  4.8 
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sections  increase  with  temperature  until  a maximum  at  4000  °K 
and  then  decrease;  5)  at  V^=4 . 0 , the  averaged  cross  sections 
are  smaller  than  at  Vj^=3 . 0 , and,  again,  increase  with 
temperature  until  4000  °K  and  then  decrease. 

The  above  demonstrates  that  as  the  directed  velocity 
increases,  the  tendency  of  continuously  decreasing  averaged 
macroscopic  cross  sections  with  increasing  temperatures  is 
altered  for  values  of  2 or  greater  (see  Fig.  4.8);  for 
these  values,  as  the  temperature  is  increased  the  averaged 
cross  sections  first  tend  to  increase,  reaching  a maximum  and 
then,  with  a further  increase  in  temperature,  they  decrease. 
This  behavior  can  be  explained  if  it  is  understood  that 
although  an  increase  in  either  the  directed  velocity  or  the 
temperature  of  the  target  atoms  contributes  to  a decrease  in 
the  averaged  macroscopic  cross  section,  the  mechanisms 
involved  are  not  exactly  the  same.  The  increase  in  temperature 
decreases  the  averaged  macroscopic  cross  sections  by  hardening 
the  neutron  spectrum,  as  explained  before.  The  directed 
velocity  of  the  target  atoms  decreases  the  effective 
microscopic  (macroscopic)  cross  section,  particularly  in  the 
low  energy  region,  and  also  hardens  the  neutron  spectrum,  as 
explained  before,  by  shifting  the  low  energy  end  of  the 
distribution  curve  to  a higher  energy  (which  contributes  to  a 
further  decrease  in  the  averaged  macroscopic  cross  section) . 
As  the  directed  velocity  of  the  target  atoms  increases,  the 
effective  microscopic  cross  section  decreases  at  low  energies. 
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As  the  temperature  increases,  for  a especific  directed 
velocity,  the  flux  tends  to  increase  at  high  energies  and 
decrease  at  low  energies.  Hence,  for  sufficiently  low 
temperatures  and  for  sufficiently  high  directed  velocities, 
the  combination  of  the  comparatively  high  flux  at  low  energies 
where  the  effective  microscopic  cross  sections  are  very  small 
results  in  an  averaged  macroscopic  cross  section  that  is 
smaller  than  the  averaged  macroscopic  cross  section  expected 
for  temperature  effects  only.  As  the  temperature  increases, 
for  the  same  directed  velocity  of  the  target  atoms,  the  flux 
tends  to  shift  towards  high  energies  and  the  weighed 
contribution  of  the  reduced  low  energy  effective  microscopic 
cross  sections  is  smaller  than  before  (at  lower  temperature) , 
leading  to  an  increase  in  the  averaged  macroscopic  cross 
section,  until  a temperature  is  reached  where  the  averaged 
macroscopic  cross  section  attains  a maximum  value;  after  this 
point,  the  spectral  hardening  due  to  the  temperature  effect 
causes  the  averaged  macroscopic  cross  section  to  decrease 
again.  The  temperature  corresponding  to  this  maximum  can  be 
considered  as  the  temperature  at  which  the  thermal  and 
directed  velocity  induced  hardening  effects  become  more 
important  than  the  target  atom  directed  velocity  effect  of 
decreasing  the  effective  microscopic  cross  section.  With 
further  increases  in  temperature,  the  spectral  hardening 
continues  to  cause  the  averaged  macroscopic  cross  section  to 


decrease. 
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The  above  described  phenomena  associated  with  the 
directed  velocity  of  the  target  atoms  and  their  effects  on  the 
averaged  macroscopic  scattering  cross  section  apply  for  all 
the  other  averaged  macroscopic  cross  sections,  i.e.,  for 
absorption  and  fission  as  well  as  for  scattering. 

4 . 2 . 2 . 3 Averaged  macroscopic  neutron  absorption  cross  section 
The  results  obtained  for  this  variable  are  presented  in 
Table  4 . 12  for  the  various  combinations  of  directed  velocity 
and  temperature  of  the  target  atoms. 

Table  4 . 12  shows  for  the  averaged  macroscopic  neutron 
absorption  cross  section,  for  completely  entrained  neutrons, 
a behavior  similar  to  the  one  presented  in  Table  4.7  for  the 
uncollided  neutrons,  although  the  details  of  the  cross  section 
behavior  as  dictated  by  the  phenomena  associated  with  the 
directed  velocity  are  different  in  the  two  cases. 


Table  4.12  - Averaged  Macroscopic  Neutron  Absorption  Cross 
Section  (1/cm)  as  a Function  of  Target  Atom  Temperature 
and  Directed  Velocity  (completely  entrained  neutrons) . 


Vm 

2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

0.0 

8.066  E-3 

7.172  E-3 

6.528  E-3 

6.060  E-3 

5.464  E-3 

1.0 

7.525  E-3 

6.759  E-3 

6.193  E-3 

5.777  E-3 

5.247  E-3 

2.0 

6.485  E-3 

5.913  E-3 

5.467  E-3 

5.131  E-3 

4.694  E-3 

3 . 0 

5.475  E-3 

5.046  E-3 

4.697  E-3 

4.429  E-3 

4.076  E-3 

4.0 

4.645  E-3 

4.311  E-3 

4.031  E-3 

3.813  E-3 

3.523  E-3 
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Table  4 . 12  shows  that  the  averaged  macroscopic  neutron 
absorption  cross  section  decreases  as  the  directed  velocity  of 
the  medium  or  the  temperature  increases.  These  are  the  general 
trends  which  are  expected  as  explained  in  Item  4. 2. 2. 2. 

As  for  the  case  of  uncollided  neutrons,  it  can  be 
observed  that  as  the  temperature  increases,  the  effect  of  the 
directed  velocity  of  the  medium  on  the  averaged  macroscopic 
neutron  absorption  cross  section  decreases  (the  ratio  between 
the  cross  sections  at  Vj^=0.0  and  Vj^=4 . 0 is  1.736  at  2000  °K, 
and  1.551  at  7000  °K)  . 

For  sufficiently  high  directed  velocities,  the  decrease 
in  the  averaged  macroscopic  neutron  absorption  cross  section 
with  an  increase  of  the  directed  velocity  of  the  target  atoms 
is  smaller  for  the  case  of  completely  entrained  neutrons  than 
for  the  case  of  uncollided  neutrons  (see  Table  4.8).  At  first 
glance,  this  seems  contradictory  because  with  increasing 
directed  velocity,  for  the  case  of  completely  entrained 
neutrons,  the  target  atom  directed  velocity  yields  both  a 
spectral  hardening  and  a decrease  in  the  effective  microscopic 
cross  section,  both  leading  to  a decrease  in  the  averaged 
macroscopic  cross  section.  On  this  basis,  one  would  expect  the 
decrease  in  the  averaged  macroscopic  cross  section  to  be 
larger  for  the  completely  entrained  case  case  than  for  the 
uncollided  case.  To  explain  the  actual  observed  behavior,  it 
should  first  be  noticed  that  the  decrease  in  the  effective 
microscopic  cross  section  due  to  the  inclusion  of  the  directed 
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velocity  of  the  target  atoms  for  the  completely  entrained  case 
diminishes  with  increasing  energies.  Hence,  the  effective 
microscopic  cross  section  attains  its  value  for  "only  thermal 
effects"  much  more  rapidly  for  the  condition  of  completely 
entrained  neutrons  than  for  the  case  of  uncollided  neutrons 
(see  Chapter  3,  Section  3.5,  Tables  3.3  and  3.4  and  Figures 
3.9  and  3.10).  Consequently,  for  very  high  directed  velocity, 
the  shifting  of  the  low  energy  end  of  the  distribution  curve 
to  higher  energies  significantly  diminishes  the  effect  of 
decreasing  the  averaged  macroscopic  cross  section  by 
decreasing  the  effective  microscopic  cross  section.  Hence,  for 
sufficiently  high  target  atom  directed  velocity,  for  the  case 
of  completely  entrained  neutrons,  the  dominant  effect  is  the 
decrease  in  the  averaged  macroscopic  cross  section  due  to  the 
directed  velocity  induced  spectral  hardening.  This  spectral 
effect  is  smaller  than  the  effect  of  decreasing  the  averaged 
macroscopic  cross  section  by  decreasing  the  effective 
microscopic  cross  section  for  the  case  of  uncollided  neutrons. 

In  the  averaging  process  for  the  completely  entrained 
case,  in  the  low  energy  region,  the  effect  of  the  decreased 
effective  microscopic  cross  section  is  counterbalanced  by  the 
shifting  of  the  neutron  flux  distribution  curve  to  higher 
energies.  In  conclusion  it  can  be  deduced  that  in  the 
averaging  process,  for  sufficiently  high  directed  velocities, 
the  decrease  in  the  effective  microscopic  cross  section  for 
the  uncollided  case  yields  a smaller  averaged  macroscopic 
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cross  section  than  the  combination  of  small  effective  cross 
section  and  directed  velocity  induced  hardening  for  the 
completely  entrained  case. 

As  in  Item  4. 2. 1.2,  the  changes  in  cross  sections  with 
temperature  show  that  as  the  velocity  of  the  medium  is 
increased,  the  temperature  effect  decreases. 

The  results  of  Table  4 . 12  are  represented  graphically  in 
Figure  4.9. 


4. 2. 2. 4 Averaged  macroscopic  neutron  fission  cross  section 

The  results  for  the  averaged  macroscopic  neutron  fission 
cross  section  are  presented  in  Table  4.13. 


Table  4.13  - Averaged  Macroscopic  Neutron  Fission  Cross 

Section  (1/cm)  as  a Function  of  Target  Atom  Temperature 
and  Directed  Velocity  (completely  entrained  neutrons) . 


V. 

2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

0.0 

6.571  E-3 

5.842  E-3 

5.318  E-3 

4.940  E-3 

4.461  E-3 

1.0 

6.123  E-3 

5.503  E-3 

5.043  E-3 

4.706  E-3 

4.282  E-3 

2.0 

5.278  E-3 

4.816  E-3 

4.454  E-3 

4.183  E-3 

3.834  E-3 

3 . 0 

4.465  E-3 

4.117  E-3 

3.832  E-3 

3.612  E-3 

3.334  E-3 

4 . 0 

3.790  E-3 

3.518  E-3 

3.289  E-3 

3.113  E-3 

2.883  E-3 

The  results  in  Table  4.13  present  a behavior  very  similar 
to  the  one  observed  in  Item  4 . 2 . 2 . 3 for  the  averaged 
macroscopic  absorption  cross  section.  The  effect  of  the  target 
atom  directed  velocity  in  reducing  the  averaged  macroscopic 
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Figure  4.9  - Plot  of  Averaged  Macroscopic 

Absorption  Cross  Section  Versus 
Temperature  (completely  entrained 
neutrons) . 
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fission  cross  section  by  reducing  the  effective  microscopic 
fission  cross  section  and  by  spectral  hardening  is  slightly 
smaller  than  in  the  case  of  the  averaged  macroscopic 
absorption  cross  section  because  the  microscopic  fission  cross 
section  drops  off  slightly  faster  with  increasing  energy  than 
the  microscopic  absorption  cross  section.  Clearly  the  spectral 
hardening  with  increasing  temperature  still  contributes  to  the 
decrease  in  the  averaged  macroscopic  fission  cross  section. 


4. 2. 2. 5 Inverse  of  the  average  neutron  velocity 

The  results  for  the  inverse  of  the  average  neutron 
velocity,  for  the  condition  of  completely  entrained  neutrons, 
are  presented  in  Table  4.14. 

Table  4.14  - Inverse  of  the  Average  Neutron  Velocity  as  a 
Function  of  Target  Atom  Temperature  and  Directed 
Velocity  (completely  entrained  neutrons) . 


2000  °K 

3000  °K 

4000  °K 

5000  °K 

7000  °K 

o 

• 

o 

0.3186 

0.2924 

0.2742 

0.2612 

0.2450 

1.0 

0.2989 

0.2772 

0.2617 

0.2505 

0.2364 

2 . 0 

0.2605 

0.2456 

0.2344 

0.2260 

0.2152 

3.0 

0.2231 

0.2133 

0.2055 

0.1995 

0.1916 

4.0 

0.1918 

0.1853 

0.1799 

0.1756 

0.1699 

Table  4.14  shows  that  both  the  increase  in  the  medium's 
temperature  and  the  increase  in  the  directed  velocity  of  the 
target  atoms  present,  as  a final  effect,  a hardened  spectrum. 
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although  this  effect  is  obtained  by  a different  combination  of 
phenomena  as  related  to  the  target  atom  directed  velocity. 

As  expected,  as  the  temperature  increases  the  effect  of 
the  directed  velocity  of  the  target  atoms  diminishes  (the 
ratio  between  the  inverse  velocities  at  V^=0.0  and  V„=4 . 0 is 
1.661  at  2000  °K  and  1.442  at  7000  °K)  and  vice  versa. 

4.3  Conclusions 


Many  conclusions  were  already  drawn  from  the  results  and 
are  included  in  the  analysis  presented  in  Section  4.2  above; 
nevertheless,  some  of  them  that  seem  to  be  most  important  are 
summarized  below: 

1)  The  directed  velocity  of  the  target  atoms  always  leads 
to  a hardening  of  the  spectrum  but  in  a way  completely 
different  from  the  hardening  due  to  an  increase  in  temperature 
(see  Figures  4.4,  4.5,  4.6,  4.7). 

2)  The  directed  velocity  of  the  target  atoms  presents  two 
main  effects:  one  is  on  the  effective  microscopic  cross 
sections  and  the  other  is  on  the  shifting  of  the  neutron 
spectrum.  The  effect  on  the  effective  microscopic  cross 
section  leads  to  different  behaviors  in  the  averaged 
macroscopic  cross  sections  for  the  conditions  of  completely 
entrained  neutrons  and  uncollided  neutrons. 

3)  In  all  the  cases  studied,  the  directed  velocity  of  the 
target  atoms  did  not  remarkably  affect  the  neutron  flux  at  the 
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high  energy  end  of  the  thermal  neutron  distribution,  leading 
only  to  a small  deviation  from  the  neutron  flux  obtained  for 
the  case  which  included  only  temperature  effects. 

4)  The  inclusion  of  the  directed  velocity  of  the  target 
atoms  tends  to  shift  the  low  energy  end  and  the  peak  point  of 
the  distribution  to  a higher  energy.  Both  the  neutron  flux  at 
the  peak  point  of  the  spectrum  and  the  low  energy  end  of  the 
spectrum  itself  increase  progressively  with  increasing  target 
atom  directed  velocity. 


CHAPTER  5 

APPROXIMATE  ANALYTICAL  EXPRESSIONS 


5.1  Introduction 


This  chapter  is  intended  to  present  approximate 
analytical  expressions  that  can  be  used  to  quickly  predict  the 
behavior  of  the  neutron  flux,  for  a system  subjected  to  the 
combined  effects  of  high  target  atom  directed  velocity  and 
temperature,  without  resorting  to  a computer  code.  The  system 
modeled  to  obtain  these  approximate  analytical  approximations 
is  similar  in  geometry  and  composition  to  the  one  presented  in 
the  preceding  chapters. 

The  cases  examined,  as  in  last  chapter,  are  for  the 
conditions  of  completely  entrained  and  uncollided  neutrons. 
However,  for  the  uncollided  condition,  the  expression  for  the 
neutron  flux  distribution  is  already  known;  a standard 
Maxwellian  expression  where  T^  (neutron  temperature)  replaces 
T should  give  very  good  approximation.  Thus,  what  remains  to 
be  done  is  to  develop  approximations  for  the  completely 
entrained  case.  The  developed  expressions  for  the  energy 
dependent  neutron  flux  are  temperature  and  directed  velocity 
dependent. 
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Also  presented  in  this  chapter  are  approximate 
expressions  to  determine  the  effective  microscopic  cross 
sections  in  the  laboratory  system  which  include  the  effects  of 
high  target  atom  directed  velocity  and  temperature.  These 
expressions,  together  with  the  expressions  for  the  energy 
dependent  flux  can  be  used  to  obtain  the  averaged  macroscopic 
cross  sections.  The  resulting  integral,  in  the  averaging 
process  to  find  the  averaged  macroscopic  cross  section,  using 
the  developed  analytical  expressions,  is  very  difficult  to 
solve  analytically,  suggesting  a numerical  integration. 

5.2  Effective  Microscopic  Cross  Section 

Approximate  analytical  expressions  for  the  effective 
microscopic  cross  sections  were  developed  in  Chapter  3 , in  an 
effort  to  obtain  a simplified  procedure  to  treat  the  problem 
of  determining  the  effects  of  high  target  atom  directed 
velocity  and  temperature  on  neutron  interaction  rates,  and  are 
summarized  below. 

5.2.1  Uncollided  Neutrons 

For  the  case  of  uncollided  neutrons,  it  was  found  that 
the  effective  microscopic  cross  section,  which  includes  the 
effects  of  high  target  atom  directed  velocity  and  temperature, 
can  be  obtained  from  the  effective  microscopic  cross  section 
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data  which  includes  only  temperature  effects  by  making  use  of 
two  different  velocity  frames. 

For  a given  neutron  speed,  v,  temperature,  T,  and 
directed  velocity,  the  developed  method  to  obtain  the 
effective  microscopic  cross  section  including  the  directed 
velocity  effects  is  as  follows: 

1)  Use  V and  to  calculate  v^,  the  averaged  neutron 
speed  in  the  moving  frame,  using  Equation  (3-2a)  or  Eguation 
(3-2b) . 

2)  Obtain  the  effective  microscopic  cross  section  in  the 
moving  frame,  considering  only  the  thermal  motion. 
The  raw  data  needed  to  obtain  extracted  from  a 
standard  cross  section  library,  like,  e.g.,  NITAWL  [34]. 

3)  Use  Equation  (3-5)  to  calculate  the  effective 
microscopic  cross  section  in  the  laboratory  system,  cjg^^(v), 
for  a neutron  with  speed  v;  the  resulting  cross  section 
includes  the  effects  of  both  the  target  atom  directed  velocity 
and  the  thermal  velocity  distribution. 

5.2.2  Completely  Entrained  Neutrons 

For  the  case  of  completely  entrained  neutrons,  it  was 
found  that  the  effective  microscopic  cross  section  which 
includes  the  effects  of  high  target  atom  directed  velocity  and 
temperature,  as  in  Item  5.2.1  above,  can  also  be  obtained  from 
the  effective  microscopic  cross  section  data  which  includes 
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only  temperature  effects  by  making  use  of  two  different 
frames . 

For  a given  neutron  speed  in  the  moving  frame 
(temperature  effects  only) , medium  temperature,  T,  and 

directed  velocity,  the  developed  method  to  obtain  the 

effective  microscopic  cross  section  including  the  directed 
velocity  effects  is  the  following: 

1)  For  a combination  of  Vj^  and  V„,  the  averaged  neutron 
speed  in  the  static  frame,  v*,  is  calculated  from  Equation  (3- 
6a)  or  Equation  (3-6b) . 

2)  The  effective  microscopic  cross  section  in  the  moving 

frame,  , for  a temperature  T is  obtained. 

3)  The  effective  microscopic  cross  section  in  the 
laboratory  frame,  a^^^(v*)  , is  calculated  using  Equation  (3-7)  . 

5 . 3 Neutron  Flux 


5.3.1  Overview 

In  Chapter  4 it  was  found  that  the  target  atom  directed 
velocity  induced  a departure  from  the  approximate  Maxwellian 
shape  of  the  neutron  flux  distribution  curve  which  is  obtained 
when  the  system  is  not  subjected  to  target  atom  directed 
velocity  effects. 

The  knowledge  that  a plot  of  In  [0(E)/E]  against  E 
(energy)  for  a Maxwellian  distribution  should  yield  a straight 
line  led  to  the  decision  to  examine  plots  of  In  [0(E) /E] 
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against  E (rather  than  plots  of  0(E)  versus  E)  for  selected 
combinations  of  directed  velocity  and  temperature.  The  purpose 
was  to  observe  the  general  behavior  of  the  curves  of  In 
[0(E)/E],  see  how  they  depart  from  the  Maxwell  ian-shaped 
curve,  and  verify  if  these  curves  were  more  easily  fitted  than 
the  0(E)  versus  E curve  itself.  The  case  of  completely 
entrained  neutrons  is  studied  and  the  results  and  conclusions 
are  presented. 

5.3.2  Completely  Entrained  Neutrons 

Figures  5.1,  5.2,  and  5.3  show  plots  of  In  [0(E)/E] 
versus  E (in  units  of  kT)  for  temperatures  of  2000  °K,  4000  °K, 
and  7 000  °K,  respectively.  The  neutron  flux  at  each  energy 
point  is  obtained  from  BRT-PCM,  and  normalized  according  to 
the  volume  averaged  total  flux,  for  each  case.  In  these  plots 
the  directed  velocities  are  varied  parametrically. 

The  general  tendency  observed  in  these  curves  is  to  show 
a similar  shape,  with  the  low  energy  end  shifted  to  higher 
energies  according  to  the  target  atom  directed  velocity. 

At  low  energies  the  behavior  of  the  curves  can  be 
approximated  by  almost  parallel  vertical  straight  lines,  the 
In  [0(E) /E]  increasing  sharply  with  energy.  The  length  of 
these  almost  vertical  part  of  the  curves  varies  considerably 
according  to  the  temperature  and  directed  velocity.  With 
further  increase  in  energy  the  increasing  behavior  is  reversed 
and  the  In  [0(E) /E]  decreases  continuously  with  increasing 
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Figure  5.1  - Plot  of  Flux  (E)/E  Versus  kT  for 
T = 2000°  (completely  entrained 
neutrons) . 


119 


E (kT) 


Figure  5.2  - Plot  of  glux  (E)/E  Versus  kT  for 
T = 4000°K  (completely  entrained 
neutrons) . 
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Fig  5.3  - Plot  of  Flux  (E)/E  Versus  kT  for 
T = 7000°K  (completely  entrained 
neutrons) . 
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energy,  approaching  a straight  line  parallel  to  the  line  of 
Vn,=0.  The  plots  also  show  that  as  the  directed  velocity 
increases,  for  a given  temperature,  the  In  [0(E)/E]  at  high 
energies  increases  with  increasing  directed  velocity. 

The  behavior  shown  by  these  curves  suggests  that  each  one 
of  them  be  approximately  represented  by  three  straight  lines. 
The  first  represents  the  very  low  energy  region,  in  order  to 
show  the  almost  vertical  line  behavior.  The  third  represents 
the  high  energy  region  where  the  curve  assumes  an  approximate 
straight  line  shape,  parallel  to  the  line  of  V|^,=0.  The  second 
straight  line  represents  the  transition  region  between  the 
first  line  and  third  line. 

The  points  of  intersection  between  the  lines  were 
determined  by  inspecting  the  results  for  each  combination  of 
temperature  and  directed  velocity  and  it  was  found  that  the 
energy  corresponding  to  the  point  where  the  first  straight 
line  intercepts  the  second  straight  line,  E^,  and  the  energy 
corresponding  to  the  point  where  the  second  straight  line 
intercepts  the  third  straight  line,  Ej,  can  be  reasonably 
assumed  to  be  temperature  independent. 

According  to  the  above,  what  is  pursued  is  to  represent 
In  [0(E)/E]  by  an  approximate  analytical  expression  like: 


In 

[0(E) /E]=  A(jE+B(j 

for 

VI 

w 

VI 

0 

w 

and 

(5-la) 

In 

[0(E)/E]=  A^E+B^ 

for 

E^<E<E2, 

and 

(5-lb) 

In 

[0(E)/E]=  A2E+B2 

for 

(5-lc) 
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where  and  Ej  are  as  defined  above  and  Eq,  the  low  energy 
cutoff  can  be  theoretically  calculated. 

The  development  of  this  approach  employed  the  following 
steps : 

1)  From  the  results  obtained  and  from  the  plots  it  is 
found  that  E^  can  be  reasonably  assumed  to  be  temperature 
independent  and  E^  was  determined  as  a function  of  V|^  by  making 
a polynomial  fitting  of  the  values  obtained  for  E^  according 
to  for  each  case. 

2)  From  the  results  obtained  and  from  the  plots  it  is 
found  that  Ej  is  temperature  independent  and,  as  for  E^,  Ej  was 
determined  as  a function  of  V|^  by  making  a polynomial  fitting 
of  the  values  obtained  for  according  to  for  each  case. 

3)  The  In  [0(E)/E]  at  E^ , B(T,V„)  , was  determined  as: 

B(T,VJ=B(T)-B(V„)  , where  (5-2) 

B(T) : function  determined  according  to  the  values 

obtained  for  B(T,V|^)  for  V^=0  (only  temperature  dependent)  , and 
B(Vn,)  : function  determined  by  correlating  the  difference 
between  B(T)  and  B(T,V|^)  , for  all  the  studied  combinations  of 
T and  and  by  suitable  fitting  of  these  differences. 

4)  The  In  [0(E)/E]  at  E2,  B'(T,V,^),  was  determined  as: 

B' (T,V„)=B*  (T)-B' (V„)  , where  (5-3) 

B'  (T)  : function  determined  by  the  fitting  of  B'  (T,V,^)  for 

V^=0  (only  temperature  dependent) , and 
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B' (V^)  : function  determined  by  fitting  the  difference 

between  B'(T)  and  B'(T,Vj^),  for  all  the  studied  combinations 
of  T and  Vj^. 

5)  The  slope  of  the  first  approximate  straight  line, 
Aq(T,V„)  , was  determined  as: 

' where  (5-4) 

Aq(T)  : function  determined  by  the  fitting  of  Ag(T,Vj,)  for 
Vj^=0  (only  temperature  dependent)  , and 

fg(AQ(V|^)  ,V„)  : function  determined  according  to  the  values 
obtained  for  Aq(V^),  for  each  case.  Ag(V^)  is  the  result 
obtained  from  the  difference  between  Ag(T)  and  Ag(T,V„)  , 
multiplied  by  a function  of  for  all  the  studied 

combinations  of  T and  V„.  fg  (Ag  (V„)  , V„)  is  expressed  by: 


fo(Ao(V„),V„)  = (Ag(V„)*(V„)i/2  (5_4a) 

6)  The  second  straight  line  was  defined  by  its  points 
B(T,V„)  and  B'(T,V„),  yielding: 

Ai(T,V„)  = [B'  (T,V„)-B(T,V„)  ]/(E2-E^)  (5-5a) 

Bi(T,V„)=[E2*B(T,V„)-Ei*B'  (T,V„)]/(E2-E^)  (5-5b) 

7)  The  slope  of  the  third  approximate  straight  line, 
A2(T,V^)  , was  determined  as: 

A2(T,VM)=A2(T)+f2(A2(V„)  ,V„)  where  (5-6) 


A2(T) : function  determined  by  the  fitting  of  A2(T,V„)  for 
Vn,=0  (only  temperature  dependent)  , and 

f2(A2(V^)  ,V^)  : function  determined  according  to  the  values 
obtained  for  A2(V^),  for  each  case.  A2(V^)  is  the  result 
obtained  from  the  difference  between  A2(T)  and  A2(T,V|^), 
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multiplied  by  a function  of  V^,  for  all  the  studied 
combinations  of  T and  V^.  f2  (A2  (V^)  , V„)  is  expressed  by; 

f2(A2(V„)  ,v„)  = (A2(V„)*(V„)  (5-6a) 

It  should  be  noticed  that  the  approximate  expressions  are 
developed  in  order  to  treat  the  directed  velocity  effects  like 
a perturbation  on  the  results  for  temperature  effects  only. 
This  approach  permits  the  quick  determination  of  whether  there 
is  a need  to  include  the  directed  velocity  effects. 

The  results  of  the  above  fittings  and  correlations  are  as 
follows : 

(V„)  =3 . 56229E-2-6 . 42956E-3*V„+2 . 4387254E-2*V„2+ 
1.057358E-3*V„^-1.5185E-4*V„'^ 

E2(Vm)=0. 106915-8. 97  668  3E-2*V„+1. 50017333E-1*V„2- 
4 . 845817E-2*V„^+5 . 739667E-3*V„^ 

B(T) =9 . 885733-5. 3 13E-1*T' +3 . 4966667E-2*T ' ^ 

B(V„)  =-0. 006829  + 0. 3 0569 05*V„+0. 23 17857*V„^-0. 04 0333 
B' (T) =9. 349333-0. 377 *T' +2 . 3 666667 *T'^ 

B'  (V„)  =-0.1*V„+0. 395833 *V„^-0. 13 *V„^+0.01416667*V„^ 

Aq(T)=4 .51533+3 .431*T'-1. 9933 3E-1*T' 2 
Aq  (V„)  =14 . 47+55 . 455*V„-2  0 . 68*V„^+2 . 255*V„^ 

Ao(T,V„)=Ao(T)+Aq(V„)*V;/2 

A2(T) =-7 .2777333+1. 1718 *T' -8. 04 6667E-2*T'^ 

A2(V„)=0. 1264-0. 04 51*V„+0. 02 64*V„2-0. 003 1*V„^ 
A2(T,V„)=A2(T)+A2(V„)*V„. 


In  the  above  equations: 
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- V|^,:  units  of  2200  m/s. 

- : units  of  eV. 

- E2:  units  of  eV. 

- T'=T/1000. 

- T:  units  of  °K. 

- Aq(T,V„),  Ag(T),  A,(T,V„),  A2(T,V„),  and  A2(T)  : units  of 

1/eV 


The  three  straight  lines  which  represent  the  In  [(p(E)/E] 
can  be  rewritten  as  follows: 


In  [0(E)/E]=  Aq(T,VJ  (E-E^)+B(T,V„)  , (5-7a) 

for  E<E^ , and 

In  [0(E)/E]=  { [B*  (T,V„)-B(T,V„)  ]/(E2-E,)  )E+ 

[E2*B(T,V„)-E^*B*  (T,V„)  ]/(E2-E^)  , (5-7b) 

for  E^<E<E2,  and 

In  [0(E)/E]=  A2(T,V„)  (E-E2)+B'  (T,V„)  , (5-7c) 

for  E>E2- 

According  to  the  above,  the  neutron  flux  at  the  energy  E, 
0(E) , including  the  effects  of  high  target  atom  velocity  and 
temperature  can  be  determined,  approximately,  by: 

0(E)=  E*exp[Ao(T,V„)  (E-E^)+B(T,V„)  ] , (5-8a) 

for  Eq<E<E^,  and 

0(E)=  E*exp{{[B'  (T,V„)-B(T,V„)  ]/(E2-E^)  )E+ 
[E2*B(T,V„)-E^*B'  (T,V„)  ]/(E2-E^)  }, 


(5-8b) 


126 


for  E^<E<E2,  and 

0(E)=  E*exp[A2(T,V„)  (E-E2)+B' (T,V„)  ],  (5-8c) 
for  E>E2- 

These  approximate  expressions  for  the  In  [0(E) /E]  were 
tested  against  the  calculated  results  obtained  from  BRT-PCM 
for  T=2000  °K  and  V„=l,  T=4000  °K  and  V„=2 , and  T=7000  °K  and 
Vn,=4 . The  results  obtained  from  the  tests  are  presented  in 
Tables  5.1,  5.2  and  5.3,  respectively. 

The  results  from  these  tables  show  a general  tendency  for 
the  approximate  expressions  to  underpredict  the  In  [0(E) /E] 
values . 

The  error  involved  is  generally  reasonable  at  the  low 
energy  end  but  can  increase  noticeably  at  the  second  energy 
point,  according  to  the  combination  of  directed  velocity  and 
temperature.  This  behavior  was  expected  because  of  the  low 
order  straight  line  approximation  used  in  the  fitting  at  the 
low  energy  end,  where  the  straight  line  was  interpolated 
between  the  real  values.  The  presented  errors  in  the  low 
energy  end  increase  with  increasing  temperature  and  directed 
velocity.  Nevertheless,  they  can  be  considered  acceptable  if 
it  is  noticed  that  the  neutron  flux  has  very  small  values  at 
these  energies,  and  tends  to  decrease  with  increasing  directed 
velocity. 

The  error  involved  is  small  at  intermediate  energies  and 
at  the  high  energy  end  and  tends  to  increase  slightly  with 
increasing  temperature  and  target  atom  directed  velocity. 
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Table  5.1  - Results  from  the  Test  of  Approximate 

Expressions  for  Tggg=2  000  °K  and  V„=l  (completely- 
entrained  neutrons) . 


ENERGY 

(eV) 

Ln  [0(E)/E] 
(calculated) 

Ln  [0(E)/E] 
(approximate) 

ERROR 

(%) 

2.703  E-2 

6.690 

6.760 

1.04 

3.574  E-2 

8 . 000 

7.310 

-8.62 

5.448  E-2 

8.500 

8.473 

-0.32 

8.442  E-2 

8.570 

8.489 

-0.94 

1.244  E-1 

8 . 510 

8.510 

0.00 

2.336  E-1 

8.040 

7.948 

-1.14 

3.812  E-1 

7.300 

7.188 

-1.53 

6.745  E-1 

5.760 

5.677 

-1.44 

9.177  E-1 

4.470 

4.424 

-1.03 

1.199  E+0 

3 . 030 

2.974 

-1.85 
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Table  5.2  - Results  from  the  Test  of  Approximate 

Expressions  for  Tggg=4000  °K  and  V„=2  (completely 
entrained  neutrons) . 


ENERGY 

(eV) 

Ln  [0(E)/E] 
(calculated) 

Ln  [0(E)/E] 
(approximate) 

ERROR 

(%) 

1.029  E-1 

4.510 

4.740 

5.10 

1.109  E-1 

6.210 

5.560 

-10.46 

1.263  E-1 

7.070 

7.110 

0.56 

1.504  E-1 

7.500 

7.230 

-3 . 60 

1.850  E-1 

7.680 

7.410 

-3.50 

2.316  E-1 

7 . 660 

7.653 

-0.09 

2.892  E-1 

7.550 

7.443 

1 

H 

• 

to 

3.569  E-1 

7.350 

7.200 

-2.04 

5.224  E-1 

6.780 

6.593 

-2.75 

8.433  E-1 

5.560 

5.423 

-2.46 

1.251  E+0 

4 . 020 

3 . 937 

-2.06 

1.570  E+0 

2.890 

2.772 

-4.40 
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Table  5.3  - Results  from  the  Test  of  Approximate 

Expressions  for  Tggg=7  000  °K  and  V|,^=4  (completely 
entrained  neutrons) . 


ENERGY 

(eV) 

Ln  [0(E)/E] 
(calculated) 

Ln  [0(E)/E] 
(approximate) 

ERROR 

(%) 

4.065  E-1 

2.610 

2.880 

10.34 

4.144  E-1 

4.500 

3.820 

-15.11 

4.289  E-1 

5.500 

5.530 

-0.54 

4.503  E-1 

6.100 

5.798 

-4.95 

4.792  E-1 

6.500 

6.161 

-5.21 

5.162  E-1 

6.660 

6.630 

-0.45 

5.619  E-1 

6.710 

6.523 

-2.78 

6.175  E-1 

6.670 

6.393 

-4.15 

7.625  E-1 

6.370 

6.054 

-4.96 

1.068  E+0 

5.530 

5.340 

-3.43 

1.324  E+0 

4.830 

4.741 

-1.84 

1.784  E+0 

3.730 

3.666 

-1.71 
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Nevertheless,  the  agreement  between  the  calculated  results  and 
the  approximate  results  can  be  considered  reasonable. 

5.3.3  Range  of  Effectiveness  for  the  Approximate  Expressions 

The  developed  approximate  expressions  are  expected  to 
give  reasonable  results  over  the  range  of  temperatures  and 
directed  velocities  studied,  i.e.,  from  2000  °K  to  7000  °K  and 
for  from  0.0  to  4.0.  This  range  can  be  expanded  by 
generating  a new  data  set,  including  the  desired  temperatures 
and/or  directed  velocities,  and  following  the  approach 
described  above  to  develop  new  expressions. 

It  should  be  noticed  that  it  might  be  possible  to  improve 
the  above  expressions  for  the  neutron  flux  by  using  a larger 
set  of  data  (speed)  points  for  each  case  and  by  using 
different  (higher  order)  or  more  functions  in  the  fitting 
process. 


5.4  Averaged  Macroscopic  Cross  Section 

The  developed  analytical  expressions  can  be  used  to 
determine  approximate  values  for  the  averaged  macroscopic 
cross  sections.  As  pointed  out  before,  the  analytical 
integration  required  to  perform  the  averaging  process  is  very 
cumbersome,  but  the  averaged  macroscopic  cross  sections  can  be 
obtained  by  calculating  the  effective  microscopic  cross 
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section  and  the  neutron  flux  at  many  energy  points  and  using 
these  values  to  perform  a numerical  integration. 

It  should  be  noticed  that  the  flux  obtained  from  the 
above  analytical  expressions  which  is  to  be  used  as  a 
weighting  factor  in  the  averaging  process  to  determine  the 
averaged  macroscopic  cross  section  is  only  energy-dependent; 
it  doesn't  include  spatial  (e.g.,  self-shielding)  effects. 
This  limitation  isn't  significant  for  the  UTVR/MHD  Disk 
Generator  or  for  most  gas  core  reactors  which  tend  to  present 
a reasonably  flat  flux  distribution;  however,  this  can  be 
important  for  certain  gas  core  reactors  (e.g.,  some  gas  core 
reactor  propulsion  concepts)  which  exhibit  a spatially 
dependent  flux. 


CHAPTER  6 

CONCLUSIONS  AND  RECOMMENDATIONS 
6 . 1 Summary 

The  basic  intent  of  this  dissertation  was  to  determine 
the  effects  of  high  target  atom  directed  velocity  and 
temperature  on  interaction  rates  for  a flowing  fissioning  gas. 

The  system  used  as  a model  in  this  study  is  a gas  cavity 
region  where  a low  density  mixture  of  fissioning  fuel  and 
working  fluid  is  in  thermodynamic  equilibrium  at  a high 
temperature  and  flowing  at  constant  velocity.  The  gas  cavity 
region  is  surrounded  by  a moderator  which  is  the  thermal 
neutron  source. 

A brief  review  of  gas  core  reactors  was  presented  first 
and  since  the  main  focus  of  this  work  is  the  thermal  energy 
region,  a concise  historical  review  of  neutron  thermal iz at ion 
was  also  presented. 

The  study  continued  with  the  development  of  the  exact 
analytical  equations  to  calculate  the  effective  microscopic 
cross  section  and  scattering  kernel,  including  the  effects  of 
high  target  atom  directed  velocity  and  temperature.  The 
equations  obtained  in  this  approach,  except  under  special 
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conditions,  are  impossible  to  solve  analytically  and  even 
their  solution  by  numerical  methods  is  very  difficult, 
suggesting  possibly  a solution  employing  Monte  Carlo  methods. 
Hence,  an  approximate  approach,  was  developed  to  solve  the 
problem,  based  on  the  exact  analytical  equations  and  using  an 
existing  computer  code  properly  modified  to  include  the 
desired  conditions  of  high  target  atom  directed  velocity  and 
temperature.  This  approximate  approach  includes  the  use  of  two 
different  reference  frames,  one  is  the  static  or  laboratory 
frame  and  the  other  is  a moving  frame  that  has  the  same 
directed  velocity  as  the  target  atoms.  A product  of  this 
development  was  the  computer  code  BRT-PCM,  which  runs  on  a 
personal  computer  and  solves  the  integral  form  of  the  neutron 
transport  equation  in  the  thermal  energy  region  to  include  the 
above  mentioned  desired  conditions.  The  developed  approximate 
approach  is  appropriate  for  the  analysis  of  two  extreme 
conditions  that  were  studied  in  detail;  1)  the  condition  of 
uncollided  neutrons  and,  2)  the  condition  of  completely 
entrained  neutrons. 

Using  the  developed  approach,  the  effects  of  the  target 
atom  directed  velocity  and  temperature,  on  the  averaged 
macroscopic  cross  sections,  for  the  uncollided  and  completely 
entrained  neutron  conditions  were  determined.  For  the 
completely  entrained  neutron  condition  the  inverse  of  the 
average  neutron  velocity  and  thermal  neutron  flux  were  also 


determined. 
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From  these  results,  a set  of  approximate  analytical 
expressions  was  derived  to  determine  the  thermal  neutron  flux 
distribution  and  the  effective  microscopic  cross  sections 
which  include  the  effects  of  high  target  atom  directed 
velocity  and  temperature.  Using  these  analytical 
approximations  the  averaged  macroscopic  cross  sections  can 
also  be  determined. 


6 . 2 Results 


The  new  results  brought  out  in  this  dissertation  can  be 
summarized,  for  each  of  the  two  major  conditions  described 
above,  as  follows: 

6.2.1  Uncollided  Neutrons 

6. 2. 1.1  Averaged  macroscopic  neutron  cross  sections 

It  was  first  shown  that  the  directed  velocity  of  the 
target  atoms  has  the  effect  of  increasing  the  hydrogen 
effective  thermal  microscopic  scattering  cross  section  for  the 
case  of  uncollided  neutrons.  According  to  the  developed 
approach,  the  effective  microscopic  cross  section  including 
the  target  atom  directed  velocity  and  temperature  effects  is 
obtained  by  scaling  the  corresponding  effective  microscopic 
cross  section  which  includes  only  thermal  effects.  Hence,  the 
results  for  the  averaged  thejrmal  macroscopic  cross  section. 
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when  the  target  atom  directed  velocity  is  included,  represent 
a combination  of  the  target  atom  directed  velocity  effects  on 
each  of  the  effective  microscopic  cross  sections  for  all  the 
isotopes  included  in  the  fuel  gas.  As  a result,  it  was  shown 
that  the  target  atom  directed  velocity  effect  is  to  increase 
the  averaged  macroscopic  scattering  cross  section  (since 
hydrogen  is  the  dominant  isotope  as  far  as  scattering  is 
concerned)  and  to  decrease  the  averaged  macroscopic  absorption 
and  fission  cross  sections.  It  should  be  noted  that  for 
sufficiently  high  temperatures,  the  hydrogen  effective 
microscopic  scattering  cross  section  decreases  with  increasing 
energy.  It  was  also  shown  that  for  materials  which  present  a 
completely  1/v  microscopic  cross  section  behavior,  the 
effective  microscopic  cross  section  is  not  affected  by  the 
inclusion  of  the  target  atom  directed  velocity. 

It  was  noted  that  due  to  the  different  energy  dependent 
behavior  presented  by  the  microscopic  absorption,  scattering 
and  fission  cross  sections  (i.e.,  the  way  they  behave  with 
energy,  e.g.,  the  microscopic  absorption  cross  section  falling 
off  faster  with  increasing  energy  than  the  microscopic 
scattering  cross  section  and  the  microscopic  fission  cross 
section  falling  off  faster  with  increasing  energy  than  the 
microscopic  absorption  cross  section) , the  observed  directed 
velocity  effects  (decrease  or  increase)  on  the  averaged 
macroscopic  cross  sections  are  different  for  each  type  of 


reaction. 


136 


It  was  also  noted  that:  1)  for  a specific  temperature,  as 
the  directed  velocity  increases,  the  target  atom  directed 
velocity  effects  on  the  averaged  macroscopic  cross  sections 
also  increase,  as  expected,  and  2)  for  a specific  directed 
velocity,  as  the  temperature  increases,  the  thermal  induced 
hardening  of  the  flux  decreases  the  flux  at  low  energies  and 
increases  the  flux  at  high  energies,  causing  the  effect  of  the 
directed  velocity  on  the  averaged  macroscopic  cross  section  to 
be  less  pronounced  because  the  neutron  flux  is  decreased  where 
the  target  atom  directed  velocity  effect  on  the  effective 
microscopic  cross  section  is  more  notable. 

For  the  uncollided  neutron  condition,  at  2000  °K,  the 
results  for  V^=2.0  present  a maximum  difference  of  11%  from 
the  results  for  V^=0.0;  this  difference  decreases  to  9.3%  at 
7000  °K.  The  results  for  V^=4 . 0 present  a difference  as  high 
as  39.5%  (for  the  averaged  macroscopic  scattering  cross 
section)  from  the  results  for  V|^=0.0,  at  2000  °K.  This 
difference  decreases  to  22.7%  at  7000  °K,  showing  that,  unlike 
for  Vj^=2.0,  an  increase  of  5000  °K  in  temperature  leads  to  a 
significant  decrease  in  the  directed  velocity  effects.  The 
percent  difference  between  the  results  for  V„=4 . 0 and  the 
results  for  V„=0.0  at  2000  °K  are  more  significant  for  the 
averaged  macroscopic  scattering  (39.5%)  and  absorption  (32.6%) 
cross  sections  than  for  the  averaged  macroscopic  fission  cross 
section  (31.1%) . At  7000  °K  the  percent  difference  between  the 
results  for  Vj^=4 . 0 and  the  results  for  V„=0.0  are  more 
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significant  for  the  averaged  macroscopic  fission  (22.7%)  and 
absorption  (14.2%)  cross  sections  than  for  the  averaged 
macroscopic  scattering  cross  section  (9.9%).  This  is  because 
of  the  above  explained  different  effects  of  the  directed 
velocity  on  the  effective  microscopic  cross  sections  which  is 
in  turn  due  to  the  different  behavior  of  each  cross  section 
with  energy. 

6. 2. 1.2  Approximate  analytical  expressions 

The  developed  approximate  analytical  approach  to 
calculate  the  effective  microscopic  cross  section  for  the  case 
of  uncollided  neutrons  presented  results  that  are  in  very  good 
agreement  with  the  only  known  previously  available  and  very 
limited  results.  These  expressions  were  tested  for  limiting 
conditions  of  temperature  and  directed  velocity  and  the 
behavior  predicted  using  the  approach  developed  for  uncollided 
neutrons  is  the  same  behavior  predicted  by  the  exact 
expressions,  under  the  same  conditions. 

6.2.2  Completely  Entrained  Neutrons 
6. 2. 2.1  Neutron  flux 

For  the  completely  entrained  neutron  condition,  it  was 
observed  that  the  addition  of  the  directed  velocity  of  the 
target  atoms  shifts  the  theirmal  neutron  flux  distribution  to 
higher  energies.  As  the  directed  velocity  of  the  target  atoms 
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is  increased  this  effect  becomes  more  pronounced.  It  was  also 
observed  that:  1)  the  inclusion  of  the  directed  velocity 
results  in  a shifting  of  the  low  energy  end  of  the  curve  to  an 
energy  that  is  much  higher  than  the  low  energy  end  for  the 
corresponding  case  which  includes  only  temperature  effects;  2) 
for  a specific  temperature,  the  neutron  flux  distribution 
exhibits  a peak  at  an  energy  that  is  higher  than  for  the  case 
of  temperature  effect  only;  this  flux  peaking  energy  increases 
with  increasing  target  atom  directed  velocity;  3)  at  low 
target  atom  directed  velocity,  the  neutron  flux  distribution 
exhibits  a peaking  which  is  slightly  higher  than  for  the  case 
of  temperature  only;  with  increasing  directed  velocity  the 
increase  in  the  peak  becomes  more  significant;  4)  as  the 
temperature  increases,  the  effect  of  the  directed  velocity  of 
the  target  atoms  is  less  pronounced  (i.e.,  it  is  somewhat 
masked  by  the  increase  in  temperature)  and  the  ratio  of  the 
neutron  flux  with  directed  velocity  effects  to  the  neutron 
flux  without  directed  velocity  effects  (temperature  effect 
only) , at  the  energy  of  the  distribution  curve  peak, 
decreases.  This  type  of  behavior  is  similar  to  that  in  most 
other  cases  where  the  effects  of  increased  temperature  versus 
the  effect  of  increased  directed  velocity  of  the  target  atoms 
were  examined;  5)  for  sufficiently  high  neutron  energies  the 
effect  of  the  directed  velocity  of  the  target  atoms  is  small 
and  the  distribution  curve  exhibits  essentially  the  shape  it 
would  have  for  temperature  effects  only,  with  only  a small 
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deviation  in  magnitude  from  the  standard  distribution  and 
the  high  energy  end  of  the  distribution  curve  is  slightly 
shifted  into  a higher  energy;  and  6)  for  a specific  directed 
velocity  of  the  target  atoms,  the  energy  corresponding  to  the 
low  energy  end  of  the  distribution  curves  is  always  the  same, 
independent  of  the  medium's  temperature. 

In  conclusion,  the  behavior  observed  for  the  neutron  flux 
with  the  inclusion  of  the  directed  velocity  of  the  target 
atoms  makes  sense  physically  if  it  is  realized  that  neutrons 
in  the  static  frame,  initially  with  energies  (velocities) 
smaller  than  the  directed  velocity  of  the  target  atoms,  will 
have  to  suffer  upscatterings  to  become  completely  entrained 
and  these  upscatterings  will  tend  to  accumulate  neutrons 
around  an  energy  peak  according  to  the  temperature  of  the 
medium  and  the  directed  velocity  of  the  target  atoms.  As  the 
energy  of  the  neutrons  increases,  the  effect  of  the  directed 
velocity  on  the  effective  microscopic  cross  sections  and  on 
the  neutron  number  density  distribution  is  diminished;  the 
upscattering  is  less  pronounced  than  for  small  neutron 
energies  and  this  explains  the  behavior  of  the  neutron  flux  at 
high  energies. 

6. 2. 2. 2 Averaged  macroscopic  neutron  cross  sections 

With  regard  to  the  behavior  of  the  averaged  macroscopic 
cross  sections  for  the  case  of  completely  entrained  neutrons, 
it  was  observed  that:  1)  for  a fixed  temperature,  an  increase 
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in  the  directed  velocity  of  the  target  atoms  always  yields  a 
decrease  in  the  averaged  macroscopic  neutron  cross  sections; 
2)  in  contrast  to  the  uncollided  case,  the  target  atom 
directed  velocity  effect  is  to  decrease  all  the  effective 
microscopic  (macroscopic)  cross  sections,  particularly  in  the 
low  energy  region;  3)  the  directed  velocity  hardens  the 
neutron  spectrum  by  shifting  the  low  energy  end  of  the 
distribution  curve  to  a higher  energy  (which  contributes  to  a 
further  decrease  in  the  averaged  macroscopic  cross  section) ; 
4)  as  the  directed  velocity  of  the  target  atoms  increases,  the 
flux  tends  to  decrease  slightly  in  the  vicinity  of  the  low 
energy  end  of  the  distribution  curve  which  diminishes  the 
directed  velocity  effect  of  decreasing  the  effective 
microscopic  cross  section  at  low  energies;  5)  the  shifting  of 
the  flux  distribution  peaking  to  higher  energies  with 
increasing  directed  velocity  contributes  to  a further  decrease 
in  the  averaged  macroscopic  cross  section;  6)  for  a fixed 
directed  velocity,  as  the  temperature  increases,  the  flux 
tends  to  shift  towards  higher  energies  and  the  weighted 
contribution  of  the  reduced  low  energy  effective  microscopic 
cross  sections  is  smaller  than  before  (at  lower  temperature) . 
Depending  on  the  behavior  of  the  microscopic  scattering  cross 
sections  with  energy,  this  effect  may  lead  to  an  increase  in 
the  averaged  macroscopic  scattering  cross  section,  until  a 
temperature  is  reached  where  the  averaged  macroscopic 
scattering  cross  section  attains  a maximum  value;  for  higher 
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temperatures  the  spectral  hardening  due  to  the  temperature 
effect  causes  the  averaged  macroscopic  scattering  cross 
section  to  decrease  again.  With  further  increases  in 
temperature,  the  neutron  spectral  hardening  continues  to  cause 
the  averaged  macroscopic  scattering  cross  section  to  decrease; 
and  7)  for  a specific  directed  velocity,  an  increase  in 
temperature  always  leads  to  a decrease  in  the  averaged 
macroscopic  absorption  and  fission  cross  sections. 

For  the  completely  entrained  neutron  condition,  the 
results  for  V^=2 . 0 present  a maximum  difference  of  19.7%  from 
the  results  for  V^=0.0,  at  2000  °K  (for  the  averaged 
macroscopic  fission  cross  section) ; this  difference  decreases 
to  14%  at  7000  °K.  The  results  for  Vj^=4 . 0 present  a difference 
as  high  as  42.4%  (for  the  averaged  macroscopic  fission  and 
absorption  cross  sections)  from  the  results  for  V=0.0,  at 
2000  °K.  This  difference  decreases  to  35.5%  at  7000  °K,  showing 
that,  as  for  the  case  of  uncollided  neutrons,  an  increase  of 
5000  °K  in  temperature  leads  to  a significant  decrease  in  the 
directed  velocity  effects.  At  low  temperatures,  the  percent 
differences  in  the  results  for  V„=4 . 0 from  the  results  for 
Vm=0.0  are  slightly  smaller  for  the  averaged  macroscopic 
scattering  cross  section  than  for  the  averaged  macroscopic 
fission  and  absorption  cross  sections;  with  increasing 
temperature  these  differences  decrease.  The  percent 
differences  for  the  averaged  macroscopic  fission  and 
absorption  cross  sections  at  some  given  relative  to  Vj,=0  are 
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approximately  the  same  for  all  the  combinations  of  temperature 
and  directed  velocity.  This  is  due  to  the  above  explained 
directed  velocity  effects  on  the  neutron  flux  distribution  and 
on  the  effective  microscopic  cross  sections  and  the  behavior 
of  each  effective  microscopic  cross  section  with  energy. 

6. 2. 2. 3 Approximate  analytical  expressions 

The  results  obtained  from  the  approximate  analytical 
expressions  for  the  effective  microscopic  cross  section  for 
the  case  of  completely  entrained  neutrons  could  not  be 
compared  with  any  other  previous  results  because  there  are  no 
data  available.  As  for  the  case  of  uncollided  neutrons,  these 
expressions  were  tested  for  limiting  conditions  of  temperature 
and  directed  velocity  and  the  behavior  predicted  using  the 
expressions  obtained  from  the  developed  approach  is  the  same 
behavior  predicted  by  the  exact  expressions,  under  the  same 
conditions . 

The  approximate  analytical  expression  developed  for  the 
thermal  neutron  flux  yield  acceptable  accuracy  over  the  range 
of  temperatures  and  directed  velocities  analyzed. 

The  flux  obtained  from  the  approximate  analytical 
expression  to  be  used  as  a weighting  factor  in  the  averaging 
process  to  determine  the  averaged  macroscopic  cross  section  is 
only  energy  dependent;  it  does  not  include  spatial  effects. 
This  limitation  isn't  significant  for  the  UTVR/MHD  Disk 
Generator  or  for  gas  core  reactors  which  present  a flux 
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distribution  which  can  be  considered  reasonably  "flat”,  but  it 
can  be  important  for  some  gas  core  reactors  systems  (e.g., 
certain  propulsion  concepts)  which  present  a flux  with 
significant  spatial  dependence. 

The  approximate  analytical  expressions  for  the  thermal 
neutron  flux  also  permit  one  to  quickly  determine  whether 
there  is  a need  to  include  the  directed  velocity  effects  in 
the  thermal  neutron  flux  and  cross  section  calculations. 

6.3  Conclusions 

This  work  brought  out  many  conclusions  about  the  high 
target  atom  directed  velocity  and  temperature  effects  on  the 
neutron  flux,  effective  microscopic  cross  sections,  and 
averaged  macroscopic  cross  sections.  High  target  atom 
temperature  effects  were  already  discussed  in  other  works. 
This  work  sought  to  analyze  the  combination  of  high  directed 
velocity  and  temperature  effects  and,  whenever  possible,  to 
and  compare  them.  Hence,  the  most  significant 
conclusions  are  those  related  to  the  high  directed  velocity 
and  some  of  these  are  summarized  below: 

a)  For  the  completely  entrained  neutron  condition,  the 
directed  velocity  effect  is  to  decrease  the  effective 
microscopic  cross  section.  This  decrease  diminishes  as  the 
energy  increases.  For  the  uncollided  neutron  condition,  the 
directed  velocity  might  induce  an  increase  or  decrease  in  the 
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effective  microscopic  cross  section,  depending  on  the  energy 
dependent  behavior  of  the  effective  microscopic  cross  section. 
This  effect  decreases  as  the  energy  increases. 

b)  For  the  completely  entrained  case,  the  inclusion  of 
the  target  atom  directed  velocity  always  leads  to  a shifting 
of  the  low  energy  end  of  the  distribution  curve  to  higher 
energies,  resulting  in  a hardened  distribution;  the  effects  in 
the  high  energy  side  of  the  distribution  are  small.  The 
directed  velocity  induced  hardening  is  completely  different 
from  the  temperature  induced  hardening. 

c)  For  the  conditions  studied,  the  directed  velocity 
effect  on  the  effective  microscopic  scattering  cross  section 
is  reversed  as  the  neutrons  pass  from  the  uncollided  condition 
to  the  completely  entrained  condition.  For  the  effective 
microscopic  absorption  and  fission  cross  sections  this  effect 
depends  on  the  energy  dependent  behavior  of  these  microscopic 
cross  sections  and  on  the  range  of  energy  studied. 

d)  The  neutron  flux  at  the  low  energy  end  of  the  spectrum 
decreases  and  the  low  energy  end  of  the  spectrum  itself 
increases  progressively  as  the  neutrons  pass  from  the 
uncollided  condition  to  the  completely  entrained  condition. 

e)  The  directed  velocity  affects  the  averaged  macroscopic 
cross  section  by  affecting  the  effective  microscopic  cross 
section  and,  for  the  case  of  completely  entrained  neutrons,  by 
also  affecting  the  spectrum  distribution.  The  averaged 
macroscopic  cross  sections  are  usually  higher  for  the 
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uncollided  condition  than  for  the  completely  entrained 
condition,  showing  that  there  is  a tendency  to  decrease  the 
averaged  macroscopic  cross  sections  as  the  neutrons  pass  from 
the  uncollided  to  the  completely  entrained  condition.  This 
implies  a decrease  of  reaction  rates  as  the  neutrons  pass  from 
the  uncollided  to  the  completely  entrained  condition. 

f)  The  directed  velocity  can  sensibly  affect  the  neutron 
i^t®^3.ction  rates,  usually  leading  to  a decrease.  For  a given 
temperature,  the  decrease  in  the  interaction  rates  tends  to 
increase  with  increasing  directed  velocity. 

g)  For  some  systems,  depending  on  the  configuration,  the 
inclusion  (or  increase)  of  the  directed  velocity  may  lead  to 
a condition  where,  due  to  the  increase  of  the  mean  free  path, 
the  neutron  escapes  from  the  gaseous  fuel,  after  a few 
<^ollisions,  without  attaining  the  condition  of  complete 
entrainment.  This  premature  escape  can  also  diminish  the 
number  of  fissions  in  the  gaseous  fuel. 

h)  The  obtained  results  are  for  systems  where  the  gaseous 
fuel  and  working  fluid  are  in  thermodynamical  equilibrium, 
flowing  at  high  directed  velocity,  and  the  gas  mixture  has 
very  low  density.  These  results  also  apply  to  any  system  which 
presents  a combination  of  high  directed  velocity,  temperature, 
and  density,  undergoing  a highly  turbulent  flow,  that  makes 
the  assumption  of  a flat  velocity  profile  reasonable. 

i)  The  results  obtained  for  the  completely  entrained  case 
and  for  the  materials  involved  in  the  modeled  system,  are  for 
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a neutron  source  at  2000  °K.  An  increase  in  the  neutron  source 
temperature  leads  to  a shifting  of  the  low  energy  end  of  the 
distribution  curve  to  higher  energies  which  diminishes  the 
directed  velocity  effects. 

j ) Although  low  energy  resonances  effects  were  not 
considered  in  this  work,  fuel  temperature  (Doppler)  broadening 
combined  with  the  directed  velocity  effect  might  lead  to  an 
increase  in  the  effective  microscopic  cross  section  at  low 
energies  resulting  in  a increase  of  the  averaged  cross 
sections  due  to  an  increase  in  the  resonance  effective 
integral  and,  consequently  in  an  enhancement  of  the 
interaction  rates. 

k)  For  each  temperature  there  is  a limiting  directed 
velocity  below  which  the  directed  velocity  effects  can  be 
reasonably  neglected.  As  the  temperature  increases,  this 
limiting  directed  velocity  also  increases.  The  same  reasoning 
applies  for  a given  directed  velocity  and  corresponding 
limiting  temperatures. 

l)  In  the  range  of  temperatures  and  directed  velocities 
analyzed,  the  results  for  the  averaged  macroscopic  cross 
sections  and  for  the  inverse  of  the  average  neutron  velocity 
presented  a small  difference  (around  5%)  when  the  target  atom 
directed  velocity  was  increased  from  Vj,=0.0  to  1.0,  for 
temperatures  above  2000  °K,  indicating  that  here  the  directed 
velocity  effects  can  be  reasonably  neglected.  Outside  these 
limits  the  target  atoms  directed  velocity  should  be  considered 
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in  the  computations  for  the  averaged  macroscopic  cross 
sections . 


6.4  Recommendations  for  Future  Work 

This  dissertation  answered  many  questions  about  the 
effects  of  high  target  atom  directed  velocity  and  temperature 
on  interaction  rates  for  a flowing,  fissioning  gas  for 
specific  conditions  but  it  was  not  intended  to  address  all  the 
questions  that  might  arise  from  this  extensive  and  complex 
topic.  In  the  development  of  this  limited  study,  many 
questions  arose  leading  to  a series  of  proposed  extensions  to 
this  dissertation  as  follows: 

1)  Analyze  the  problem  to  include  low  (thermal)  energy 
resonances  effects  if  the  gaseous  fuel  includes  isotopes  that 
present  resonance  cross  section  in  the  low  energy  range.  This 
analysis  might  be  accomplished  by  including,  after  suitable 
modifications,  an  auxiliary  program,  like  DOPCRS  [35],  that 
performs  Doppler  broadening  of  fission,  absorption  and 
scattering  cross  sections  for  resonance  isotopes,  in  BRT-PCM. 

2)  Analyze  the  behavior  of  the  neutron  flux  and  averaged 
macroscopic  cross  sections  as  the  number  densities  of  the  more 
significant  isotopes  are  changed  significantly  from  those  used 
in  this  study.  This  analysis  is  important  to  extend  the  work 
to  include  systems  which  are  different  from  the  system 
modeled.  For  example,  some  gas  core  reactors  used  in  certain 
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propulsion  system  concepts,  have  much  higher  hydrogen  number 
densities  (by  up  to  two  orders  of  magnitude)  than  the  system 
modeled  in  this  work.  This  analysis  should  be  performed  for  a 
number  of  fixed  temperatures  in  order  to  observe  the  effects 
of  the  directed  velocity  as  the  number  density  is  changed. 
This  can  be  readily  studied  using  BRT-PCM. 

3)  Analyze  a system  that  doesn't  present  a flat  velocity 
profile.  The  assumption  of  a flat  profile  is  reasonable  for  a 
very  low  density  gaseous  system  flowing  at  high  velocity.  One 
suggested  approach  to  handle  the  problem  if  the  condition  of 
a flat  velocity  profile  can't  be  met  is  to  employ  Monte  Carlo 
methods.  This  can  turn  out  to  be  particularly  important  in  the 
boundary  layer  region. 

4)  The  analysis  presented  is  for  the  two  extreme  cases  of 
uncollided  and  completely  entrained  neutrons.  It  is  expected 
that  the  thermal  neutron  flux  and  the  effective  cross  sections 
for  the  intermediate  cases  between  the  uncollided  and 
completely  entrained  states  will  assume  intermediate 
distributions/values  between  those  for  these  two  extreme 
cases.  A complete  analysis  for  these  intermediate  cases  should 
be  performed  which  includes  the  transfer  of  momentum  between 
the  neutrons  and  target  atoms.  This  would  probably  require  the 
use  of  Monte  Carlo  methods.  Alternatively,  a suggested 
approximate  method  to  be  developed  to  analyze  some  but  not  all 
of  these  intermediate  conditions  is;  1)  for  a given  directed 
velocity,  V^,  divide  into  a set  of  discrete  intervals  AV^; 
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2)  start  with  an  analysis  of  the  uncollided  condition;  3)  then 
for  this  small  AV^  analyze  the  completely  entrained  condition; 

4)  assuming  that  an  equilibrium  condition  is  reached  at  AV^,, 
use  the  data  obtained  in  the  preceding  step  to  repeat  the 
calculation  for  the  uncollided  condition.  This  process  is 
repeated  as  many  times  as  necessary  until  the  sum  of  AV„ 
equals  Actually  what  is  simulated  is  a continuous 
transition  between  a series  of  two  extreme  states  that  are  AVj^ 
apart  (and  which  are  each  assumed  to  be  quasi-equilibrium 
states)  until  the  final  completely  entrained  case  at  V,^  is 
reached.  The  above  approximate  method  can  be  applied  using  the 
approach  developed  in  this  work. 

5)  At  this  time  there  are  no  experimental  data  available 
against  which  to  compare  the  predicted  results.  It  is 
recommended  that  experiments  be  developed  to  obtain  data 
against  which  the  predictions  can  be  compared. 

6)  The  analyses  performed  in  this  work  were  all  done  in 
cylindrical  geometry.  Nevertheless,  the  conclusions  related  to 
the  target  atom  directed  velocity  effects  are  valid  for  any 
system  possessing  a similar  range  of  directed  velocity  and 
temperature.  The  numerical  results  for  specific  systems 
involving  different  geometries  and  isotopes  for  various 
combinations  of  high  target  atom  directed  velocity  and 
temperature  can  be  obtained  using  BRT-PCM  after  modifications 
are  made  to  include  new  subroutines  in  BRT-PCM  to  allow  it  to 
handle  these  additional  geometries. 


APPENDIX  A 
BRT-I 

A.l  Introduction 

The  BATTELLE  REVISED  THERMOS  computer  code  (BRT-I)  [33] 
is  a revised  version  of  THERMOS  which  includes  many 
modifications  to  the  original  THERMOS  [28]  computer  code, 
developed  by  Henry  C.  Honeck  in  1961.  The  THERMOS  code 
computes  the  scalar  neutron  spectrum  as  a function  of  position 
in  a one  dimensional  slab  or  cylindrical  lattice  by  solving 
numerically  the  integral  transport  equation  with  isotropic 
scattering. 


A. 2 Brief  Review  of  THERMOS  Code 


The 

purpose 

of  THERMOS 

is  to 

calculate  the 

thermal 

neutron 

flux  in 

a reactor 

lattice 

with  emphasis 

on  the 

accurate  prediction  of  thermal  neutron  reaction  rates.  The 
quantity  of  interest  is  then  the  scalar  neutron  density, 
N(r,v) . 

THERMOS  is  written  entirely  in  FORTRAN  II  and  is  composed 
of  a main  program,  XPREP,  and  the  major  subroutines  GEOM, 
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ITER,  and  EDIT.  There  are  two  GEOM  subroutines,  one  for  slab 
and  the  other  for  cylindrical  geometry.  A flow  diagram  for 
THERMOS  is  shown  in  Fig.  A.l.  This  diagram  is  reproduced  from 
the  original  THERMOS  manual  [28]  and  in  the  diagram,  solid 
lines  denote  flow  of  control  and  dashed  lines  denote  flow  of 
data. 

The  code  assumes  that  the  emission  (birth  rate)  density 
H(r,v,n)=  / dv'  J dD.'  P(r,v,  v' , Aig)  N(r,v',H)  (A-1) 

is  isotropic  and  can  be  written  as 

H(r,v,)=  / dv'  P(r,v,v')  N(r,v').  (A-2) 

The  justification  for  this  assumption  is  that  the  birth 
rate  density  will  be  nearly  isotropic  if  either  the  scattering 
process  is  nearly  isotropic  or  the  neutron  density  is  nearly 
isotropic.  Both  conditions  can  be  considered  reasonable  if  it 
is  recognized  that  the  scattering  kernel  for  thermal  neutrons 
is  more  isotropic  than  for  high  energy  neutrons  due  to  thermal 
motion  and  binding,  and  the  gradient  of  the  flux  in  a lattice 
is  generally  small. 

The  basic  neutron  balance  equation  is  written  as 
N(r,v)=  H' (r*  ,v)T(r,r' ,v)dr'  (A-3) 

where  N(r,v)dv  is  the  neutron  number  density  at  position  r 
having  speed  v in  dv,  and  H'(r',v)  is  (1/v)  (H(r ' ,v) ) , where 
H(r',v)dv  is  the  birth  rate  density  of  neutrons  at  position  r' 
with  speed  v in  dv. 

THERMOS  uses  a transport  point  kernel 
Tpt(r,r' ,v)  = (l/47T(r-r' )2)exp{-  I J^/ds  S*(s,v)l  } 


(A-4) 
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where  Tpj(r,r',v)  represents  the  flux  of  uncollided  neutrons 
at  r produced  by  a unit  source  of  neutrons  of  speed  v located 
at  r'.  The  speed  range  is  divided  into  an  epithermal  region 
where  both  the  energy  and  spatial  distribution  of  neutrons  are 
known,  and  a thermal  region.  The  boundary  between  these  two 
regions,  v*,  can  be  chosen  by  the  user  and  has  0.78  eV  as  the 
default  value  . 

The  slowing  down  source  used  by  THERMOS  is 


and  the  computer  code  solves  numerically  the  integral  equation 
for  the  neutron  density  defined  as 


where  the  primes  on  S and  P denotes  division  by  v. 

The  ability  to  use  any  tabulated  scattering  kernel  is  one 
of  the  features  of  the  code.  The  free  gas  kernel  is  so 
frequently  used  and  easy  to  compute  that  the  Brown  and  St. 
John  form  is  built  into  both  the  library  tape  program  and 
THERMOS  itself.  This  modified  free  gas  kernel  uses  the  Brown 
and  St.  John  approximation  for  the  scattering  cross  section 


where  v^.  is  the  relative  speed  between  the  neutron  and  the 
target  atom  and  and  x^  are  parameters  adjusted  such  that 
S®(v)  fits  the  measured  scattering  cross  section.  The 
scattering  kernel  obtained  from  Equation  (A-7)  is  given  by 


S(r',v)=  dv'P(r ' ,v,v' )N(r ' ,v) 


(A-5) 


N(r,v)=  dr'  T(r,r',v)  [S'(r',v)  + 

Jp  dv'P'(r',v,v')  N(r',v')] 


(A-6) 


(A-7) 
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P(V,V')=  E Xn'  (V/V)  ^ 

{exp(-xJ^2^'2)E^  + exp[  (B^Vm,)  (v'2- 

v'(B^2  + x^(M^+1))/(B^2  + xJ)]E2)  I 


(A-8) 


The  scattering  kernel  shown  in  Equation  (A-8) , integrated  over 
v'  yields: 


M^=  atomic  mass  (amu) 

The  code  uses  T as  the  material  temperature  in  units  of 
the  reference  temperature  Tq=293.6°K. 

The  expression  for  the  scattering  kernel  is  valid  for  v' 
greater  than  or  equal  to  v;  the  rest  of  the  kernel  is  computed 
from  detailed  balance. 

The  purpose  of  this  brief  presentation  on  THERMOS  is  to 
show  its  main  features  and  structure  since  BRT-1  originated 


(A-9) 


where 


from  THERMOS. 
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A. 3 Brief  Review  of  BRT-I  Code 


Like  THERMOS,  BRT-I  solves  the  integral  neutron  transport 
equation.  BRT-I  computes  the  space  dependent  thermal  neutron 
density,  flux  and  current  spectra  over  an  arbitrary  thermal 
energy  range  (the  current  default  range  is  from  0.0  eV  to 
0.683  eV)  in  either  slab  or  cylindrical  geometry. 

In  this  code,  the  neutron  density  is  computed  from  the 
collision  probability  form  of  the  integral  transport  theory 
iTiatrix  equation  using  either  a combination  of  power  iteration, 
overrelaxation  and  extrapolation  or  straight  power  iteration. 
The  neutron  currents  are  computed  from  either  the  gradient  of 
the  scalar  flux  or  the  uncollided  flux  matrix.  The  flux  and 
current  spectra  are  used  to  weight  point  thermal  cross 
sections  over  an  arbitrary  thermal  energy  range  for  use  in 
niultigroup  transport  or  diffusion  theory  codes. 

Although  BRT-I  is  based  on  THERMOS,  many  improvements 
where  made  to  the  original  code.  Important  revisions  include: 

a)  An  improved  library  tape  which  contains  the  first 
moment  scatter  matrices  for  a number  of  elements  and  contains 
the  free  atom  scattering  cross  section  for  all  elements. 

b)  Transport  kernel  modifications  to  include  the 
treatment  of  void  regions. 

c)  Addition  of  the  cosine  current  calculation  of  the 
transport  kernel  to  allow  for  the  input  of  a right  hand 


albedo. 
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d)  The  ability  to  handle  up  to  8 materials  mixtures 
(THERMOS  was  restricted  to  only  5 mixtures) , 30  space  points 
in  both  slab  and  cylinder  geometry  (THERMOS  allows  only  20 
space  points)  , and  30  isotopes  (THERMOS  can  handle  only  10 
isotopes) . 

e)  An  option  that  allows  one  to  use  an  approximate  method 
of  correcting  for  linear  anisotropic  scattering.  The 
anisotropic  correction  to  the  scattering  kernel  is  applied 
only  to  the  diagonal  terms  of  the  scattering  matrices  (for 
tabulated  materials) . For  materials  without  a first  moment 
scattering  kernel  (i.  e.  for  manually  input  data) , the 
expression 

a®^(v)  = (2/3A)a®°(v) 

is  used.  Actually  2/3A  is  the  value  of  the  mean  cosine  of 
neutron  scattering  in  the  laboratory  system  for  isotropic 
scattering  in  the  center  of  mass  system. 

f)  An  option  that  allows  the  user  to  obtain  multiple 
energy  group  structure  edits,  with  a maximum  of  five  broad 
group  edits  over  the  energy  range. 

g)  Greater  variety  of  boundary  conditions  are  available 
in  BRT-I  in  lieu  of  THERMOS'  restriction  to  only  vacuum  or 
reflecting  boundaries. 

h)  Inclusion  of  neutron  current  calculation  routines  and 
the  use  of  neutron  current  to  weight  the  transport  and  first 
moment  scatter  cross  sections. 
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i)  Improvement  in  the  calculation  of  the  diffusion 
coefficient  due  to  the  use  of  improved  neutron  current 
calculation. 

Beyond  the  options  already  mentioned,  BRT-I  also  offers 
the  user  the  ability  to  input  a desired  speed  mesh  and  to 
input  data  for  isotopes  not  contained  in  the  BRT-I  library; 
for  these  isotopes  the  user  can  input  tabular  cross  sections 
and  scattering  kernel  or  use  the  built-in  free  gas  model  that 
includes  the  temperature  for  each  isotope;  slowing  down  source 
data  can  also  be  input  for  these  added  isotopes. 

The  output  of  BRT-I  is  divided  into  three  main  sections. 
Distributions  and  averages  over  both  space  and  energy  are 
listed  in  the  first  section.  In  the  second  section, 
microscopic  and  macroscopic  averages  are  listed  for  the 
spatial  edit  regions  specifically  defined  in  the  input, 
independent  of  the  physical  region. 

The  third  section  is  mainly  for  editing  purposes,  and 
lists  point  and  region  averaged  absorption  and  fission  cross 
sections  and  reaction  rates. 

The  output  of  BRT-I  is  very  detailed.  It  presents,  first, 
the  most  important  data  from  the  input,  and  then  the  results 
of  the  computations.  The  most  relevant  data  presented  in  the 
output  are: 

a)  Tables  of  neutron  density  and  neutron  current  for  each 
spatial  point  and  velocity  point. 
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b)  For  each  velocity  interval  associated  with  each 
velocity  point,  and  each  region,  the  output  shows  the 
absorption,  fission  and  scattering  interaction  rates 
integrated  over  all  spatial  points  in  the  region;  these 
interaction  rates  are  also  shown  integrated  over  all 
velocities.  The  volume-averaged  values  for  each  region  are 
also  presented. 

c)  For  each  spatial  point,  the  output  shows  the  neutron 
number  density,  neutron  flux,  absorption,  fission  and 
scattering  interaction  rates  integrated  over  all  velocities, 
niultiplied  by  the  volume  associated  with  that  point  and  also 
averaged  over  the  total  volume  of  the  region. 

d)  The  microscopic  constants  for  each  isotope  are  given 
according  to  the  region  and  mixture.  The  smeared  values  for 
each  region  and  the  smeared  values  for  the  ensemble  of  all 
regions  (cell)  are  also  presented. 

e)  Finally,  macroscopic  constants  like  inverse  velocity, 
diffusion  coefficient,  and  scattering,  absorption,  and  fission 
cross  sections  are  presented  for  each  region  and  for  the 
"cell". 

The  flow  diagram  for  BRT-I  is  shown  in  Fig.  A. 2 where  the 
dashed  lines  indicate  conditional  calls  to  subroutines,  the 
solid  lines  demonstrate  the  unrestricted  normal  flow  of  the 
program,  and  the  numbers  show  the  sequence  in  which  the 
subroutines  are  called. 
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Figure  A. 2 - Schematic  Flow  Diagram  of  BRT-I. 


APPENDIX  B 

PRELIMINARY  THERMAL  SPECTRUM  STUDIES 
USING  THE  COMPUTER  CODE  BRT-I 


B.l  Introduction 


This  early  computer  work  was  performed  with  two  main 
objectives:  (1)  to  select  and  become  acquainted  with  a 
computer  code  which,  we  expect,  can  be  modified  to  include  the 
new  or  modified  thermal  energy  cross  section  formulations 
developed  as  a result  of  this  research  and  (2)  to  perform 
analysis  with  the  unmodified  version  of  this  code  to  observe 
the  thermal  spectrum  and  cross  section  behavior  as  selected 
parameters  are  changed  and  also  to  test  the  options  offered  by 
the  code. 

The  chosen  code  is  the  BATTELLE  REVISED  THERMOS  (BRT-I) 
[28,33]  (see  App.  A),  which  computes  the  space  dependent 
thermal  energy  neutron  density,  flux,  and  current  spectra  over 
an  arbitrary  thermal  energy  range  (the  current  default  range 
is  from  0.0  eV  to  0.683  eV)  in  either  slab  or  cylindrical 
geometry. 

Many  calculations  were  performed  with  the  unmodified 
version  of  the  code  on  a generic  gas  core  (cavity)  reactor, 
with  cylindrical  geometry,  to  examine  the  effects  on  neutron 
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temperature,  spectrum,  and  cell  averaged  cross  sections,  of 
changes  in  selected  parameters.  The  chosen  system  has  a core 
with  a 3 0 cm  radius  and  a surrounding  moderator/reflector 
thickness  of  55  cm.  The  selected  parameters  were  pressure  (or 
number  density,  if  the  gas  temperature  is  kept  constant)  of 
the  hydrogen;  temperature  of  the  gas;  and  temperature  of  the 
external  moderator. 

The  code  has  a built-in  option  that  permits  the  use  of 
the  Brown  and  St.  John  (BSJ)  free  gas  kernel  approximation  as 
a way  of  entering  the  scattering  cross  section  and  scattering 
kernel  data  for  an  isotope  (see  App.  A)  . The  parameters  in  the 
BSJ  approximation  are  obtained  by  performing  fits  for  each 
isotope  and  energy  range  to  the  scattering  cross  section  data 
for  the  isotope.  When  this  option  is  used,  all  the  other  cross 
section  data  for  this  isotope  must  be  read  in. 

The  standard  BRT-I  library  contains  cross  section  data  at 
specific  temperatures.  Anticipating  the  need  of  working  at 
temperatures  that  are  much  higher  than  those  included  in  this 
library,  some  cases  were  run  to  test  the  above  option.  The 
purpose  was  to  verify  how  results  obtained  with  the  use  of  the 
BSJ  approximation  agree  with  results  obtained  when  the  cross 
section  data  for  the  isotope  were  acquired  directly  from  the 
standard  BRT-I  library.  The  scattering  kernels  on  this  library 
were  generated  with  the  FLANGE  code  which  processes  thermal 
neutron  cross  section  data,  including  resonance  parameters, 
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angular  distributions,  and  S(a,^,T)  scattering  laws  from  the 
ENDF/B  files. 

The  gas  core  reactor  thermal  spectrum  is  expected  to 
extend  to  energies  well  above  the  current  BRT-I  library 
default  upper  limit  of  0.683  eV;  hence,  an  option  in  the  code 
was  used  which  permits  user  specification  of  the  energy 
limits.  Initially,  the  default  energy  range  was  extended  to 
1.83  eV.  To  do  this,  the  new  energy  structure  and  all  the 
required  cross  section  data  had  to  be  read  in  as  input. 

The  initial  extension  of  the  upper  energy  limit  to  1.83 
eV  was  selected  because  this  is  the  upper  energy  limit  of  the 
thermal  region  in  the  AMPEX  Master  Library  which  is  available 
at  the  University  of  Florida.  The  AMPEX  package  contains  a 
module  called  NITAWL  [34]  which  is  used  to  prepare  working 
libraries  from  the  AMPEX  Master  Library.  Cross  section  data 
from  a NITAWL  generated  working  library  were  used  as  input  to 
BRT  when  the  energy  limit  was  extended  to  1.83  eV.  For  the 
scattering  kernel,  the  BSJ  free  gas  kernel  approximation  is 
used  as  input  and  its  parameters  are  obtained  by  performing 
fits  to  the  NITAWL  library  scattering  cross  section  data,  over 
the  new  energy  range,  for  each  isotope. 

The  results  obtained  using  this  new  energy  range  and  the 
NITAWL  library  data  were  compared  to  the  results  obtained 
before,  with  the  default  energy  limit  and  the  BRT-I  library 
data,  for  the  same  configuration  under  the  same  conditions. 
The  comparison,  however,  is  valid  only  if  the  same  energy 
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range  is  used.  There  is  an  option  in  BRT-I  which  permits  the 
editing  of  the  results  in  energy  over  any  subrange  that  is 
contained  within  the  input  or  default  energy  range  (see  App. 
A).  Use  of  this  option  to  edit  results  for  a 0.65  eV  upper 
energy  limit,  in  the  cases  when  the  energy  limit  was  extended 
to  1.8  3 eV,  allowed  us  to  perform  the  above  mentioned 
comparison. 

After  verifying  that  the  thermal  energy  range  could 
reliably  be  extended  using  the  BSJ  free  gas  approximation  and 
the  read  in  cross  section  data,  the  code  was  run  again, 
changing  the  temperature  of  the  gases  in  the  core,  first  for 
hydrogen  and  then  for  the  other  "gases",  to  examine  the 
resulting  changes  in  neutron  temperature,  spectrum,  and  cell 
averaged  cross  sections. 

For  the  preliminary  calculations,  materials  and  number 
densities  that  are  typically  encountered  in  gaseous  core 
reactors  were  employed.  The  gas  in  the  core  can  include 
hydrogen,  uranium  235,  uranium  238,  helium,  and  fluorine,  and 
the  reflector/moderating  medium  is  beryllium  or  beryllium 
oxide. 

To  keep  the  work  organized,  the  runs  were  labeled 
"BRT.n",  where  n is  the  number  of  each  run. 
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B.2  Description  and  Results 

1)  In  the  first  four  cases  the  effects  of  increasing  the 
pressure  (or  number  density,  if  the  temperature  is  kept 
constant)  of  hydrogen  in  the  core  are  examined  (in  these 
cases,  deuterium  was  inadvertently  used  instead  of  hydrogen) . 

The  temperature  (°K)  and  nuiaber  densities  of  the 
materials  (except  deuterium)  are  presented  in  Table  B.l.  It 
should  be  noticed  that,  in  most  cases,  the  highest 
temperatures  available  for  each  material  in  the  current  BRT-I 
library  were  used. 

Table  B.l  - Temperature  and  Number  Density  of 
Materials . 


NUCLIDE 

T 

N 

(°K) 

(atoms/b. cm) 

Core 

uranium  235 

290 

4.283  E-5 

uranium  238 

290 

3.224  E-6 

helium 

1200 

6.764  E-4 

fluorine 

1200 

2.763  E-4 

Moderator 

beryllium 

1000 

1.200  E-1 

A brief  description  of  the  cases  is  as  follows: 

- BRT.l:  In  this  case,  hydrogen  (deuterium)  was  not 

included  in  the  core.  This  is  a reference  case  against  which 
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results  from  other  cases,  containing  hydrogen  (deuterium)  are 
compared. 

- BRT. 2/3/4: In  these  cases  hydrogen  (deuterium)  at  570°K 
was  used  and  the  number  densities  of  hydrogen  (deuterium) 
change  according  to  the  pressure. 

The  hydrogen  pressures  and  atom  densities  for  the  above 
cases  are  listed  in  Table  B.2  and  the  results  are  presented  in 
Table  B.3,  where: 

M:  Moderator  Region 

C:  Core  Region 

Tj^:  Neutron  Temperature  (°K) 

E^:  Neutron  Energy  at  the  Average  Velocity  (eV) . 

The  above  symbols  and  the  units  used  in  the  Table  B.3  are  also 
used  in  the  remaining  Tables  in  this  Appendix  B. 


Table  B.2  - Hydrogen  Pressures  and  Number  Densities. 


CASE 

Hydrogen  Pressure 
(atm) 

Hydrogen  Atom 
Density (atoms/b. cm) 

BRT.  1 

0.000 

0.000 

BRT.  2 

1.666 

2.142  E-5 

BRT.  3 

10.000 

1.287  E-4 

BRT.  4 

100.000 

1.287  E-3 

Conclusions: 
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a)  From  the  results  in  Table  B.3,  it  is  determined  that, 
initially,  due  to  the  increased  absorption  by  hydrogen 
(deuterium)  addition  there  is  a hardening  of  the  core 
spectrum.  However  there  is  a point  after  which  the  hardening 
tendency  is  reversed;  with  further  hydrogen  (deuterium) 
additions  the  scattering  effect  begins  to  dominate  and  there 
is  rethermalization  in  the  core  to  lower  energies  because  the 


Table  B.3  - Results  of  Cases  BRT.l,  BRT.2,  BRT.3,  BRT.4. 


VARIABLE 

REGION 

BRT.l 
(no  H) 

BRT.2 
(1.66  atm) 

BRT.3 
(10  atm) 

BRT.4 
(100  atm) 

S®  (cm''') 

C 

1.317  E-2 

1.117  E-2 

1.118  E-2 

1.131  E-2 

S®  (cm'') 

M 

6.709  E-4 

5.670  E-4 

5.670  E-4 

5.671  E-4 

(cm'') 

C 

1.114  E-2 

9.416  E-3 

9.426  E-3 

9.530  E-3 

vS^(cm'') 

C 

2.707  E-2 

2.288  E-2 

2.291  E-2 

2.316  E-2 

S®  (cm'') 

c 

2.085  E-3 

2.152  E-3 

2.538  E-3 

6.724  E-3 

E®  (cm'') 

M 

7.584  E-1 

7.489  E-1 

7.489  E-1 

7.489  E-1 

T„  (°K) 

C 

933 

1275 

1273 

1250 

E„  (eV) 

C 

0.1020 

0.1399 

0.1396 

0.1370 

hydrogen  (deuterium)  is  at  lower  temperature  (570°K)  than  the 
beryllium  (1000°K)  . 

b)  The  effects  on  the  moderator  parameters  were  notable 
after  the  initial  introduction  of  hydrogen  (deuterium)  was 
made,  but  further  addition  of  hydrogen  (deuterium)  had  a much 
smaller  effect  on  the  moderator  than  on  the  core. 
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2)  The  next  three  cases  used  hydrogen  instead  of 
deuterium,  at  the  same  temperature  and  pressures  as  above;  the 
beryllium  is  now  at  290°K.  The  other  materials  have  the  same 
number  densities  and  temperatures  as  for  Item  1. 

The  results  are  presented  in  Table  B.4. 

- Conclusion: 

The  tendency  is  that  the  changes  in  the  core  constants 
occur  continuously  in  the  same  direction  as  the  number  density 
of  hydrogen  is  increased.  This  is  because  rethermalization  in 
the  core  is  now  to  higher  energies,  since  the  hydrogen  is  at 
a higher  temperature  (570°K)  than  the  beryllium  (290°K)  . 


Table  B.4  - Results  of  Cases  BRT.5,  BRT.6,  BRT.7. 


BRT.5 

BRT.6 

BRT.7 

VARIABLE 

REGION 

(no  H) 

(10  atm) 

(100  atm) 

C 

1.898  E-2 

1.881  E-2 

1.789  E-2 

M 

9.520  E-4 

9.522  E-4 

9.535  E-4 

S® 

C 

2.120  E-3 

7.280  E-3 

5.108  E-2 

Z® 

M 

6.977  E-1 

6.977  E-1 

6.977  E-1 

C 

485 

495 

555 

C 

0.0532 

0.0543 

0.0609 

3)  BRT.  8/9/10  have  the  same  data  as  Item  2 but  the 
beryllium  was  replaced  by  beryllium  oxide.  The  temperatures 
remained  the  same  as  for  Item  2 . 


The  number  densities  for  BeO  are  : 
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- beryllium:  7.0  E-2  (atoms/b.cm) 

- oxygen:  7.0  E-2  (atoms/b.cm) 

The  results  are  given  in  Table  B.5. 

- Conclusions:  Looking  at  the  results  it  is  determined 
that,  in  these  cases,  the  behavior  of  the  constants  is  quite 
similar  to  Item  2.  Again,  there  are  very  small  changes  in  the 
moderator  constants,  and  the  tendency  of  changes  in  the  core 
parameters,  as  the  number  density  of  hydrogen  is  increased,  is 
not  reversed. 


Table  B.5  - Results  of  Cases  BRT.8,  BRT.9,  BRT.IO. 


VARIABLE 

REGION 

BRT.8 
(no  H) 

BRT.9 
(10  atm) 

BRT. 10 
(100  atm) 

E® 

C 

1.890  E-2 

1.877  E-2 

1.809  E-2 

E® 

M 

6.185  E-4 

6.188  E-4 

6.196  E-4 

E^ 

C 

1.607  E-2 

1.593  E-2 

1.513  E-2 

vE^ 

C 

3.906  E-2 

3.872  E-2 

3.676  E-2 

E® 

c 

2.118  E-3 

7.288  E-3 

5.156  E-2 

E® 

M 

9.708  E-1 

9.708  E-1 

9.708  E-1 

T. 

C 

489 

496 

543 

En 

C 

0.0536 

0.0544 

0.0596 

4)  The  use  of  BRT-I  is  usually  restricted  to  the 
materials  and  temperatures  from  its  library,  but  BRT-I  offers 
an  option  by  which  the  user  can  enter  his  own  cross  section 
data  in  a tabular  form.  This  includes  the  ability  to  employ  a 
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built-in  Brown  and  St.  John  (BSJ)  free  gas  scattering  kernel 
approximation . 

To  test  this  option,  above  the  BRT.IO  run  is  taken  as  a 
benchmark  and  an  attempt  is  made  to  duplicate  the  results 
obtained  in  Item  3,  using  the  current  BRT-I  library  data  for 
all  the  materials,  except  hydrogen.  The  hydrogen  cross  section 
data  were  input  manually.  The  hydrogen  data  entered  were  in 
fact  actually  obtained  from  the  current  BRT-I  library  and  the 
hydrogen  scattering  cross  sections  from  the  BRT-I  library  were 
used  to  generate  the  parameters  needed  for  the  BSJ  free  gas 
scattering  kernel . The  code  then  uses  this  kernel  to 
calculate,  internally,  the  hydrogen  scattering  cross  sections. 
This  option  permits  the  user  to  designate  the  temperature  of 
the  input  isotope.  In  this  particular  case  the  temperatures 
and  number  densities  for  all  the  materials  were  kept  the  same 
as  in  BRT.IO  in  order  to  directly  compare  the  results. 

The  BRT-I  code  built-in  Brown  and  St.  John  free  gas 
scattering  kernel  approximation  is  described  by: 
a®(v^)=  A/exp(-k^*(v^)2)  + . . .+A/exp(-k^*(v^)2)  ^ 

In  the  above  expression: 

a®  = microscopic  scattering  cross  section 
Vj,  = relative  velocity 
Aj,kj  = fitting  parameters 

m = number  of  terms  used  in  the  approximation 
To  use  this  approximation  the  scattering  cross  section 
data  for  the  material  (obtained,  e.  g. , from  a cross  section 
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library)  are  fit  to  the  above  expression,  adjusting  the 
parameters  Aj  and  kj  so  as  to  optimize  the  fit.  This  can  be 
roughly  done  by  hand,  graphing  the  data  on  a log-log  scale  and 
making  a linear  fitting  that  results  in  only  one  term  in  the 
approximation  (m=l) . Alternatively,  if  more  terms  in  the 
approximation  are  desired,  a computer  program  can  be  employed 
to  determine  the  values  for  the  parameters  Aj  and  k,. . 

Cases  BRT.ll  to  BRT.15  were  used  to  test  many  different 
fittings  (m=l,  m=2 , m=3)  of  the  scattering  cross  section  data 
from  the  BRT-I  library  to  the  BSJ  free  gas  kernel 
approximation  to  determine  which  one  would  give  the  best 
agreement  with  the  results  obtained  when  using  hydrogen  data 
directly  from  the  BRT-I  library  (BRT.IO). 

It  must  be  mentioned  that  BRT-I  has  a built-in  slowing 
down  source  generator.  However,  when  the  data  for  the 
materials  are  read  in,  as  for  hydrogen  in  the  cases  now  been 
investigated,  parameters  needs  to  be  input  to  allow  BRT-I  to 
calculate  the  slowing  down  source  contribution  from  the  read- 
in  material.  In  the  case  of  hydrogen  this  is  particularly 
important  because  of  its  significant  slowing  down  (scattering) 
properties . 

A brief  description  of  the  cases  is: 

- BRT.ll:  (m=l)  In  this  case  the  parameters  A-  and  kj  were 
determined  by  making  a linear  fitting,  using  a graph  and  hand 
calculations.  Thirty  hydrogen  scattering  cross  section  data 
points,  from  the  BRT-I  library,  were  used. 
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- BRT.12:  (m=2)  In  this  case,  trying  to  achieve  a better 
fitting,  a computer  program  was  used  to  determine  the  Aj  and 
kj  values.  The  available  thirty  values  (for  thirty  different 
energies)  of  the  hydrogen  scattering  cross  section  from  the 
BRT-I  library  were  again  used. 

- BRT.13:  (m=2)  This  case  is  like  BRT.12  but  instead  of 
using  all  thirty  fine  group  data  points,  only  twelve  of  them 
were  used  to  make  the  fitting. 

- BRT.14:  (m=2)  This  is  like  the  above  case  but  a 
different  set  of  twelve  data  points  are  used  to  make  the 
fitting. 

- BRT.15:  (m=2)  This  case  was  run  mainly  to  compare 
results  obtained  from  the  use  of  free  hydrogen  cross  section 
with  those  obtained  by  using  data  for  hydrogen  in  H2O.  The 
BRT-I  library  data  for  hydrogen  is  for  hydrogen  in  water;  to 
have  more  precise  results  in  subsequent  work  (Item  5) , these 
BRT-I  library  data  were  scaled  using  the  hydrogen  data  from 
ENDF/B-IV  for  the  free  hydrogen  atom.  The  fitting  was  then 
done  to  the  scaled  free  hydrogen  data. 

- Results: 

The  results  of  BRT.15  are  quite  different,  as  expected, 
from  the  results  of  the  other  cases  because  free  hydrogen 
cross  section  data  are  different  from  the  hydrogen  in  H2O  data 
in  the  BRT-I  library.  The  other  results  (BRT.ll,  BRT.12, 
BRT.13,  BRT.14)  were  all  reasonable  and  the  ones  that 
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demonstrated  better  agreement  with  the  BRT.IO  benchmark 
results  are  presented  in  Table  B.6. 

- Conclusions:  From  the  above  it  can  be  verified  that  the 
results  obtained  from  the  case  when  the  hydrogen  data  are  read 
in  and  the  BSJ  free  gas  kernel  is  used  are  in  what  can  be 
considered  to  be  a good  agreement  (approximately  3%  error) 


Table  B.6  - Results  of  Cases  BRT.IO,  BRT.ll,  BRT.14. 


VARIABLE 

REGION 

BRT.IO 

(benchmark) 

BRT.ll 

(m=l) 

BRT.14 

(m=2) 

C 

1.809  E-2 

1.842  E-2 

1.842  E-2 

M 

6.196  E-4 

6.185  E-4 

6.187  E-4 

C 

1.513  E-2 

1.541  E-2 

1.540  E-2 

C 

3.676  E-2 

3.745  E-2 

3.743  E-2 

E® 

c 

5.156  E-2 

5.597  E-2 

5.099  E-2 

E® 

M 

9.708  E-1 

9.708  E-1 

9.708  E-1 

C 

543 

524 

525 

c 

0.0596 

0.0575 

0.0576 

with  the  results  obtained  when  the  current  BRT-I  library  data 
are  used  directly.  The  small  differences  (up  to  2 to  3%)  that 
appear  in  the  results  were  expected,  since  the  BSJ 
approximation  is  being  used.  Actually,  several  factors 
contribute  to  the  difference  in  the  results.  The  most 
significant  are  listed  below  : 

(1)  The  first  and  second  moment  scattering  matrices  are 
included  in  the  current  BRT-I  library.  With  the  BSJ  free  gas 
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approximation,  the  anisotropic  scattering  data  are 
approximated.  This  can  possibly  explain  the  softening  in  the 
neutron  spectrum  when  the  BSJ  approximation  is  used. 

(2)  The  current  BRT-I  library  uses  scattering  kernels 
generated  from  ENDF/B  S(q:,/3,T)  scattering  kernels.  When  the 
BSJ  free  gas  approximation  is  used,  the  scattering  kernel  is 
calculated  internally  using  the  fitting  parameters  obtained 
from  the  fitted  scattering  cross  section  data.  Although  the 
fittings  are  performed  from  cross  sections  which  are  based 
upon  S(a,)S,T)  data,  there  is  always  an  approximation  involved 
that  leads  to  small  differences  in  the  results. 

The  results  obtained  in  this  item  demonstrate  that  this 
approach  is  reliable  and  can  be  used  to  find  solutions  for 
other  problems.  It  should  be  noted  that  cross  section  data 
available  from  a library  other  than  the  BRT-I  library  can  be 
used  as  input  to  BRT-I . 

5)  Now,  after  testing  the  use  of  the  built-in  BSJ  free 
gas  kernel  in  Item  4,  BRT-I  is  used  for  fuel  gas  nuclide 
temperatures  beyond  the  limits  of  the  current  BRT-I  library. 

For  the  cases  in  this  item  instead  of  using  BeO,  as  in 
Item  4,  Be  is  used,  and,  since  very  high  temperatures  are 
expected  in  an  actual  gas  core  reactor,  the  temperature  of  the 
moderator  is  increased  to  900°K. 

Since  the  moderator  material  and  temperature  are  changed 
from  BeO  to  Be  and  from  290°K  to  900°K,  BRT.IO  cannot  be  used 
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as  a benchmark  and  another  benchmark  case,  BRT.17,  was 
established  using  all  cross  section  data  from  the  BRT-I 
library. 

Checking,  again,  the  reliability  of  the  BSJ  free  gas 
kernel , a case  was  run  where  the  data  for  the  hydrogen  are 
read  in  and  the  BSJ  free  gas  kernel  is  used.  The  results 
obtained  in  this  case  were  compared  with  the  results  of  the 
benchmark,  case  BRT.17.  Next,  cases  were  run  using  the  BSJ 
kernel  in  which  the  temperature  of  hydrogen  is  increased,  to 
observe  the  changes  that  occur  in  the  thermal  spectrum  and 
thermal  constants.  To  examine  how  the  changes  in  temperature 
of  the  other  gaseous  core  nuclides  affect  the  results,  the 
temperature  of  hydrogen  was  then  kept  constant  while  the 
temperature  for  the  other  "gases"  was  changed. 

The  number  densities  are  all  kept  as  in  Item  1 (see  Table 

B.  1)  . 

A brief  description  of  the  cases  follows: 

- BRT.17:  This  case  is  a benchmark;  the  cross  section 
data  for  all  the  isotopes  are  from  the  BRT-I  library.  The 
temperature  of  the  hydrogen  is  570°K. 

- BRT.18:  A BSJ  kernel  is  used  for  free  hydrogen  at 

570°K.  The  cross  section  data  for  all  other  isotopes  are  from 
the  BRT-I  library. 

- BRT.19:  Similar  to  case  18  but  the  temperature  of  the 
hydrogen  is  now  1000°K. 

- BRT.20:  Like  case  19,  but  with  hydrogen  at  2000°K. 


175 


- BRT.21:  The  data  for  all  the  gaseous  core  nuclides  are 
read  in  using  the  same  temperatures  as  for  Item  1 (see  Table 
B.l)  but  with  hydrogen  at  1000°K;  the  BSJ  free  gas  kernel  is 
used  for  all  the  "gases”. 

- BRT.22:  Now  all  the  "gases"  (including  hydrogen)  are  at 
1000°K  and  the  input  for  all  of  them  is  a BSJ  approximation. 
By  comparing  the  results  of  cases  21  and  22,  the  effects  of 
changing  the  temperatures  of  all  the  "gases"  except  hydrogen 
can  be  examined. 

- BRT.23:  Same  as  case  21  above,  but  with  all  the  "gases" 
at  2 000°K,  to  examine  how  the  thermal  spectrum  and  thermal 
constants  behave  when  all  "gases"  are  at  this  temperature. 

The  results  are  presented  in  Table  B.7. 

- Conclusions: 

a)  Comparing  cases  BRT.17  and  BRT.18  it  is  observed  that 
there  is  a reasonable  difference  in  the  results,  more  than  the 
difference  in  the  results  between  BRT.IO  and  BRT.ll  (up  to  4 
to  7%  versus  up  to  2 to  3%) . This  is  due  to  the  use  of  free 
hydrogen  data  in  BRT.18  and  not  a result  of  the  use  of  the  BSJ 
free  gas  kernel.  The  hydrogen  used  in  BRT.17  was  hydrogen  in 
H2O.  In  contrast,  both  BRT.IO  and  BRT.ll  employed  hydrogen  in 

H2O. 

b)  Cases  BRT.19  and  BRT.20,  as  expected,  show  us  that  as 
the  temperature  of  hydrogen  is  continuously  increased,  the 
neutron  spectrum  in  the  core  tends  to  be  hardened  due  to 
increased  rethermalization  to  higher  energies,  since  the 
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hydrogen  in  cases  BRT.19  and  BRT.20  is  at  higher  temperatures 
than  the  moderator. 

c)  Case  BRT.22,  when  compared  to  BRT.21  indicates  that 
changes  in  the  temperatures  of  the  "gases"  other  than  hydrogen 
(from  the  Table  B.l  values  to  1000°K)  don't  significantly 
affect  the  results. 

d)  The  results  of  cases  BRT.22  and  BRT.23,  when  compared 
with  cases  BRT.19  and  BRT.20,  demonstrate  very  good  agreement, 
indicating  that  the  use  of  the  built-in  BSJ  free  gas  kernel 
for  all  the  "gases"  is  reliable. 

6)  Work  is  now  started  with  an  energy  range  that  extends 
above  the  default  limit  set  by  the  current  BRT-I  library;  the 
default  upper  energy  limit  of  0.683  eV  changed  to  1.83  eV.  The 
new  energy  intervals  are  read  in  as  an  input,  and  because  of 
the  different  energy  limits,  the  current  BRT-I  library  cannot 
be  used;  all  the  cross  section  data  must  be  read  in. 

All  the  energy  intervals  and  cross  section  data  used  as 
an  input  to  BRT-I  and  also  used  to  determine  the  parameters 
for  the  built-in  BSJ  free  gas  kernel  are  obtained  from  the 
NITAWL  library. 

Changing  the  default  BRT-I  upper  energy  limit  requires 
that  all  the  cross  section  data  must  be  read  in.  For  the 
scattering  data  this  can  be  done  without  difficulty  for  all 
the  gases  by  means  of  the  BSJ  free  gas  kernel.  However,  a 
rigorous  treatment  for  beryllium  requires  the  manual  input  of 
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all  the  data  for  the  scattering  matrix  and  this  can  be  rather 
cumbersome . 

When  the  spectrum  obtained  using  a free  gas  model  is 
compared  with  the  spectrum  obtained  using  the  S(a,)S,T) 
scattering  law  data  for  beryllium  (or  beryllium  oxide) , it  is 
noticed  that  the  latter  presents  a spectrum  that  is  noticeably 
harder  than  the  free  gas  spectrum;  however,  the  free  gas  model 
becomes  more  accurate  for  higher  temperatures  [19],  and  for 
temperatures  high  enough  (how  high  depends  on  the  material) 
all  models  should  exhibit  gas-like  behavior.  This  fact  and 
also  the  knowledge  that  high  temperatures  will  be  worked  with 
throughout  this  research  indicates  that  for  beryllium,  use  of 
the  same  free  gas  treatment  as  employed  for  the  other  nuclides 
in  the  gas  core  may  be  a reasonable  approximation. 

As  mentioned  in  Item  B.l,  a particular  feature  of  BRT-I 
permits  the  edit  of  results  for  different  energy  intervals.  In 
this  case,  an  edit  of  the  results  over  the  energy  interval 
from  0 eV  to  0.65  eV  and  from  0 eV  to  1.83  eV  is  used  (in  the 
NITAWL  energy  structure,  0.65  eV  represents  the  closest  fit  to 
the  default  BRT-I  library  cutoff  of  0.683  eV) . This  allows  the 
results  obtained  in  the  0 eV  to  0.65  eV  energy  range,  using 
the  extended  upper  energy  limit  and  the  NITAWL  library  data, 
to  be  compared  with  those  of  a previous  case,  BRT.22,  where 
the  default  BRT-I  upper  energy  limit  of  0.683  eV  and  BRT-I 
library  data  are  used.  The  comparison  is  thus  performed 
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essentially  over  the  corresponding  energy  range  (0  to  0.65  eV 
versus  0 to  0.683  eV) . 

In  the  cases  studied,  the  temperature  of  all  the  gases  is 
1000°K  and  the  temperature  of  beryllium  is  kept  at  900°K. 

The  number  densities  are  the  same  as  in  Item  1 (see  Table 
B.l).  The  number  densities  of  the  gases  were  not  changed  from 
the  Table  B.l  values  to  correspond  to  1000°K  because  the 
intention  here  is  to  demonstrate  the  ability  to  extend  the 
BRT-I  energy  range. 

A brief  description  of  the  cases  is: 

- BRT.24:  In  this  case,  the  input  data  for  all  the 

materials,  including  Be,  are  read  in  and  the  BSJ  free  gas 
kernel  is  used.  The  data  used  to  obtain  BSJ  fits  for  Be  and 
all  the  other  materials  are  extracted  from  the  NITAWL  library. 

- BRT.25:  Like  case  BRT.24  for  the  gases,  but,  for  Be, 
the  microscopic  scattering  cross  section  data  are  read  in,  and 
an  option  in  the  code  is  used  that  assumes  P(V'-»  V)  = a®(V), 
if  V'=  V and  0 if  V is  different  from  V. 

The  results  are  shown  in  Table  B.8. 

- Conclusions: 

a)  Case  BRT.24  (0.65  eV) , had  all  the  cross  section  data 
and  energy  intervals  read  in  and  BSJ  free  gas  kernel 
approximations  were  used  for  all  the  materials.  Case  BRT.22 
(0.683  eV)  has  a very  similar  energy  range  but  used  the  cross 
section  data  for  beryllium  from  the  BRT-I  library  and  the 
built-in  energy  intervals.  Comparing  both,  it  is  observed  that 
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the  results  are  in  reasonable  agreement.  Thus  it  can  be  said 
that  the  use  of  the  new  energy  structure  with  read  in  cross 
section  data  and  the  use  of  the  BSJ  free  gas  kernel  for  Be  at 
900°K  are  reasonable. 


Table  B.8  - Results  of  Cases  BRT.24  (0.65  eV) , 
BRT.24  (1.83  eV) , BRT.25. 


VARIABLE 

REGION 

BRT.24 
(0.65  eV) 
(1000°K) 

BRT.24 
(1.83  eV) 
(1000°K) 

BRT.25 
(0.65  eV) 
(1000°K) 

C 

1.242  E-2 

1.176  E-2 

1.124  E-2 

M 

5.637  E-4 

5.533  E-4 

8.548  E-4 

C 

1.029  E-2 

9.752  E-3 

9.327  E-3 

C 

3.737  E-2 

3.633  E-2 

3.559  E-2 

E® 

M 

8.185  E-1 

8.160  E-1 

7.684  E-1 

C 

1111 

1211 

1317 

En 

C 

0.1218 

0.1329 

0.1444 

b)  Case  BRT.24  (1.83  eV)  shows  that  increasing  the 
thermal  energy  range  results  in  an  increase  in  the  average 
neutron  temperature  and  energy  in  the  core.  Actually,  this 
shows  that  there  is  significant  upscattering  above  0.65  eV 
which  can't  be  accounted  for  when  the  thermal  energy  cutoff  is 
at  0.65  eV  or  at  0.683  eV.  When  this  upscattering  above  0.65 
eV  is  accounted  for,  as  in  the  case  when  the  energy  cutoff  is 
1.83  eV,  the  neutron  spectrum  is  harder.  Similar  reasoning 
explains  why  the  BRT.22  results  (with  a cutoff  of  0.683  eV) 
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presents  a slightly  harder  spectrum  than  the  BRT.24  results 
with  a cutoff  of  0.65  eV. 

c)  Case  BRT.25  is  an  experiment  to  examine  an  option 
offered  by  the  code.  Comparing  its  results  with  the  results  of 
BRT.22  or  BRT.24  it  is  concluded  that  this  option  is  not 
suitable  for  the  present  case.  This  was  expected  since  the 
transfer  cross  section  data  for  Be  are  in  significant 
disagreement  with  the  assumption  presupposed  in  this  option. 


APPENDIX  C 

RELATED  EARLIER  WORK  BY  OBLOW 
C.l  Purpose 

Significant  work  related  to  the  proposed  research  was 
done  by  Oblow.  His  work  represented  an  effort  to  extend  the 
theory  of  neutron  thermalization  to  those  high  temperatures  at 
which  gas  fueled  nuclear  reactors  are  expected  to  operate. 
Major  emphasis  was  placed  on  determining  neutron  spectra  and 
neutron  cross  section  behavior  in  high  temperature  monoatomic 
gaseous  media. 


C.2  Methodology 

Assuming  an  infinite  homogeneous  medium  composed  of  a 
monoatomic  gas  in  thermodynamic  equilibrium  at  a specific 
temperature,  Oblow  used  the  unified  general  theory  of  neutron 
interaction  in  a monoatomic  gas,  developed  by  Blackshaw  [31], 
and  applied  the  concept  of  effective  cross  sections  to  the 
above  mentioned  theory. 

Having  laid  the  theoretical  foundations  for  a description 
of  high  temperature  gas  core  reactor  physics,  Oblow  developed 
suitable  numerical  methods  for  evaluating  the  derived 
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theoretical  expressions.  These  methods  together  with 
stationary  target  cross  section  data  can  then  be  used  to 
predict  the  behavior  of  neutrons  in  high  temperature  gaseous 
media . 


C.3  Results 


The  following  summarizes  key  results  and  conclusions  from 
Oblow's  work: 

a)  In  the  case  of  the  high  temperature  Doppler  effect,  a 
resonance  peak  shift  toward  lower  energies  was  observed  in  the 
effective  cross  section  behavior  of  low  energy  resonances  at 
high  temperatures  (0.296  eV  Pu^^’  resonance  at  70,000°K).  The 
conventional  Doppler  function  employing  the  Bethe-Placzek 
approximation  doesn't  show  this  shifting,  resulting  in  large 
errors  (maximum  error  of  the  order  of  50%  at  70,000°  K,  in  the 
low  energy  valley  below  the  resonance  peak)  in  the  predicted 
effective  cross  section  behavior. 

Oblow  concluded  that  the  shift  effect  is  due  to  target 
speeds  in  high  temperature  media  being  comparable  to  neutron 
speeds.  These  conditions  favor  a shift  in  any  effective  cross 
section  curve  to  lower  energies;  however,  in  the  case  of  high 
energy  resonances  (e.  g.  , 6.67  eV  resonance)  this  effect 
is  masked  by  Doppler  broadening. 
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b)  As  a result  of  the  approach  of  the  scattering  cross 
section  to  1/v  behavior  at  high  temperatures,  the  hydrogen 
effective  scattering  cross  section  is  significantly  increased. 

c)  The  study  of  the  elastic  transfer  process  with 
resonance  scattering  demonstrated  that  resonance  effects  are 
large  in  both  the  transfer  energy  distribution  and  the  average 
cosine  of  the  transfer  angular  distribution.  For  the 
particular  isotope  studied  (Pu^^°  with  a resonance  at  1.056  eV) 
the  resonance  was  found  to  increase  backscattering 
probabilities  in  the  angular  distributions  of  transfers  to  the 
energy  of  the  resonance  peak  and  forward  scattering 
probabilities  for  all  other  energy  transfers. 

d)  In  the  study  of  elastic  transfer  for  highly  subthermal 
neutrons  with  constant  scattering  cross  section,  it  was  shown 
that  neutrons  have  high  probabilities  for  transferring  over  a 
wide  energy  range  at  high  temperatures,  even  in  a heavy  gas. 

e)  The  behavior  of  neutrons  from  a monoenergetic 
subthermal  source  showed  that  in  a weakly  absorbing  medium 
subthermal  neutrons  completely  rethermalize  but  as  the 
absorption  increases  the  rethermalization  decreases  and,  at 
the  highest  level  of  absorption  tested,  rethermalization 
effects  were  limited  to  single  scattering  perturbations.  The 
introduction  of  resonance  scattering  and  absorption  is  seen  to 
cause  the  flux  to  increase  substantially  both  in  the  resonance 
and  on  its  high  energy  side.  The  inclusion  of  hydrogen 
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enhanced  the  effects  of  rethermalization  more  than  the 
inclusion  of  resonance  scattering. 

f)  In  the  case  of  rethermalization  from  a subthermal 
Maxwellian  source  distribution,  it  was  verified  that  the 
insertion  of  hot  hydrogen  causes  a larger  flux  perturbation 
(neutron  rethermalization) , over  the  whole  thermal  region, 
than  resonance  scattering. 

g)  Average  cross  sections  were  seen  to  be  very  sensitive 
to  both  Maxwellian  source  and  media  temperature  changes. 

C.4  Relationship  with  the  Present  Work 

As  discussed  above,  Oblow  based  his  theoretical 
development  of  exact  equations  related  to  neutrons 
interactions  with  target  atoms  submitted  to  very  high 
temperature  on  works  previously  done  by  Murray  and  Blackshaw. 
He  developed  numerical  solutions  to  solve  those  equations  and 
analyzed  the  results  obtained  by  using  the  numerical  solution 
for  selected  cases.  The  present  research  tends  to  follow  the 
guidelines  of  Oblow 's  work  but  includes  the  effects  of  high 
directed  velocity  of  the  target  atoms. 


APPENDIX  D 

THE  MODIFIED  BRT-I  COMPUTER  CODE 


D.l  BRT-PC 


The  code  selected  to  be  used  in  this  work,  BRT-I  (see 
Appendix  A) , is  written  in  FORTRAN-IV  and  runs  on  mainframe 
computer  systems . 

In  the  development  of  this  work  the  original  BRT-1  code 
was  used  extensively  (see  Appendix  B) . Following  the  approach 
developed  in  Chapter  3 , the  code  is  expected  to  be  used  even 
more  extensively,  in  its  modified  form,  in  the  analysis  of 
selected  cases  to  determine  the  effects  of  high  directed 
velocity  and  temperature  of  the  target  atoms  on  thermal 
neutron  interaction  rates.  In  order  to  facilitate  the 
implementation  of  these  required  modifications  and  also  in 
order  to  avoid  relying  on  a mainframe  system,  a personal 
computer  version  of  BRT-I  was  created. 

The  personal  computer  adapted  version  of  BRT-I,  herein 
called  BRT-PC,  requires  the  manual  input  of  all  data  and  in 
its  present  form  is  capable  of  solving  problems  only  in 
cylindrical  geometry.  Except  for  the  limitation  or  restriction 
to  cylindrical  geometry,  the  structure  and  the  capabilities  of 
the  original  BRT-I  code  were  kept  intact. 
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D.2  BRT-PCM 


The  solution  method  developed  in  Chapter  3 to  handle  the 
effects  of  target  atom  directed  velocity  on  thermal 
interaction  rates  requires  the  use  of  BRT-I  to  first  solve  the 
integral  neutron  transport  equation  in  the  moving  frame  and 
then  transform  the  results  back  to  the  static  frame. 

Since  what  is  desired  here  are  the  data  relative  to  the 
static  frame,  suitable  modifications  need  to  be  made  in  BRT- 
PC.  These  modifications  include  changes  in  the  input  file, 
inclusion  of  the  directed  velocity  in  the  code  input,  changes 
in  some  subroutines,  and  also  the  addition  of  subroutines  to 
change  the  results  from  the  moving  frame  back  to  the  static 
frame.  These  modifications  in  the  code  led  to  the  creation  of 
a new  code  called  BRT-PCM.  This  computer  code  is  very 
sensitive  to  the  input  data  and  the  physics  aspects  involved 
in  the  particular  problem  to  be  solved.  Information  contained 
in  the  rest  of  this  Appendix  will  help  to  indicate  which  input 
data  and  what  values  are  appropriate  to  be  used  in  BRT-PCM. 
The  schematic  flow  diagram  of  BRT-PCM  is  shown  in  Figure  D.l. 

As  stated  in  Chapter  3 , the  developed  solution  method  is 
suited  to  handle  two  different  extreme  cases:  uncollided 
neutrons  and  completely  entrained  neutrons.  These  two 
different  conditions  dictate  the  need  for  different 
treatments.  For  the  completely  entrained  neutrons  condition 
this  is  accomplished  in  BRT-PCM  by  the  subroutine  DIRVE.  For 
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Figure  D.l  - Schematic  Flow  Diagram  of  BRT-PCM. 
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the  uncollided  neutron  condition,  BRT-PCM  can  be  used  to 
generate  the  data  at  zero  target  atom  atom  directed  velocity. 
These  data  are  then  used  in  the  computations  for  the  data 
(averaged  macroscopic  cross  sections)  when  the  target  atom 
directed  velocity  effects  are  included. 

The  main  differences  between  BRT-PC  and  BRT-PCM  are 
outlined  below: 

a)  Input  file. 

The  main  differences  between  the  input  file  of  BRT-PC  and 
the  input  file  of  BRT-PCM  are  as  follows: 
a.l)  Speed  mesh 

Although  BRT-PC  uses  the  option  contained  in  BRT-I  that 
permits  a user-defined  speed  mesh,  this  speed  mesh  is  usually 
obtained  from  another  code,  like  NITAWL,  which  provides  cross 
sections  as  a function  of  speed  and  no  restrictions  are  placed 
on  the  choice  of  the  speed  mesh  upper  and  lower  limits.  The 
inclusion  of  the  directed  velocity  of  the  target  atoms  puts 
some  restrictions  on  the  speed  mesh  upper  and  lower  limits  to 
be  used  as  input  data  for  BRT-PCM. 
a. 2)  Cross  sections 

BRT-PC  uses  pointwise  cross  section  input  data  extracted 
from  an  available  library  according  to  the  speeds  in  the 
library.  Since  the  cross  sections  are  speed-dependent,  they 
must  be  selected  according  to  the  input  speeds.  However,  BRT- 
PCM  may  have  a speed  mesh,  including  upper  and  lower  limits, 
which  does  not  coincide  with  that  of  an  available  library.  In 
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this  case,  the  cross  sections  for  the  speed  points  in  the  BRT- 
PCM  speed  mesh  can  be  calculated  from  an  expression  obtained 
by  fitting  a curve  of  the  cross  section  versus  speed  data  from 
a known  library  (like  NITAWL)  and  then  extrapolating  or 
interpolating  as  necessary, 

a. 3)  Target  atom  directed  velocity 

The  BRT-PCM  input  file  has  a specific  data  entry  input 
for  the  target  atom  directed  velocity,  V„;  the  original  BRT-I 
and  BRT-PC  have  no  such  input  quantity, 
a. 4)  Problem  type  selection 

BRT-PCM  includes  the  options  of  either  no  directed 
velocity  or  directed  velocity  of  target  atoms  (completely 
entrained  neutrons) ; the  user  must  specify  the  selected 
problem  type  in  the  input  file, 
b)  Main  program. 

BRT-PCM  has  one  new  subroutine,  DIRVE,  that  is  inserted 
after  the  solution  of  the  neutron  integral  transport  equation 
(before  the  subroutine  EDIT) . This  subroutine  translates  the 
obtained  results  from  the  moving  frame  back  to  the  static 
frame.  The  main  program  also  has  to  read  the  new  variables 
from  the  modified  input  file  and  present  a way  to  control  the 
operational  sequence  according  to  the  specific  case  selected 
for  analysis  (see  Item  a. 4)  above). 
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D.3  Comparison  between  the  Results  from  BRT-PCM  and  BRT-I 

As  explained  above,  a few  modifications  were  made  in  the 
original  BRT-I  computer  code  to  create  the  personal  computer 
version,  BRT-PCM.  To  determine  how  the  modifications 
introduced  affected  the  results  obtained  with  BRT-PCM,  both 
BRT-I  and  BRT-PCM  codes  were  run  with  the  same  input  data. 
Since  BRT-I  is  not  capable  of  including  the  directed  velocity 
of  the  medium,  both  runs  were  performed  with  equal  to  zero. 

The  input  data  are  for  3 0 space  points  and  3 0 speed 
points,  with  all  the  isotopes  at  1000  °K  the  system  modeled  is 
similar  in  geometry  to  the  one  presented  in  Appendix  B, 
Section  B.l.  The  output  results  for  the  selected  key  variables 
are  presented  in  the  table  below. 


Table  D.l  - Results  of  BRT-PCM  and  BRT-I. 


VARIABLE 

REGION 

BRT-PCM 

BRT- 

-I 

ERROR  (%) 

E®  (cm’'') 

C 

1.168  E-2 

1.176 

E-2 

- 0.67 

S®  (cm*'') 

M 

5.509  E-4 

5.533 

E-4 

- 0.43 

(cm*') 

C 

9.682  E-3 

9.752 

E-3 

- 0.72 

(cm*'') 

C 

2.353  E-2 

2.370 

E-2 

- 0.72 

E®  (cm*') 

c 

3.619  E-2 

3.633 

E-3 

- 0.38 

E®  (cm*'') 

M 

8.153  E-1 

8.160 

E-1 

- 0.85 

1/v 

C 

4.335  E-1 

4 .360 

E-1 

- 0.58 

1/v 

M 

4.833  E-1 

4.855 

E-1 

- 0.45 
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In  the  table  above  all  the  cross  sections  are  region  and 
energy  averaged  values  and, 

C:  core  region, 

M:  moderator  region, 

1/v:  inverse  of  the  average  speed,  in  units  of  2200  cm/s. 

From  Table  D.l  it  can  be  verified  that  the  difference 
(error  involved)  between  the  results  obtained  from  BRT-PCM  and 
the  results  obtained  from  BRT-I  is  always  less  than  1%  and 
demonstrates  the  adequacy  of  using  BRT-PCM. 

D.4  Comparison  between  the  Results  from  BRT-PCM  for 
Different  Speed  and  Space  Structures 

The  original  BRT-I  is  usually  run  with  speed  and  space 
mesh  structures  consisting  of  30  points  each.  The  BRT-PCM 
version  is  also  capable  of  handling  the  same  mesh  structures 
but  this  requires  a large  amount  of  computer  time  on  the 
Personal  Computer  and,  since  the  present  work  requires  many 
runs  to  examine  many  different  cases,  it  was  decided  to  run 
BRT-PCM  using  speed  and  space  mesh  structures  with  only  20 
points  each.  The  results  obtained  from  this  coarse  structure 
are  not  expected  to  be  as  accurate  as  the  results  obtained 
with  a 3 0 point  structure  and  the  errors  involved  for  a 
typical  problem  are  shown  in  Table  D.2.  The  input  data  related 
to  isotopes,  number  densities  and  temperature  (2000  °K)  are 
exactly  the  same  for  both  cases.  The  system  modeled  in  both 
cases  is  similar  in  geometry  to  the  one  presented  in  Appendix 
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B,  Section  B.l;  the  only  difference  is  in  the  number  of  points 
for  the  space  and  speed  structure.  Both  cases  were  run  within 
the  same  speed  range,  from  0.3174  to  8.5744  (in  units  of  2200 
m/s)  , with  the  speed  points  equally  distributed  according  to 
the  number  of  speed  points  for  each  case  (20  and  3 0 speed 
points) . The  space  point  distribution  for  each  case  follows: 

- 20-point  structure: 

- fuel  gas  region  (core) : points  1-7 

- moderator  region:  points  8-20 

- 30 -point  structure: 

- fuel  gas  region  (core) : points  1-10 

- moderator  region:  points  11-30 


Table  D.2  - Results  of  BRT-PCM  for  20-Point  and  30-Point 
Structures . 


VARIABLE 

REGION 

20  POINTS 

30  POINTS 

ERROR  (%) 

S®  (cm’’') 

C 

8.093  E-3 

7.463  E-3 

8.44 

E®  (cm''') 

M 

3.869  E-4 

3.773  E-4 

2.54 

(cm"') 

C 

6.590  E-3 

6.081  E-3 

8.36 

vE^  (cm''') 

C 

1.601  E-2 

1.478  E-2 

8.36 

E®  (cm'') 

c 

3.282  E-2 

3.155  E-2 

4.01 

E®  (cm'') 

M 

8.066  E-1 

8.029  E-1 

0.46 

1/v 

c 

3.193  E-1 

3.013  E-1 

5.96 

1/v 

M 

3.359  E-1 

3.275  E-1 

2.58 
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Table  D.2  shows  the  results  for  the  fuel  gas  (core)  and 
moderator  region.  The  data  from  Table  D.2  show  for  the 
averaged  neutron  scattering  cross  section  approximately  0.4% 
error  in  the  moderator  region  and  approximately  4.0%  error  in 
the  core  region  and  for  the  averaged  macroscopic  neutron 
absorption  cross  section  approximately  2.5%  error  in  the 
moderator  region  and  approximately  8.4%  error  in  the  core 
region. 

To  determine  which  contribution  to  the  errors  presented 
in  Table  D.2  is  due  to  the  course  speed  mesh  and  which  is  due 
to  the  course  space  mesh,  one  case  was  run  with  a 3 0-point 
speed  structure  and  a 20-point  space  structure.  The  results 
are  presented  in  Table  D.3 


Table  D.3  - Results  of  BRT-PCM  for  20  and  30  Space  Point 
Structures  (30  Speed  Points) . 


VARIABLE 

REGION 

20  POINTS 

30  POINTS 

ERROR  (%) 

S®  (cm'’’) 

C 

7.793  E-3 

7.463  E-3 

4.42 

S®  (cm''') 

M 

3.823  E-4 

3.773  E-4 

1.32 

(cm'') 

C 

6.347  E-3 

6.081  E-3 

4.37 

vS^  (cm''') 

c 

1.542  E-2 

1.478  E-2 

4.33 

S®  (cm''') 

c 

3.220  E-2 

3.155  E-2 

2.06 

2®  (cm'') 

M 

8.048  E-1 

8.029  E-1 

0.23 

1/v 

c 

3.104  E-1 

3.013  E-1 

3.02 

1/v 

M 

3.319  E-1 

3.275  E-1 

1.34 
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To  compare  the  results,  the  geometry  was  kept  the  same  as  for 
the  previous  cases,  the  speed  distribution  is  the  same  used 
for  the  above  30-point  structure  case  and  the  space 
distribution  is  the  one  used  for  the  above  20-point  structure 
case . 

Comparing  the  errors  presented  in  Tables  D.2  and  D.3,  it 
can  be  observed  that  the  errors  involved  for  the  case  when 
only  the  space  point  structure  is  modified  are  approximately 
50  % smaller  than  the  errors  for  the  case  when  both  the  speed 
point  and  space  point  structures  are  modified.  If  it  is 
assumed  that  the  space  and  speed  effects  are  independent  of 
each  other,  it  can  be  concluded  that  each  contributes 
approximately  50  % to  the  total  errors  in  Table  D.2. 


D.5  Source  Listing  of  BRT-PCM 


Q ******  battelle  revised  thermos  ****** 

C 

C ******  UNIVERSITY  OF  FLORIDA  VERSION  FOR  IBM  AT  ****** 

C 

C MAIN  PROGRAM 

C 

C 

C ********************************************************* 

C PC  MODIFIED  VERSION  OF  BRT-I  TO  INCLUDE  THE  TARGET  ATOM 

C DIRECTED  VELOCITY 

C ********************************************************* 

PROGRAM  BRT 
C 


NANIS0(2) 

HOLID(30) 

LINE 
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CHARACTER* 4 
CHARACTER* 18 
CHARACTER* 7 2 


C 


C 

C 


C 

100 


DIMENSION  T(30,30,30),  P(30,30,8), 

+ XTM(30,30),  XAM(30,30),  S(30,30),  Q(30,30), 

+ F(30,30),  Fl(30,30),  CONCTA ( 30 , 8 ) , XAT(30,30), 

+ XFT(30,30),  XST(30,30),  XS1(30,30),  RI(30), 

+ RO(30),  VNUT(30),  WSTBA(30),  DV(30),  V(30), 

+ RAD(30),  MTBL(30),  VOL(30),  WSTBB(30), 

+ NRTBL (30) 

DATA  NANISO  / 4H  NO,  4H  YES  / 


WRITE  (*,*)  'BRTPCM' 

NTIN  = 01 
NTOUT  =02 
NTEMP  =03 

OPEN (NTIN , FILE= ' INPUT ' , STATUS= ' OLD ' ) 

OPEN (NTOUT , FILE= ' BRT . OUT ' , STATUS= ' UNKNOWN ' ) 
OPEN(NTEMP,  FILE='TMP',  STATUS= ' UNKNOWN • ) 


READ  (NTIN, ' (A) ' , END=999)  LINE 
WRITE  (NTOUT, ' (A) ' ) LINE 


C 

READ  (NTIN, 295, ERR=999)  IDENTA,  IBY,  VM,  ICASE 
C 

C INPUT  OF  DIRECTED  VELOCITY,  VM,  AND  ICASE.  ICASE  CAN  BE 

C 0,  IF  VM=0,  0 FOR  UNCOLLIDED  NEUTRONS  AND  -1  FOR 

C ENTRAINED  NEUTRONS 
C 

WRITE  (NTOUT,  300) 

WRITE  (NTOUT,  315)  VM,  ICASE 
IF  (IBY)  130,  140,  130 
130  WRITE  (NTOUT, 310) IDENTA,  IDENT 
IDENT=IDENTA 
IF  (IBY)  270,270,180 
C 

140  DO  150  1=1,  30 

T(I,  1,  1)  = 0.0 
150  CONTINUE 
C 


IDENT  = IDENTA 

READ  (NTIN, 290)  NX,  IX,  MX,  ISOX,  ISOXE,  ICX,  ICXE,  NANI, 
+ NRSMR,  NBUCK 

C READ  (NTIN, 295)  VM,  CASE 
JYES  = 1 

IF  (NBUCK. NE.O)  JYES=2 
lYES  = NANI  + 1 

WRITE  (NTOUT, 330) IDENT,  NX,  IX,  MX,  ISOX,  ISOXE,  ICX, 

+ ICXE,  NANISO(IYES) , NANISO(JYES) 

IF  (NBUCK. NE.O)  READ  (NTIN, 320)  BF,DZ,DY 
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IF  (ISOX+ISOXE+ICX+ICXE)  160,  160,  170 
C 


160 

C 

STOP  ' 

ERROR  100' 

170 

ISOXT 

= ISOX  + ICX 

ISOXET 

= ISOXE  + ICXE 

ISOT 

= ISOXT  + ISOXET 

IF  (ISOT.GT.30)  STOP  'ERROR  103' 

READ  (NTIN,340)  (MTBL(N) , N=l,  NX) 

C ********************************************************* 


WRITE  (*,*)  'TAPISO' 

CALL  TAPISO  (DV,  V,  CONCTA,  ISOX,  ISOXE,  IX,  ISOXT,  MX, 
+ NTIN) 

WRITE (*,*)  'CARISO' 

CALL  CARISO  (XTM,  XAM,  IX,  V,  DV,  P,  CONCTA,  ISOX,  XAT, 
+ XFT,  XST,  VNUT,  HOLID,  WSTBA,  WSTBB,  ISOXT, 

+ MX,  NTIN,  BF,  NTOUT,  DY,  DZ , ICX,  NANI,  NX, 

+ ICXE,  MTBL,  NBUCK,  IDENT,  ISOT) 

WRITE  (*,*)  'SOURCE' 

CALL  SOURCE  (NRSMR,  IX,  NX,  MTBL,  IDENT,  DV,  V, 

+ NTIN,  NTOUT,  S) 

C ********************************************************* 
180  CONTINUE 

READ  ( NTIN , 3 5 0 ) LEANT , NXA , NGEOM , NCUR , NPSMR , MICROP , NAX , CONA 
ISAME=LEAKT 

C NGEOM=l  REGULAR  CYLINDRICAL  GEOMETRY  WITH  ARGONNE 
C CORRECTIONS 

C NGEOM=2  FAST  COSINE  CURRENTS  CYLINDRICAL  GEOMETRY  FROM 
C SRL 

C NGEOM=3  REGULAR  SLAB  GEOMETRY  WITH  ARGONNE  CORRECTIONS 
C NGEOM=4  FAST  COSINE  CURRENTS  SLAB  GEOMETRY  FROM  SRL 
GOTO  (190,200,210,220),  NGEOM 
190  WRITE  (*,*)  'GEOM' 

CALL  GEOM (NTIN,  NTOUT,  CONA,  IX,  NX,  IDENT,  XTM,  RI , RO, 
+ RAD,  VOL,  T,  Q,  NAX) 

C LABELS  210,  220,  230  WERE  CHANGED 
200  GO  TO  230 

WRITE  (*,*)  'GEORl' 

C200  CALL  GEORl 
210  GO  TO  230 

WRITE  (*,*)  'GEOX' 

C210  CALL  GEOX 
220  GO  TO  230 

WRITE  (*,*)  'GEOXl' 

C220  CALL  GEOXl 
C230  CALL  ETIMEI  (NQ2) 

230  IF  (ISAME)  240,240,250 
240  LEAKT=1 

GO  TO  260 
250  LEAKT=-1 
260  WRITE  (*,*)  'ITER' 

CALL  ITER  (NPSMR,  P,  T,  V,  RAD,  NTOUT,  DV,  IX,  NX,  XAM, 
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+ MTBL,  LEAKT,  VOL,  IDENT,  NGEOM,  ISAME,  NCUR, 

+ RI,  RO,  NANI,  XTM,  S,  Q,  F,  FI) 

C CALL  ETIMEI  (NQ3) 

IF  (ICASE)  261,270,262 

261  WRITE  (*,*)  'DIRVE' 

C CALL  DIRVE  (VM,  V,  DV,  XST,  XAT,  XFT,  XAM,  XTM,  XSl,  IX, 
C + NX,  ISOX,  MX,  P,  F,  FI) 

C 

GOTO  270 

262  WRITE  (*,*)  'DIRVU' 

CALL  DIRVU  (VM,  V,  DV,  XST,  XAT,  XFT,  XAM,  XTM,  XSl,  IX, 
+ NX,  ISOX,  MX,  P,  F,  FI) 

C 

270  WRITE  (*,*)  'EDIT' 

CALL  EDIT  (NCUR,  WSTBA,  WSTBB,  ISOX,  VNUT,  V,  RAD,  NTOUT, 
+ DV,  IX,  IBY,  XAT,  MTBL,  LEAKT,  VOL,  XFT, 

+ IDENT,  CONCTA,  NTIN,  NX,  ISOT,  XST,  RI , RO,  F, 

+ FI,  HOLID,  NRTBL,  XSl) 

WRITE  (*,*)  'SMR' 

IF  (NPSMR.NE.O)  CALL  SMR  (NRTBL,  ISOX,  F,  MICROP,  VNUT, 
+ NPSMR,  V,  NTOUT,  XSl,  IX,  NANI,  NX,  XAM, 

+ XAT,  MTBL,  VOL,  XFT,  IDENT,  CONCTA,  XTM, 

+ P) 

C CALL  ETIMEI  (NQ4) 

C CALL  SUMARY 

GO  TO  100 
C 

290  FORMAT  (1015) 

295  FORMAT  (215 , ElO . 5 , 15) 

300  FORMAT(/,  52X,  '****************************'/ 

+ , 52X,  '*  BATTELLE  REVISED  THERMOS  *'/ 

+ , 52X,  '*  U OF  F IBM  AT  VERSION  *'/ 

+ , 52X,  '****************************') 

310  FORMAT  (6H0IDENTI6,  15H  SAME  XSECT  ASI6) 

315  FORMAT  (//,'VM  = ' , ElO . 5 , lOX , ' CASE  = ',15) 

320  FORMAT  (7E10.5) 

330  FORMAT  ('IDENT  = ',16/ 'SPACE  POINTS  = ',13/ 'GROUPS  = ', 
+ 13/ 'MIXTURES  =' ,13/ 'LIBRARY  ISOTOPES  IN  CELL  = ',13/ 

+ 'LIBRARY  ISOTOPES  IN  EDIT  = ',13/  'ADDED  ISOTOPES  IN 
+ CELL  = ',I3/'ADDED  ISOTOPES  IN  EDIT  = ', 13/ ' ANISOTROPY 
+ CORRECTION  TO  SCATTERING  KERNEL  = ',  A4/ ' TRANSVERSE 
+ BUCKLING  CORRECTION  = ' , A4 ) 


340 

FORMAT 

(5011) 

350 

FORMAT 

(715, ElO. 5) 

999 

STOP 

END 
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D.6  Source  Listing  of  DIRVE 


SUBROUTINE  DIRVE  (VM,  V,  DV,  XST,  XAT,  XFT,  XAM,  XTM, 

+ XSl,  IX,  NX,  ISOX,  MX,  P,  F,  FI) 

DIMENSION  V (30) , VE(30),  DV(30),  C(30),  XST(30,30), 

+ XAT(30,30),  XFT(30,30) ,XAM(30,30) , E(30), 

+ XTM(30,30),  XS1(30,30),  F(30,30),  Fl(30,30), 

+ P(30,30,8) 

C 

C ********************************************************* 

C THIS  SUBROUTINE  CHANGES  THE  VARIABLES  FROM  THE  MOVING 
C FRAME  TO  THE  STATIC  FRAME  FOR  THE  CASE  OF  COMPLETELY 
C ENTRAINED  NEUTRONS 

C ********************************************************* 


c 

c 

c 

c 


100 

c 

c 


110 

c 

c 


120 

c 

c 

130 

C 

C 

140 

150 


NTOUT=  02 

OPENING  A FILE  TO  STORE  THE  XSl  CONVERTED  TO  THE  MOVING 
FRAME 

OPEN  (10,  FILE  = 'CLD',  STATUS= ' UNKNOWN ' ) 

CALL  OF  XSl  FROM  FILE  09  TO  BE  MODIFIED 
REWIND  09 
DO  100  J=l,ISOX 

READ  (09,440)  (XSl ( J , I ) , 1=1 , IX) 


CROSS  SECTIONS  IN 
DO  110  J=l,ISOX 
WRITE  (NTOUT,550) 
WRITE  (NTOUT,500) 
WRITE  (NTOUT,560) 
WRITE  (NTOUT,520) 
WRITE  (NTOUT,560) 
WRITE  (NTOUT,540) 
WRITE  (NTOUT,560) 
CONTINUE 


MOVING  FRAME 
J 

(XST(J,I) ,1=1, IX) 
(XAT(J,I) ,1=1, IX) 
(XFT(J,I) ,1=1, IX) 


CONVERTION  OF  NEUTRON  SPEEDS  TO  THE  STATIC  FRAME 

DO  120  1=1, IX 

VE(I)=V(I) 

CONTINUE 
DO  160  1=1, IX 

SPEED  IN  THE  MOVING  FRAME  SMALLER  THAN  THE  DIRECTED 
VELOCITY  OF  THE  MEDIUM 
IF  (V(I)-VM)  130,140,140 
V(I)=VM+(V(I)**2/(3.*VM) ) 

GOTO  150 

SPEED  IN  THE  MOVING  FRAME  GREATER  THAN  THE  DIRECTED 
VELOCITY  OF  THE  MEDIUM 
V(I)=V(I)+(VM**2/(3.*V(I) ) ) 

C(I)=VE(I)/V(I) 
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E(I)=((V(I)*2200) **2)/1.916E08 
160  CONTINUE 
C 

C CHANGE  OF  SPEED  MESH 
D=0.01 

170  DV(1)=D*V(1) 

DO  190  1=2, IX 

DV(I)=2*(V(I)-V(I-1) )-DV(I-l) 

IF  (DV(I) ) 180, 180, 190 
180  D=D+0.01 

GO  TO  170 
190  CONTINUE 
C 

C CONVERTION  OF  CROSS  SECTIONS  TO  THE  STATIC  FRAME 
DO  200  J=l,ISOX 
DO  200  1=1, IX 
XAT(J,I)=C(I) *XAT(J,I) 

XFT(J,I)=C(I) *XFT(J,I) 

XST(J,I)=C(I) *XST(J,I) 

XTM(J,I)=C(I) *XTM(J,I) 

XAM(J,I)=C(I) *XAM(J,I) 

XS1(J,I)=C(I) *XS1(J,I) 

200  CONTINUE 
C 

DO  220  J=l,ISOX 

220  WRITE  (10,440)  (XSl ( J , I ) , 1=1 , IX) 

C OUTPUT  OF  N(N,I),  V(I) 

WRITE  (NTOUT,400) 

NC=1 

ND=10 

230  IF  (NX-ND)  240,250,250 
240  ND=NX 

250  WRITE  (NTOUT,420)  (K,K=NC,ND) 

WRITE  (NTOUT,430) 

DO  260  1=1, IX 

260  WRITE  (NTOUT,410)  I,  (F (K, I) , K=NC, ND) 

IF  (NX-ND)  280,280,270 
270  NC=NC+10 
ND=ND+10 
GOTO  230 
280  CONTINUE 

C MODIFIED  VALUE  OF  V(I)  IN  STATIC  FRAME 
WRITE  (NTOUT,450) 

WRITE  (NTOUT,470)  (V (I) , 1=1 , IX) 

C MODIFIED  VALUE  OF  DV(I)  IN  STATIC  FRAME 
WRITE  (NTOUT,460) 

WRITE  (NTOUT,480)  (DV(I) , 1=1, IX) 

C ENERGY  POINTS  IN  THE  STATIC  FRAME 
WRITE  (NTOUT,570) 

WRITE  (NTOUT,580)  (E (I) , 1=1, IX) 

C 

C MODIFIED  VALUE  OF  CROSS  SECTIONS  IN  STATIC  FRAME 


c 

290 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 


201 


DO  290  J=l,ISOX 
WRITE  (NTOUT,550) 
WRITE  (NTOUT,490) 
WRITE  (NTOUT,560) 
WRITE  (NTOUT,510) 
WRITE  (NTOUT,560) 
WRITE  (NTOUT,530) 
WRITE  (NTOUT,560) 


J 

(XST(J,I) ,1=1, IX) 
(XAT(J,I) ,1=1, IX) 
(XFT(J,I) ,1=1, IX) 


CONTINUE 

RETURN 

FORMAT  (/,50X,  ' N IN  STATIC  FRAME') 

FORMAT  (I5,2X,10(2X,E10.5) ) 

FORMAT  (1H03X,1HI6X,10(2HN(I2,3H,I)5X) ) 

FORMAT  (3X,4H — ,10(12H  )) 

FORMAT  (7E10.5) 

FORMAT  (/, 5 OX, 'NEUTRON  VELOCITY  IN  STATIC  FRAME') 
FORMAT  (/,50X,'DV  IN  STATIC  FRAME') 

FORMAT  (10(2HV=(E10.5) ,2X) ) 

FORMAT  (10(3HDV=(E10.5) ,2X) ) 

FORMAT  (/,50X,  'SIGMA  S IN  STATIC  FRAME') 

FORMAT  (/,50X,  'SIGMA  S IN  MOVING  FRAME') 

FORMAT  (/,50X,  'SIGMA  A IN  STATIC  FRAME') 

FORMAT  (/,50X,  'SIGMA  A IN  MOVING  FRAME') 

FORMAT  (/,50X,  'SIGMA  F IN  STATIC  FRAME') 

FORMAT  (/,50X,  'SIGMA  F IN  MOVING  FRAME') 

FORMAT  (/,55X, 'ISOTOPE=' ,13) 

FORMAT  (10 (2X,E10.5) ) 

FORMAT  (/,50X, 'ENERGY  IN  STATIC  FRAME') 

FORMAT  (10(2HE=(E10.5) ,2X) ) 

END 
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